笔记

第1章 波函数与Schrödinger方程

1.1 波函数及其统计诠释

波粒二象性的分析

实物粒子的波粒二象性(particle-wave duality)
物质波

de Broglie(1923)提出实物粒子(静质量 m0m\ne0 的粒子)具有波粒二象性的假设,即与动量为 pp 和能量为 EE 的粒子相应的物质波(matter wave)的波长 λ\lambda 和频率 ν\nu

λ=hp,ν=Eh \lambda = \frac{h}{p}, \kern 12pt \nu = \frac{E}{h}

并称之为物质波

证明实物粒子波动性的实验
  • Davisson-Germer电子衍射实验(1927)

  • G.P.Thomson电子衍射实验(1927)

  • Jonsson实验(1961)

  • Arndt的 C60\text{C}_{60} 分子干涉实验(1999)

实物粒子波粒二象性的理解

以电子为例,它是粒子和波动两重性矛盾的统一,但这个波不再是经典概念下的波,粒子也不再是经典概念中的粒子

粒子性

指微观粒子与物质相互作用时的“颗粒性”(corpuscularity)或“原子性”(atomicity),具有集中的能量EE与动量p \vec{p}\

但与经典的粒子不同,微观粒子没有确定的轨道,而应采用“概率”的概念。

波动性

指微观粒子在空间传播时的“相干(coherent)叠加性”,有“干涉”、“衍射”、“偏振”等现象,具有波长λ\lambda和波矢k \vec{k}\

但与经典的波不同,没有某种实际物理量(如质点的位移、电场、磁场)的波动分布.

状态的描述

经典力学

经典力学中的质点由 r\vec{r},p(v)\vec{p}(\vec{v}) 描写:

  • 每一时刻该二量具有完全确定的值,且随时间连续变化;

  • 质点的其他力学量(Ek,V,LE_k,V,\vec{L},) 都可表示为 r\vec{r},p\vec{p}, 的函数( r\vec{r},p\vec{p}, 完全决定了质点的性质);

  • 质点状态的变化遵从牛顿定律:若已知 r0\vec{r}_0p0\vec{p}_0 ,则任时刻的 r(t)\vec{r}(t)p(t)\vec{p}(t) 唯一确定

{v(t)=0tFmdt+v0p(t)=0tFdt+p0r(t)=0tv(t)dt+r0 \left{\begin{matrix} \vec{v}(t) = \int_0^t \frac{\vec{F}}{m}\mathrm{d}t + \vec{v}_0 \\ \vec{p}(t) = \int_0^t \vec{F}\mathrm{d}t + \vec{p}_0 \\ \vec{r}(t) = \int_0^t \vec{v}(t)\mathrm{d}t + \vec{r}_0 \end{matrix}\right.

  • r(t)\vec{r}(t) 描写粒子的运动轨道。
量子力学

量子力学中微观粒子的状态由波函数(wave function)描写:

  • 微观粒子不可能同时具有确定的 r\vec{r},p\vec{p}, ,也就是没有确定的轨道

  • 对于一般状态的微观粒子,应该用一般的时间和空间的复函数 ψ(r,t)\psi(\vec{r},t) 来描写,它称为波函数(亦称态矢量)。波函数是在空间的一个分布(在给定时间 tt ,它是坐标的函数 ψ(r)\psi(\vec{r}) ),是微观粒子波粒二象性的表现。

  • 波函数 ψ(r,t)\psi(\vec{r},t) 可以决定微观粒子的一切力学量和行为;

  • 波函数 ψ(r,t)\psi(\vec{r},t) 的变化遵从Schrödinger方程

波函数的统计诠释

概率波

M.Born(1926)提出的概率波把微观粒子的“原子性”与波的“相干叠加性”统一在了一起。

量子力学假定之一:一个微观粒子的状态总可以用一个波函数 ψ(r,t)\psi(\vec{r},t) 来完全描述,波函数是粒子坐标和时间的复函数,模平方 ψ(r,t)2|\psi(\vec{r},t)|^2 代表粒子空间分布的概率密度。波函数本身称为概率波幅(probability amplitude).

ψ(r,t)2ΔxΔyΔz|\psi(\vec{r},t)|^2\Delta x\Delta y\Delta z 表示在 r\vec{r}, 点处的体积 ΔxΔyΔz\Delta x\Delta y\Delta z 中找到粒子的概率.

波函数的归一

对于概率分布来说,重要的是相对概率分布,所以将波函数乘以一个常数,它仍然描写量子体系的同一个状态。即对于任意非零常数 CC ,波函数 ψ(r,t)\psi(\vec{r},t)Cψ(r,t)C\psi(\vec{r},t) 描述的相对概率完全相同,这表明波函数有一个常数因子不确定性

根据波函数的统计诠释,很自然要求微观粒子(不产生,不湮没)在空间各点的概率之总和为 11 ,即要求波函数 ψ(r,t)\psi(\vec{r},t) 满足归一化条件

()ψ(r,t)2d3r=1 \int_{(全)} \left|\psi(\vec{r},t)\right|^2 \mathrm{d}^3r = 1

一般的,若波函数 ψ(r,t)\psi(\vec{r},t) 满足平方可积条件,即存在有限正常数 AA ,使得

()ψ(r,t)2d3r=A \int_{(全)} \left|\psi(\vec{r},t)\right|^2 \mathrm{d}^3r = A

则有

()1Aψ(r,t)2d3r=1 \int_{(全)} \left|\frac{1}{\sqrt{A}}\psi(\vec{r},t)\right|^2 \mathrm{d}^3r = 1

应当注意,即使加上归一化条件,波函数仍有相位(phase)不定性,即假设 ψ(r,t)\psi(\vec{r},t) 是归一化的波函数,则 eiδψ(r,t)\mathrm{e}^{\mathrm{i}\delta}\psi(\vec{r},t) 也是归一化的,且二者描述的是同一个概率波。

对于某些理想(非物理)的情况,波函数是不能归一的,例如平面波(自由粒子的波函数):ψ(r,t)=Aei(prEt)\psi(\vec{r},t) = A\mathrm{e}^{\frac{\mathrm{i}}{\hbar}\left(\vec{p}\cdot\vec{r}-Et\right)} ,其有以下两种非常规的归一化方式。

箱归一化

平面波是理想模型,实际上应该用“波包”来描述自由粒子,即粒子分布在有限空间,例如分布在 L2xL2-\frac{L}{2} \le x \le \frac{L}{2} 内,这时的波函数

ψp(x)=1Leipx \psi_p(x) = \frac{1}{\sqrt{L}} \mathrm{e}^{\frac{\mathrm{i}}{\hbar}px}

其称为箱归一化的平面波,满足

L2L2ψp(x)2dx=1 \int_{-\frac{L}{2}}^{\frac{L}{2}} \left|\psi_p(x)\right|^2 \mathrm{d}x = 1

对于三维的情况,用 Ω\Omega 表示自由粒子分布的体积,则波函数

ψp(r)=1Ωeipx \psi_{\vec{p}}(\vec{r}) = \frac{1}{\sqrt{\Omega}} \mathrm{e}^{\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{x}}

满足

Ωψp(r)2d3r=1 \int_{\Omega} \left|\psi_{\vec{p}}(\vec{r})\right|^2 \mathrm{d}^3r = 1

δ\delta 函数“归一化”(规格化)

为处理连续谱本征函数的“归一化”,Dirac引进了 δ\delta 函数,其定义为:

δ(x)={0,x0,x=0 \delta(x) = \begin{cases} 0, & x \ne 0 \ \infty, & x = 0 \end{cases}

εεδ(x)dx=+δ(x)dx=1(ε>0) \int_{-\varepsilon}^{\varepsilon} \delta(x) \mathrm{d}x = \int_{-\infty}^{+\infty} \delta(x) \mathrm{d}x = 1 \kern 12pt \left( \varepsilon > 0 \right)

或等价的表示为:对于在 x=x0x = x_0 邻域连续的任意函数 f(x)f(x)

+f(x)δ(xx0)dx=f(x0) \int_{-\infty}^{+\infty} f(x)\delta(x-x_0) \mathrm{d}x = f(x_0)

δ\delta 函数有如下性质:

δ(x)=12π+e±ikxdk \delta(x) = \frac{1}{2\pi} \int_{-\infty}^{+\infty} \mathrm{e}^{\pm\mathrm{i}kx} \mathrm{d}k

δ(x)=12π+e±ipxdp \delta(x) = \frac{1}{2\pi\hbar} \int_{-\infty}^{+\infty} \mathrm{e}^{\pm\frac{\mathrm{i}}{\hbar}px} \mathrm{d}p

δ(x)=δ(x) \delta(-x) = \delta(x)

δ(ax)=1aδ(x) \delta(ax) = \frac{1}{|a|} \delta(x)

xδ(xa)=aδ(xa) x\delta(x-a) = a\delta(x-a)

δ\delta 规格化的平面波为

ψp(x)=12πeipx \psi_p(x) = \frac{1}{\sqrt{2\pi\hbar}} \mathrm{e}^{\frac{\mathrm{i}}{\hbar}px}

其满足 δ\delta 函数规格化条件

+ψp(x)ψp(x)dx=12π+ei(pp)xdx=δ(pp) \int_{-\infty}^{+\infty} \psi^*{p'}(x), \psi_p(x), \mathrm{d}x = \frac{1}{2\pi\hbar} \int{-\infty}^{+\infty} e^{\frac{\mathrm{i}}{\hbar}(p-p')x} \mathrm{d}x = \delta(p-p')

对于三维的情况,波函数

ψp(r)=1(2π)32eipr \psi_{\vec{p}}(\vec{r}) = \frac{1}{\left(2\pi\hbar\right)^{\frac{3}{2}}} \mathrm{e}^{\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{r}}

其满足 δ\delta 函数规格化条件

+ψp (r)ψp(r)d3r=δ(pp )=δ(pxpx)δ(pypy)δ(pzpz) \int_{-\infty}^{+\infty} \psi^*{\vec{p}\ '}(\vec{r}), \psi{\vec{p}}(\vec{r}), \mathrm{d}^3r = \delta(\vec{p}-\vec{p}\ ') = \delta(p_x-p'_x) \delta(p_y-p'_y) \delta(p_z-p'_z)

多粒子波函数

对于 NN 个粒子组成的体系,它的波函数表示为

ψ(r1,r2,,rN) \psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N)

其中 r1(x1,y1,z1),r2(x2,y2,z2),,rN(xN,yN,zN)\vec{r}_1(x_1,y_1,z_1),\vec{r}_2(x_2,y_2,z_2),\cdots,\vec{r}_N(x_N,y_N,z_N) 分别表示各粒子的空间坐标,此时该波函数描述的是抽象的 3N3N 维位形空间(configuration space)中的概率波,

ψ(r1,r2,,rN)2d3r1d3r2d3rN \left|\psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N)\right|^2 \mathrm{d}^3r_1 \mathrm{d}^3r_2 \cdots \mathrm{d}^3r_N

表示粒子 11 出现在 (r1,r1+dr1)(\vec{r}_1,\vec{r}_1+\mathrm{d}\vec{r}_1) 中,同时粒子 22 出现在 (r2,r2+dr2)(\vec{r}_2,\vec{r}_2+\mathrm{d}\vec{r}_2) 中……同时粒子 NN 出现在 (rN,rN+drN)(\vec{r}_N,\vec{r}_N+\mathrm{d}\vec{r}_N) 中的概率,对应的归一化条件表示为

()ψ(r1,r2,,rN)2d3r1d3r2d3rN=1 \int_{(全)} \left|\psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N)\right|^2 \mathrm{d}^3r_1 \mathrm{d}^3r_2 \cdots \mathrm{d}^3r_N = 1

概率诠释对波函数的要求
平方可积

一般来说 ψ(r)\psi(\vec{r}) 应处处取为有限值,但在平方可积的条件下:

()ψ(r)2d3r=有限正常数 \int_{(全)} \left|\psi(\vec{r})\right|^2 \mathrm{d}^3r = 有限正常数

可以存在有限个孤立奇点。

归一化条件

一个真实的波函数需要满足归一化条件

()ψ(r)2d3r=1 \int_{(全)} \left|\psi(\vec{r})\right|^2 \mathrm{d}^3r = 1

但在量子力学中并不排除使用某些不能归一化的理想的波函数,如平面波 ψ(r)eipr\psi(\vec{r}) \sim e^{\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{r}}δ\delta 波包 ψ(r)δ(r)\psi(\vec{r}) \sim \delta(\vec{r})

单值性

要求 ψ(r)2|\psi(\vec{r})|^2 单值,即粒子的概率分布是确定的,但不能由此要求 ψ(r)\psi(\vec{r}) 单值。

连续性

波函数 ψ(r)\psi(\vec{r}) 及其各阶微商的连续性,要根据体系所处的势场 V(r)V(\vec{r}) 的性质来分析。

内积

波函数 ψ\psiϕ\phi 的内积(inner product)定义为

(ψ,ϕ)=+ψ(x)ϕ(x)dx \left(\psi,\phi\right) = \int_{-\infty}^{+\infty} \psi^*(x)\phi(x) \mathrm{d}x

内积是态矢空间中两个态矢量的“点乘”,是一个复数,其有以下性质:

(ψ,ψ)0 \left(\psi,\psi\right) \ge 0

(ψ,ϕ)=(ϕ,ψ)=(ϕ,ψ) \left(\psi,\phi\right) = \left(\phi,\psi\right)^* = \left(\phi^,\psi^\right)

(ψ,C1ϕ1+C2ϕ2)=C1(ψ,ϕ1)+C2(ψ,ϕ2) \left(\psi,C_1\phi_1+C_2\phi_2\right) = C_1\left(\psi^,\phi_1\right) + C_2\left(\psi^,\phi_2\right)

(C1ψ1+C2ψ2,ϕ)=C1(ψ1,ϕ)+C2(ψ2,ϕ) \left(C_1\psi_1+C_2\psi_2,\phi\right) = C_1^\left(\psi_1^,\phi\right) + C_2^\left(\psi_2^,\phi\right)

特别的,内积没有对称性,即一般

(ψ,ϕ)=(ϕ,ψ)(ϕ,ψ) \left(\psi,\phi\right) = \left(\phi,\psi\right)^* \ne \left(\phi,\psi\right)

(ψ,ϕ)=0\left(\psi,\phi\right) = 0 时,称 ψ\psiϕ\phi 正交

使用 ()dτ\int_{(全)} \mathrm{d}\tau 代表对体系的全部坐标空间进行积分,例如

对于一维粒子

()dτ=+dx \int_{(全)} \mathrm{d}\tau = \int_{-\infty}^{+\infty} \mathrm{d}x

对于三维粒子

()dτ=+++dxdydz \int_{(全)} \mathrm{d}\tau = \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} \mathrm{d}x \mathrm{d}y \mathrm{d}z

对于 NN 个粒子组成的体系

()dτ=++dx1dy1dz1dxNdyNdzN \int_{(全)} \mathrm{d}\tau = \int_{-\infty}^{+\infty} \cdots \int_{-\infty}^{+\infty} \mathrm{d}x_1 \mathrm{d}y_1 \mathrm{d}z_1 \cdots \mathrm{d}x_N \mathrm{d}y_N \mathrm{d}z_N

在内积的定义下,有

(ψ,ψ)=()dτψψ=()dτψ2 \left(\psi,\psi\right) = \int_{(全)} \mathrm{d}\tau \psi^*\psi = \int_{(全)} \mathrm{d}\tau |\psi|^2

这样就可以简单的表示归一化条件为

(ψ,ψ)=1 \left(\psi,\psi\right) = 1

坐标表象和动量表象上的波函数

ψ(r)2|\psi(\vec{r})|^2 表示粒子在坐标空间中的概率密度相似, φ(p)2|\varphi(\vec{p})|^2 表示粒子的动量分布的概率密度,(归一化后)粒子动量在 (p,p+dp)(\vec{p},\vec{p}+\mathrm{d}\vec{p}) 范围中概率为 φ(p)2d3p|\varphi(\vec{p})|^2 \mathrm{d}^3p

粒子的量子态,既可以用 ψ(r)\psi(\vec{r}) 描述,也可以用 φ(p)\varphi(\vec{p}) 来描述(还可以有其他的描述方式)。它们彼此间有确定的变换关系,彼此完全等价,描述的都是同一个量子态,只不过表象(representation)不同而已。称 ψ(r)\psi(\vec{r}) 是粒子态在坐标表象中的表示,而 φ(p)\varphi(\vec{p}) 则是同一个状态在动量表象中的表示。

波函数 ψ\psiφ\varphi 之间满足Fourier变换,在一维情形下

φ(p)=1(2π)12+ψ(x) eipx dx \varphi(p) = \frac{1}{(2\pi\hbar)^{\frac{1}{2}}} \int_{-\infty}^{+\infty} \psi(x)\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}px}\ \mathrm{d}x

ψ(x)=1(2π)12+φ(p) eipx dp \psi(x) = \frac{1}{(2\pi\hbar)^{\frac{1}{2}}} \int_{-\infty}^{+\infty} \varphi(p)\ \mathrm{e}^{\frac{\mathrm{i}}{\hbar}px}\ \mathrm{d}p

在三维情形下

φ(p)=1(2π)32+ψ(r) eipr d3r \varphi(\vec{p}) = \frac{1}{(2\pi\hbar)^{\frac{3}{2}}} \int_{-\infty}^{+\infty} \psi(\vec{r})\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{r}}\ \mathrm{d}^3r

ψ(r)=1(2π)32+φ(p) eipr d3p \psi(\vec{r}) = \frac{1}{(2\pi\hbar)^{\frac{3}{2}}} \int_{-\infty}^{+\infty} \varphi(\vec{p})\ \mathrm{e}^{\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{r}}\ \mathrm{d}^3p

可以推得,两种表象上的波函数的归一化是等价的,即

(ψ,ψ)=(φ,φ) \left(\psi,\psi\right) = \left(\varphi,\varphi\right)

(ψ,ψ)=1(φ,φ)=1 \left(\psi,\psi\right) = 1 \Longleftrightarrow \left(\varphi,\varphi\right) = 1

力学量用算符表达

算符代表对波函数的某种作用或运算。

平均值假定

粒子处于波函数 ψ(r)\psi(\vec{r}) 所描述的状态下,力学量(又叫可观测量)都有确定的概率分布,因而有确定的平均值(又叫期待值)。在任意状态 ψ\psi 上,对力学量 AA 进行足够多次的测量,所得结果的平均值为

Aˉ=(ψ,A^ψ)(ψ,ψ) \bar{A} = \frac{(\psi,\hat{A}\psi)}{(\psi,\psi)}

其中 A^\hat{A} 是力学量 AA 对应的算符,若波函数已归一化,则

Aˉ=(ψ,A^ψ) \bar{A} = (\psi,\hat{A}\psi)

坐标算符与动量算符

在波函数 ψ\psi 已归一化的条件下,位置 xx 的平均值为

xˉ=+ψ(r)2xd3r=+ψ(r)xψ(r)d3r \bar{x} = \int_{-\infty}^{+\infty} \left|\psi(\vec{r})\right|^2, x, \mathrm{d}^3r = \int_{-\infty}^{+\infty} \psi^*(\vec{r}), x, \psi(\vec{r}), \mathrm{d}^3r

可以得到坐标表象下的坐标算符为

x^=x \hat{x} = x

同理

y^=y,z^=z,r^=r \hat{y} = y, \kern 12pt \hat{z} = z, \kern 12pt \hat{\vec{r}} = \vec{r}

如果状态用动量表象波函数 φ(p)\varphi(\vec{p}) 来表示,则粒子动量的平均值为

pˉ=+φ(p)2pd3p=+φ(p)pφ(p)d3p \bar{\vec{p}} = \int_{-\infty}^{+\infty} \left|\varphi(\vec{p})\right|^2, \vec{p}, \mathrm{d}^3p = \int_{-\infty}^{+\infty} \varphi^*(\vec{p}), \vec{p}, \varphi(\vec{p}), \mathrm{d}^3p

可以得到动量表象下的动量算符为

p^=p,p^x=px,p^y=py,p^z=pz \hat{\vec{p}} = \vec{p}, \kern 12pt \hat{p}_x = p_x, \kern 12pt \hat{p}_y = p_y, \kern 12pt \hat{p}_z = p_z

通过表象的转换,可以推得坐标表象下的动量算符为

p^=i,p^x=ix,p^y=iy,p^z=iz \hat{\vec{p}} = -\mathrm{i}\hbar\nabla, \kern 12pt \hat{p}_x = -\mathrm{i}\hbar \frac{\partial}{\partial x}, \kern 12pt \hat{p}_y = -\mathrm{i}\hbar \frac{\partial}{\partial y}, \kern 12pt \hat{p}_z = -\mathrm{i}\hbar \frac{\partial}{\partial z}

动量表象下的坐标算符为

r^=ip,x^=ipx,y^=ipy,z^=ipz \hat{\vec{r}} = \mathrm{i}\hbar \frac{\partial}{\partial \vec{p}}, \kern 12pt \hat{x} = \mathrm{i}\hbar \frac{\partial}{\partial p_x}, \kern 12pt \hat{y} = \mathrm{i}\hbar \frac{\partial}{\partial p_y}, \kern 12pt \hat{z} = \mathrm{i}\hbar \frac{\partial}{\partial p_z}

注:梯度、散度、旋度的介绍以及其在各种坐标系下的表示,可参考魏斌老师的课件。这里给出柱坐标与球坐标下的梯度算符,

柱坐标 (r,ϕ,z)(r,\phi,z)

f=frer+1rfθeθ+fzez \nabla f = \frac{\partial f}{\partial r}\vec{e}r + \frac{1}{r}\frac{\partial f}{\partial \theta}\vec{e}\theta + \frac{\partial f}{\partial z}\vec{e}_z

球坐标 (r,θ,φ)(r,\theta,\varphi)

f=frer+1rfθeθ+1rsinθfφeφ \nabla f = \frac{\partial f}{\partial r}\vec{e}r + \frac{1}{r}\frac{\partial f}{\partial \theta}\vec{e}\theta + \frac{1}{r\sin\theta}\frac{\partial f}{\partial \varphi}\vec{e}_\varphi

力学量算符

对于有经典对应的力学量,例如动能、势能和轨道角动量,由经典力学中的函数形式假定量子力学中的算符形式,可以由坐标算符与动量算符通过运算得到,即

A=A(r,p)A^=A(r^,p^) A = A(\vec{r},\vec{p}) \Longrightarrow \hat{A} = A(\hat{\vec{r}},\hat{\vec{p}})

如一维谐振子的能量算符

H=(px)22m+12kx2H^=(p^x)22m+12kx^2 H = \frac{(p_x)^2}{2m} + \frac12kx^2 \Longrightarrow \hat{H} = \frac{(\hat{p}_x)^2}{2m} + \frac12k\hat{x}^2

如粒子的轨道角动量算符

L=r×pL^=r^×p^=ijkx^y^z^p^xp^yp^z \vec{L} = \vec{r} \times \vec{p} \Longrightarrow \hat{\vec{L}} = \hat{\vec{r}} \times \hat{\vec{p}} =\begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \ \hat{x} & \hat{y} & \hat{z} \ \hat{p}_x & \hat{p}_y & \hat{p}_z \end{vmatrix}

L^x=y^p^zz^p^yL^y=z^p^xx^p^zL^z=x^p^yy^p^x \hat{L}_x = \hat{y}\hat{p}_z - \hat{z}\hat{p}_y \\kern 12pt\ \hat{L}_y = \hat{z}\hat{p}_x - \hat{x}\hat{p}_z \\kern 12pt\ \hat{L}_z = \hat{x}\hat{p}_y - \hat{y}\hat{p}_x

在坐标表象下,上述算符的表达式为

H^=22md2dx2+12kx2 \hat{H} = -\frac{\hbar^2}{2m} \frac{\mathrm{d}^2}{\mathrm{d}x^2} + \frac12kx^2

L^=r×(i)=ijkxyzixiyiz \hat{\vec{L}} = \vec{r} \times (-\mathrm{i}\hbar\nabla) =\begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \ x & y & z \ -\mathrm{i}\hbar \frac{\partial}{\partial x} & -\mathrm{i}\hbar \frac{\partial}{\partial y} & -\mathrm{i}\hbar \frac{\partial}{\partial z} \end{vmatrix}

L^x=i(yzzy)L^y=i(zxxz)L^z=i(xyyx) \hat{L}_x = -\mathrm{i}\hbar (y\frac{\partial}{\partial z}-z\frac{\partial}{\partial y}) \\kern 12pt\ \hat{L}_y = -\mathrm{i}\hbar (z\frac{\partial}{\partial x}-x\frac{\partial}{\partial z}) \\kern 12pt\ \hat{L}_z = -\mathrm{i}\hbar (x\frac{\partial}{\partial y}-y\frac{\partial}{\partial x})

对于已归一化的波函数,力学量 AA 在坐标表象与动量表象下的平均值表达式分别为

Aˉ=+ψ(r)A(r,i)ψ(r)d3r \bar{A} = \int_{-\infty}^{+\infty} \psi^*(\vec{r}), A(\vec{r},-\mathrm{i}\hbar\nabla), \psi(\vec{r}), \mathrm{d^3}r

Aˉ=+φ(p)A(ip,p)φ(p)d3p \bar{A} = \int_{-\infty}^{+\infty} \varphi^*(\vec{p}), A(\mathrm{i}\hbar \frac{\partial}{\partial \vec{p}},\vec{p}), \varphi(\vec{p}), \mathrm{d^3}p

1.2 Schrödinger方程

Schrödinger方程的建立与能量本征方程

单粒子运动的Schrödinger方程

在势场 U(r)U(\vec{r}) 中的例子的波函数 ψ(r,t)\psi(\vec{r},t) 满足以下Schrödinger波动方程

iψt=[22m2+U(r)]ψ \mathrm{i}\hbar\frac{\partial\psi}{\partial t} = \left[-\frac{\hbar^2}{2m}\nabla^2 + U(\vec{r})\right]\psi

假设势能 UU 不显含 tt ,上述方程可以使用分离变量法求解,即令

ψ(r,t)=ψ(r)f(t) \psi(\vec{r},t) = \psi(\vec{r})f(t)

代入原方程,分离变量,可得:

if(t)dfdt=1ψ(r)[22m2+U(r)]ψ(r)=E \frac{\mathrm{i}\hbar}{f(t)} \frac{\mathrm{d}f}{\mathrm{d}t} = \frac{1}{\psi(\vec{r})} \left[-\frac{\hbar^2}{2m}\nabla^2 + U(\vec{r})\right] \psi(\vec{r}) = E

其中 EE 是既不依赖于 tt ,也不依赖于 r\vec{r}, 的常数,首先考虑方程

if(t)dfdt=E \frac{\mathrm{i}\hbar}{f(t)} \frac{\mathrm{d}f}{\mathrm{d}t} = E

解得

f(t)eiEt f(t) \sim \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}Et}

则Schrödinger波动方程的特解为

ψ(r,t)=ψE(r) eiEt \psi(\vec{r},t) = \psi_E(\vec{r})\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}Et}

其中 ψE(r)\psi_E(\vec{r}) 满足以下方程

[22m2+U(r)]ψE(r)=EψE(r) \left[-\frac{\hbar^2}{2m}\nabla^2 + U(\vec{r})\right] \psi_E(\vec{r}) = E\psi_E(\vec{r})

对于此不含时Schrödinger方程(又称为定态Schrödinger方程),在某些条件下(特别是束缚态边条件),只有某些离散的 EE 值所对应的解才是物理上可以接受的,这些 EE 值称为体系的能量本征值(energy eigenvalue),而相应的解 ψE(r)\psi_E(\vec{r}) 称为能量本征函数(energy eigenfunction),该方程也称为势场 U(r)U(\vec{r}) 中粒子的能量本征方程。不同的能量本征值相应的本征函数是正交归一化的(设 EE 取离散值),即

(ψE,ψE)=δEE={1,E=E0,EE (\psi_E,\psi_{E'}) = \delta_{EE'} = \begin{cases} 1, & E = E' \ 0, & E \ne E' \end{cases}

引入Hamilton算符 H^\hat{H} (对于一个粒子在势场 U(r)U(\vec{r}) 中运动的情况, H^=22m2+U(r)\hat{H} = -\frac{\hbar^2}{2m} \nabla^2 + U(\vec{r}) ),则可以得到Schrödinger方程的普遍表达

iψt=H^ψ \mathrm{i}\hbar \frac{\partial\psi}{\partial t} = \hat{H}\psi

H^\hat{H} 不显含 tt 时,体系的能量是守恒量,此时的能量本征方程

H^ψ=Eψ \hat{H}\psi = E\psi

多粒子体系Schrödinger方程

设体系由 NN 个粒子组成,粒子质量分别为 mi (i=1,2,3,,N)m_i\ (i=1,2,3,\cdots,N) ,第 ii 个粒子收到的外势场为 Ui(ri)U_i(\vec{r}_i) ,粒子之间的相互作用为 V(r1,r2,,rN)V(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N) ,体系的波函数用 ψ(r1,r2,,rN,t)\psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N,t) 表示,则含时Schrödinger方程表示为

itψ(r1,r2,,rN,t)=[i=1N(22mii2+Ui(ri))+V(r1,r2,,rN)]ψ(r1,r2,,rN,t) \mathrm{i}\hbar\frac{\partial}{\partial t} \psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}N,t) = \left[\sum{i=1}^N \left(-\frac{\hbar^2}{2m_i}\nabla_i^2 + U_i(\vec{r}_i)\right) + V(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N) \right] \psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N,t)

其中

i2=2xi2+2yi2+2zi2 \nabla_i^2 = \frac{\partial^2}{\partial x_i^2} + \frac{\partial^2}{\partial y_i^2} + \frac{\partial^2}{\partial z_i^2}

不含时Schrödinger方程表示为

Eψ(r1,r2,,rN,t)=[i=1N(22mii2+Ui(ri))+V(r1,r2,,rN)]ψ(r1,r2,,rN,t) E\psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}N,t) = \left[\sum{i=1}^N \left(-\frac{\hbar^2}{2m_i}\nabla_i^2 + U_i(\vec{r}_i)\right) + V(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N) \right] \psi(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N,t)

在该体系中,Hamilton算符

H^=i=1N(22mii2+Ui(ri))+V(r1,r2,,rN) \hat{H} = \sum_{i=1}^N \left(-\frac{\hbar^2}{2m_i}\nabla_i^2 + U_i(\vec{r}_i)\right) + V(\vec{r}_1,\vec{r}_2,\cdots,\vec{r}_N)

Schrödinger方程的讨论

定域的概率守恒
粒子的空间概率密度

ρ(r,t)=ψ(r,t)2=ψ(r,t)ψ(r,t) \rho(\vec{r},t) = |\psi(\vec{r},t)|^2 = \psi^*(\vec{r},t)\psi(\vec{r},t)

概率流密度矢量

j(r,t)=i2m(ψψψψ)=12m(ψp^ψψp^ψ) \vec{j}(\vec{r},t) = -\frac{\mathrm{i}\hbar}{2m} (\psi^\nabla\psi - \psi\nabla\psi^) = \frac{1}{2m} (\psi^\hat{\vec{p}}\psi - \psi\hat{\vec{p}}\psi^)

连续性方程

ρt+j=0 \frac{\partial\rho}{\partial t} + \nabla\cdot\vec{j} = 0

该式对任意闭区域 τ\tau 的积分为

ddtτρdτ=SjdS \frac{\mathrm{d}}{\mathrm{d}t} \int_\tau \rho \mathrm{d}\tau = -\oint_S \vec{j} \cdot \mathrm{d}\vec{S}

该等式左边表示在闭区域 τ\tau 中找到例子的总概率(或粒子数)在单位时间内的增量,而右边则便是单位时间内通过 τ\tau 的封闭表面 SS 而流入 τ\tau 内的概率(粒子数),所以该式表达了概率(粒子数)守恒

在该积分表达式中,如果令 τ\tau\to\infty (即取全空间),由于任何真实的波函数应满足平方可积的条件,可以证明等式右侧的积分趋于零(也可以认为是在无穷远处不存在粒子的注入或流出,即不存在净粒子流),故

ddt()ρdτ=0()ρdτ=Const \frac{\mathrm{d}}{\mathrm{d}t} \int_{(全)} \rho \mathrm{d}\tau = 0 \Longrightarrow \int_{(全)} \rho \mathrm{d}\tau = \mathrm{Const}

这表明粒子在全空间的总概率守恒,即粒子既未产生,也未湮没。

初值问题

由于Schrödinger方程只含波函数 ψ(r,t)\psi(\vec{r},t) 对时间的一次微商,只要在初始时刻( t=0t=0 )体系的状态 ψ(r,0)\psi(\vec{r},0) 给定,则以后任何时态 tt 的状态 ψ(r,t)\psi(\vec{r},t) 原则上就完全确定了。

以下给出自由粒子的初值问题的解法:

对于自由粒子,其满足如下Schrödinger方程:

iψt=22m2ψ \mathrm{i}\hbar\frac{\partial\psi}{\partial t} = -\frac{\hbar^2}{2m}\nabla^2\psi

描述自由粒子的一般状态的波函数,具有波包的形式,可以视为许多平面单色波的叠加,即

ψ(r,t)=1(2π)32+φ(p) ei(prEt) d3p \psi(\vec{r},t) = \frac{1}{(2\pi\hbar)^{\frac{3}{2}}} \int_{-\infty}^{+\infty} \varphi(\vec{p})\ \mathrm{e}^{\frac{\mathrm{i}}{\hbar}(\vec{p}\cdot\vec{r}-Et)}\ \mathrm{d}^3p

式中 E=p22mE = \frac{p^2}{2m} ,其满足上述Schrödinger方程,其初态波函数为

ψ(r,0)=1(2π)32+φ(p) eipr d3p \psi(\vec{r},0) = \frac{1}{(2\pi\hbar)^{\frac{3}{2}}} \int_{-\infty}^{+\infty} \varphi(\vec{p})\ \mathrm{e}^{\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{r}}\ \mathrm{d}^3p

其中 φ(p)\varphi(\vec{p})ψ(r,t)\psi(\vec{r},t) 的Fourier展开的波幅,它并不依赖于 tt ,上式的逆变换为

φ(p)=1(2π)32+ψ(r,0) eipr d3r \varphi(\vec{p}) = \frac{1}{(2\pi\hbar)^{\frac{3}{2}}} \int_{-\infty}^{+\infty} \psi(\vec{r},0)\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}\vec{p}\cdot\vec{r}}\ \mathrm{d}^3r

φ(p)\varphi(\vec{p}) 由初态 ψ(r,0)\psi(\vec{r},0) 完全确定,可得

ψ(r,t)=1(2π)3+d3r+d3p eip(rr)iEt ψ(r,0) \psi(\vec{r},t) = \frac{1}{(2\pi\hbar)^3} \int_{-\infty}^{+\infty} \mathrm{d}^3r' \int_{-\infty}^{+\infty} \mathrm{d}^3p\ \mathrm{e}^{\frac{\mathrm{i}}{\hbar} \vec{p}\cdot(\vec{r}-\vec{r}) - \frac{\mathrm{i}}{\hbar} Et}\ \psi(\vec{r},0)

这样,体系的初始状态 ψ(r,0)\psi(\vec{r},0) 完全决定了以后任何时刻 tt 的状态 ψ(r,t)\psi(\vec{r},t)

定态与非定态

若在初始状态( t=0t=0 )体系处于某一个能量本征态 ψ(r,0)=ψE(r)\psi(\vec{r},0) = \psi_E(\vec{r}) ,则

ψ(r,t)=ψE(r) eiEt \psi(\vec{r},t) = \psi_E(\vec{r})\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}Et}

该波函数所描述的态,称为定态(stationary state)(体系的能量有确定值的状态,各种力学性质不随时间而改变);由若干个能量不同的本征态叠加所形成的态称为非定态(nonstationary state)。

ψ(r,t)=ECE ψE(r) eiEt \psi(\vec{r},t) = \sum_{E} C_E\ \psi_E(\vec{r})\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}Et}

处于定态下粒子具有以下特征:

  • 粒子在空间的概率密度 ρ(r)\rho(\vec{r}) 以及概率流密度 j(r)\vec{j}(\vec{r}) 不随时间改变。

  • 任何(不显含 tt 的)力学量的平均值不随时间改变。

  • 任何(不显含 tt 的)力学量的测量值概率分布不随时间改变。

关于Schrödinger方程的一些说明
  • Schrödinger方程是量子力学的一个基本假定,不能从其他更根本的假定来证明,其正确性由在各种具体情况下从方程得出的结论和实验结果比较来验证.。

  • Schrödinger方程是线性偏微分方程,满足“状态叠加原理”对波函数的要求,其解波函数是一个复函数。

  • Schrödinger方程是非相对论粒子的、且不发生实物粒子产生和湮灭的情况下,波函数满足的方程。

1.3 量子态叠加原理

量子态叠加原理

如果 ψ1,ψ2,,ψn\psi_1,\psi_2,\cdots,\psi_n 是体系的可能状态,那么 ψ=nCnψn\psi = \sum_n C_n\psi_nCnC_n 为复常数)也是体系的可能状态。

对于一个指定的量子体系,如果我们找到了它的“完备的基本状态”,那么任何状态都可以由这些基本状态叠加而得到。

测量与波函数坍缩

假设粒子处于非定态

ψ(r,t)=nCn ψn(r) eiEnt \psi(\vec{r},t) = \sum_{n} C_n\ \psi_n(\vec{r})\ \mathrm{e}^{-\frac{\mathrm{i}}{\hbar}E_nt}

即很多能量本征值 En (n=1,2,3,)E_n\ (n=1,2,3,\cdots) 的本征态 ψn\psi_n 的叠加,则在测量例子能量时,求和中包含的所有能量本征值 EnE_n 都有可能出现,出现的概率分别为 Cn2|C_n|^2 (应满足归一化条件 nCn2=1\sum_n |C_n|^2 = 1 ) 。当测量结果为某个能量本征值 EnE_n 时,粒子的状态就变为相应的能量本征态 ψn\psi_n ,按照von Neumann的看法,量子力学中把此称为量子态坍缩,即在测量的过程中,粒子的状态由叠加态坍缩为某一能量本征态。

在任意状态 ψ\psi 上,对任意力学量 AA 进行足够多次的测量,所得结果的平均值(期望值)为

Aˉ=(ψ,A^ψ)(ψ,ψ) \bar{A} = \frac{(\psi,\hat{A}\psi)}{(\psi,\psi)}

若满足归一化条件 (ψ,ψ)=1(\psi,\psi)=1 ,则

Aˉ=(ψ,A^ψ) \bar{A} = (\psi,\hat{A}\psi)

若归一化的 ψ(r)\psi(\vec{r}) 不是算符 A^\hat{A} 的本征函数,只要 AA 是可观察的力学量,对于

ψ(r)=nCn ψn(r)(nCn2=1) \psi(\vec{r}) = \sum_{n} C_n\ \psi_n(\vec{r}) \kern 24pt \left(\sum_n |C_n|^2 = 1\right)

若在每个本征态有 A^ψn=Anψn\hat{A}\psi_n = A_n\psi_n ,则

Aˉ=(ψ,A^ψ)(ψ,ψ)=([nCn ψn(r)],A^[nCn ψn(r)])([nCn ψn(r)],[nCn ψn(r)])=nCn2AnnCn2=nCn2An \bar{A} = \frac{(\psi,\hat{A}\psi)}{(\psi,\psi)} = \frac{([\sum_{n} C_n\ \psi_n(\vec{r})],\hat{A}[\sum_{n} C_n\ \psi_n(\vec{r})])}{([\sum_{n} C_n\ \psi_n(\vec{r})],[\sum_{n} C_n\ \psi_n(\vec{r})])} = \frac{\sum_{n} |C_n|^2 A_n}{\sum_{n} |C_n|^2} = \sum_{n} |C_n|^2 A_n

第2章 一维势场中的粒子

2.1 一维运动问题的一般分析

问题分类与基本概念

一维定态的分类:束缚态与非束缚态

本章主要是用Schrödinger方程来处理一维粒子的能量本征态问题,一般分为两类问题:

  • 束缚态问题束缚态(bound state)是指粒子局限在有限的空间中,即粒子在无穷远处出现的概率等于零的状态,即当 x±x\to\pm\infty 时,有 ψ(x)0\psi(x)\to0 ;而非束缚态(或称散射态)指粒子可以出现在无穷远处的状态,即当 x+x\to+\inftyxx\to-\infty 时, ψ(x)0\psi(x)\ne0 。束缚态问题会给出势场函数 V(x)V(x) ,需要求出波函数 ψ(x)\psi(x) 以及能量本征值 EE (通常是离散的 EnE_n )。

  • 散射问题:会给出势场函数 V(x)V(x) 与能量 EE ,需要求出波函数 ψ(x)\psi(x)

在求解上述两个问题的能量本征方程时,要根据具体物理问题的边界条件来定解。(束缚态问题还有着 limx ψ(x)=0\lim_{x\to\infty}\ \psi(x)=0 的无穷远处条件)

简并度

如果对一个给定的能量 EE ,只有一个线性独立的波函数存在,则称该能级是非简并的;否则称它是简并的,其线性独立的波函数的个数称为它的简并度

宇称

定义一维粒子的空间反射算符 PP

Pψ(x)=ψ(x) P \psi(x) = \psi(-x)

其对应的本征方程为

Pψ(x)=πψ(x) P \psi(x) = \pi \psi(x)

定义宇称(parity)为空间反射算符的本征值 π\pi ,可以证明,空间反演算符只有 ±1\pm1 两个本征值,

Pψ(x)=ψ(x)={ψ(x),π=+1,偶(正)宇称ψ(x),π=1,奇(负)宇称 P\psi(x) = \psi(-x) = \begin{cases} \psi(x), & \pi=+1, & 偶(正)宇称 \ -\psi(x), & \pi=-1, & 奇(负)宇称 \end{cases}

即空间反射不变的波函数具有偶(正)宇称(even parity);变号的波函数具有奇(负)宇称(odd parity);还有一些波函数没有确定的宇称,它们不是空间反射算符的本征态。

一维定态Schrödinger方程解的一般性质

质量为 mm 的粒子,沿 xx 方向运动,势能为 V(x)V(x) ,则定态Schrödinger方程表示为

[22md2dx2+V(x)]ψ(x)=Eψ(x) \left[-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}+V(x)\right] \psi(x) = E \psi(x)

在量子力学中,如果不作特别的声明,都认为 V(x)V(x) 取实值,即 V(x)=V(x)V(x) = V^*(x)

定理1 共轭定理

ψ(x)\psi(x) 是定态Schrödinger方程的一个解,对应的能量本征值为 EE ,则 ψ(x)\psi^*(x) 也是该方程的一个解,对应的能量也是 EE

证明

对定态Schrödinger方程取复共轭,可得

[22md2dx2+V(x)]ψ(x)=Eψ(x) \left[-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}+V(x)\right] \psi^(x) = E \psi^(x)

显然 ψ(x)\psi^*(x) 也是定态Schrödinger方程的解,且对应的能量本征值为 EE

推论

假设对应于能量的某个本征值 EE ,定态Schrödinger方程的解无简并,则可取为实解(除了一个无关紧要的常数因子外)。

定理2

对应于能量的某个本征值 EE ,总可以找到定态Schrödinger方程的一组实解,凡属于 EE 的任何解,均可表示为这一组实解的线性叠加

对于能级有简并的情况,要用到此定理;通过定理1和定理2,可以说明定态Schrödinger方程的基本解组可全取为实解。

证明

假设 ψ(x)\psi(x) 是定态Schrödinger方程的一个解:

如果 ψ(x)\psi(x) 是实解,则可把它归入实解的集合中去;

如果 ψ(x)\psi(x) 是复解,则由定理1可知, ψ(x)\psi^*(x) 也是方程的解,且同属于能量本征值 EE 。根据线性微分方程解的叠加性定理,

φ(x)=ψ(x)+ψ(x),χ(x)=i[ψ(x)ψ(x)] \varphi(x) = \psi(x) + \psi^(x), \kern 1em \chi(x) = -\mathrm{i} [\psi(x) - \psi^(x)]

也是方程同属于能量 EE 的解,且彼此独立。 φ(x)\varphi(x)χ(x)\chi(x) 均为实解, 而 ψ(x)\psi(x)ψ(x)\psi^*(x) 均可表示为其线性叠加,即

ψ=12(φ+iχ),ψ=12(φiχ) \psi = \frac12 (\varphi + \mathrm{i} \chi), \kern 1em \psi^* = \frac12 (\varphi - \mathrm{i} \chi)

定理3 反射定理

设势能函数 V(x)V(x) 具有空间反射不变性,即 V(x)=V(x)V(x)=V(-x),那么若 ψ(x)\psi(x) 是定态Schrödinger方程对应于能量本征值 EE 的解,则 ψ(x)\psi(-x) 也是该方程对应于能量 EE 的解。

证明

xxx\longrightarrow -x 时,有

d2dx2d2d(x)2=d2dx2,V(x)V(x)=V(x) \frac{\mathrm{d}^2}{\mathrm{d}x^2} \longrightarrow \frac{\mathrm{d}^2}{\mathrm{d}(-x)^2} = \frac{\mathrm{d}^2}{\mathrm{d}x^2}, \kern 1em V(x) \longrightarrow V(-x) = V(x)

则定态Schrödinger方程转化为

[22md2dx2+V(x)]ψ(x)=Eψ(x) \left[-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}+V(x)\right] \psi(-x) = E \psi(-x)

显然 ψ(x)\psi(-x) 也是定态Schrödinger方程的解,且对应的能量本征值为 EE

推论

如果对应于某能量 EE ,定态Schrödinger方程的解无简并,则解必有确定的宇称

因为此时 ψ(x)\psi(x)ψ(x)\psi(-x) 代表同一个解,它们最多可以差一个常数因子 π\pi ,即 Pψ(x)=ψ(x)=πψ(x)P\psi(x) = \psi(-x) = \pi\psi(x)

定理4

设势能函数 V(x)V(x) 具有空间反射不变性,即 V(x)=V(x)V(x)=V(-x),则对应于任何一个能量本征值 EE ,总可以找到定态Schrödinger方程的一组解 (每个解都有确定的宇称),而属于能量本征值 EE 的任何解,都可用它们来展开.

对于能级有简并的情况,能量本征态并不一定就具有确定宇称,此时,可以用该定理来处理;通过定理3和定理4,可以说明定态Schrödinger方程的基本解组可全取为具有确定宇称的解。

证明

假设 ψ(x)\psi(x) 是定态Schrödinger方程的一个解:

如果 ψ(x)\psi(x) 有确定的宇称,则可把它归入有确定的宇称的解集中去;

如果 ψ(x)\psi(x) 无确定的宇称,则由定理3可知, ψ(x)\psi(-x) 也是方程的解,且同属于能量本征值 EE ,但不同于 ψ(x)\psi(x) 。根据线性微分方程解的叠加性定理,

f(x)=ψ(x)+ψ(x),g(x)=ψ(x)ψ(x) f(x) = \psi(x) + \psi(-x),\kern 1em g(x) = \psi(x) - \psi(-x)

也是方程同属于能量 EE 的解,且彼此独立。 f(x)f(x)g(x)g(x) 均具有确定宇称: f(x)=f(x)f(-x)=f(x) , g(x)=g(x)g(-x)=-g(-x) ;而 ψ(x)\psi(x)ψ(x)\psi(-x) 均可表示为其线性叠加,即

ψ(x)=12[f(x)+g(x)],ψ(x)=12[f(x)g(x)] \psi(x) = \frac12[f(x)+g(x)], \kern 1em \psi(-x) = \frac12[f(x)-g(x)]

定理5

对于阶梯方位势

V(x)={V1,x<aV2,x>a V(x) = \begin{cases} V_1, & x<a \ V_2, & x>a \end{cases}

(V2V1)(V_2-V_1) 有限,则能量本征函数 ψ\psi 及其导数 ψ(x)\psi'(x) 必定是连续的;但若 V2V1|V_2-V_1|\to\infty ,则该定理不成立。

证明

根据定态Schrödinger方程

ψ(x)=2m2[EV(x)]ψ(x) \psi''(x) = -\frac{2m}{\hbar^2} \left[E-V(x)\right] \psi(x)

V(x)V(x) 连续的区域,由 ψ(x)\psi''(x) 存在可以推出 ψ(x)\psi(x)ψ(x)\psi'(x) 是连续的。

V(x)V(x) 发生阶梯形跳跃处, V(x)ψ(x)V(x)\psi(x) 发生跃变,但变化是有限的,在 xax\sim a 邻域对上述方程积分,得

limε0+aεa+εψ(x)dx=ψ(a+0+)ψ(a0+)=2m2limε0+aεa+εdx[EV(x)]ψ(x) \lim_{\varepsilon\to0^+} \int_{a-\varepsilon}^{a+\varepsilon} \psi''(x) \mathrm{d}x = \psi'(a+0^+) - \psi'(a-0^+) = -\frac{2m}{\hbar^2} \lim_{\varepsilon\to0^+} \int_{a-\varepsilon}^{a+\varepsilon} \mathrm{d}x [E-V(x)] \psi(x)

由于 [EV(x)]ψ(x)[E-V(x)]\psi(x) 是有限的,当 ε0+\varepsilon\to0^+ 时,上式右边积分趋于零,因此

ψ(a+0+)=ψ(a0+) \psi'(a+0^+) = \psi'(a-0^+)

ψ(x)\psi'(x)V(x)V(x) 的跳跃点 x=ax=a 处是连续的,因而 ψ(x)\psi(x) 也是连续的。

定理6 Wronskian定理

ψ1(x)\psi_1(x)ψ2(x)\psi_2(x) 均为定态Schrödinger方程属于同一能量 EE 的解,则

ψ1ψ2ψ2ψ1=Const(x无关) \psi_1\psi'_2 - \psi_2\psi'_1 = \mathrm{Const}(与x无关)

其中 ψ1ψ2ψ2ψ1\psi_1\psi'_2 - \psi_2\psi'_1 称为 ψ1(x)\psi_1(x)ψ2(x)\psi_2(x) 的Wronskian行列式,即

W[ψ1,ψ2](x)=ψ1(x)ψ2(x)ψ1(x)ψ2(x) W\psi_1,\psi_2 = \begin{vmatrix} \psi_1(x) & \psi_2(x) \ \psi'_1(x) & \psi'_2(x) \end{vmatrix}

证明

由定态Schrödinger方程可得

ψ1=2m2[EV(x)]ψ1ψ2=2m2[EV(x)]ψ2 \psi''_1 = -\frac{2m}{\hbar^2} \left[E-V(x)\right] \psi_1 \ \kern1em \ \psi''_2 = -\frac{2m}{\hbar^2} \left[E-V(x)\right] \psi_2

ψ1\psi_1 ×\times 下式 - ψ2\psi_2 ×\times 上式,可得

ψ1ψ2ψ2ψ1=0 \psi_1\psi''_2 - \psi_2\psi''_1 = 0

ψ1ψ2ψ2ψ1=ψ1ψ2+ψ1ψ2ψ1ψ2ψ2ψ1=(ψ1ψ2ψ2ψ1)=0 \psi_1\psi''_2 - \psi_2\psi''_1 = \psi_1\psi''_2 + \psi'_1\psi'_2 - \psi'_1\psi'_2 - \psi_2\psi''_1 = (\psi_1\psi'_2 - \psi_2\psi_1)' = 0

积分,得

ψ1ψ2ψ2ψ1=Const(x无关) \psi_1\psi'_2 - \psi_2\psi'_1 = \mathrm{Const}(与x无关)

推论

对于束缚态,当 xx\to\infty 时, ψ0\psi\to0 ,所以该定理中的常数必为 00 ,因此对于同属于能量 EE 的任何两个束缚态波函数 ψ1\psi_1ψ2\psi_2

ψ1ψ2=ψ2ψ1 \psi_1\psi'_2 = \psi_2\psi'_1

定理7 不简并定理

设粒子在规则势场 V(x)V(x) (无奇点)中运动,如存在束缚态,则必定是非简并的。

注:对于常见的不规则势阱,在绝大多数情况下(如无限深方势阱、 δ\delta 势阱等),该定理也成立;但对于某些不规则势阱,如一维氢原子( V(x)1xV(x) \propto -\frac{1}{|x|}),除基态外,其他束缚态简并度均为 22 ,其特征是波函数的节点(指 ψ(x)=0\psi(x)=0 的点)出现在 V(x)V(x) 的奇异点处,两个简并态具有不同宇称。

证明

ψ1\psi_1ψ2\psi_2 是定态Schrödinger方程属于同一能量 EE 的两个束缚态解,则

ψ1ψ2=ψ2ψ1 \psi_1\psi'_2 = \psi_2\psi'_1

在不包含 ψ1(x)\psi_1(x)ψ2(x)\psi_2(x) 节点的区域中,等式左右两边同除以 ψ1ψ2\psi_1\psi_2 ,得

ψ2ψ2=ψ1ψ1 \frac{\psi'_2}{\psi_2} = \frac{\psi'_1}{\psi_1}

(lnψ1ψ2)=0 \left(\ln{\frac{\psi_1}{\psi_2}}\right)' = 0

积分得

lnψ1ψ2=lnC(C是与x无关的常数) \ln{\frac{\psi_1}{\psi_2}} = \ln C \kern 1em (C是与x无关的常数)

ψ1=Cψ2 \psi_1 = C \psi_2

这表明 ψ1\psi_1ψ2\psi_2 代表同一个量子态,即能级不简并。

2.2 束缚态问题:一维无限深势阱和有限深势阱、一维谐振子

一维无限深方势阱

模型描述与结论

一维无限深方势阱表示为

V(x)={0,0<x<a,0<x,x>a V(x) = \begin{cases} 0, & 0<x<a \ \infty, & 0<x,x>a \end{cases}

在该势阱中的质量为 mm 的粒子,能量是量子化的,即构成的能谱是离散的,体系的能量本征值

En=2π2n22ma2(n=1,2,3,) E_n = \frac{\hbar^2\pi^2n^2}{2ma^2} \kern 1em (n=1,2,3,\cdots)

对应的能量本征函数为

ψn(x)={2asin(nπxa),0<x<a0,x<0,x>a \psi_n(x) = \begin{cases} \sqrt{\frac{2}{a}} \sin\left(\frac{n\pi x}{a}\right), & 0<x<a \ 0, & x<0,x>a \end{cases}

注:若将一维无限深方势阱表示为

V(x)={0,x<a2,x>a2 V(x) = \begin{cases} 0, & |x|<\frac{a}{2} \ \infty, & |x|>\frac{a}{2} \end{cases}

则能量本征值不变,能量本征函数变为

ψn(x)={{a2cos(nπxa),n=1,3,5,,(偶宇称)a2sin(nπxa),n=2,4,6,,(奇宇称)x<a20,x<a2 \psi_n(x) = \begin{cases} \begin{cases} \sqrt{\frac{a}{2}} \cos\left(\frac{n\pi x}{a}\right) , & n=1,3,5,\cdots,(偶宇称) \ \sqrt{\frac{a}{2}} \sin\left(\frac{n\pi x}{a}\right) , & n=2,4,6,\cdots,(奇宇称) \end{cases} & |x|<\frac{a}{2} \ 0, & |x|<\frac{a}{2} \end{cases}

模型求解

在势阱内 (0<x<a)(0<x<a) ,能量本征方程为

d2dx2ψ(x)+2mE2ψ(x)=0 \frac{\mathrm{d}^2}{\mathrm{d}x^2} \psi(x) + \frac{2mE}{\hbar^2} \psi(x) = 0

其中粒子的能量 E>0E>0 ,令

k=2mE>0 k = \frac{\sqrt{2mE}}{\hbar} > 0

则能量本证方程可表示为

d2dx2ψ(x)+k2ψ(x)=0 \frac{\mathrm{d}^2}{\mathrm{d}x^2} \psi(x) + k^2 \psi(x) = 0

解得

ψ(x)=Asin(kx+δ) \psi(x) = A\sin(kx+\delta)

其中 AAδ\delta 为待定常数。因为势壁无限高,从物理上考虑,粒子不能透过势壁;按波函数的统计诠释,要求在阱壁上及阱外波函数为 00 。这样就得到了边界条件

ψ(0)=0,ψ(a)=0 \psi(0) = 0, \kern 1em \psi(a) = 0

ψ(0)=Asin(δ)=0\psi(0) = A\sin(\delta) = 0 ,可取 δ=0\delta = 0 ,则 ψ(x)=Asin(kx)\psi(x) = A\sin(kx) ,由 ψ(a)=Asin(ka)=0\psi(a) = A\sin(ka) = 0 ,可知

ka=nπk=nπa(n=1,2,3,) ka = n\pi \Longrightarrow k = \frac{n\pi}{a} \kern 1em (n=1,2,3,\cdots)

注: n=0n=0 给出的波函数 ψ(x)=0\psi(x)=0 ,无物理意义;而 nn 取负值与 nn 取对应的正值得到的波函数只相差一个常数 1-1 ,描述的是同一个量子态。

k=2mE=nπaEn=2π2n22ma2(n=1,2,3,) k = \frac{\sqrt{2mE}}{\hbar} = \frac{n\pi}{a} \Longrightarrow E_n = \frac{\hbar^2\pi^2n^2}{2ma^2} \kern 1em (n=1,2,3,\cdots)

则能量本征函数

ψn(x)=Asin(nπax) \psi_n(x) = A\sin\left(\frac{n\pi}{a}x\right)

归一化后可得

ψn(x)=2asin(nπax) \psi_n(x) = \sqrt{\frac{2}{a}}\sin(\frac{n\pi}{a}x)

讨论
能量本征值

En=2π2n22ma2n2(n=1,2,3,) E_n = \frac{\hbar^2\pi^2n^2}{2ma^2} \propto n^2 \kern 1em (n=1,2,3,\cdots)

能级的分布是不均匀的,能级越高,密度越小

ΔEn2π2ma2n,ΔEnEn=2nn0 \Delta E_n \approx \frac{\hbar^2\pi^2}{ma^2} n, \kern 1em \frac{\Delta E_n}{E_n} = \frac{2}{n} \overset{n\to\infty}{\longrightarrow} 0

nn 充分大时,可以认为能量连续。

最低能级不为零:

E1=2π22ma2>0 E_1 = \frac{\hbar^2\pi^2}{2ma^2} > 0

这可以用不确定性关系来解释:

粒子局限在无限深方势阱中,位置不确定度 Δxa\Delta x \sim a ,则动量不确定度 ΔpΔxa\Delta p \sim \frac{\hbar}{\Delta x} \sim \frac{\hbar}{a} ,故能量不能为零,

Ep22m(Δp)22m22ma20 E \sim \frac{p^2}{2m} \sim \frac{(\Delta p)^2}{2m} \sim \frac{\hbar^2}{2ma^2} \ne 0

概率密度

ρn(x)=ψn(x)2=2asin2(nπax) \rho_n(x) = |\psi_n(x)|^2 = \frac{2}{a}\sin^2\left(\frac{n\pi}{a}x\right)

无限深方势阱中较低几条能级的波函数

由该图也可以看出,除端点 (x=0,a)(x=0,a) 外,基态(能量最低态, n=1n=1 )波函数无节点,第 kk 激发态( k=n1k=n-1 )有 kk 个节点。

一维有限深对称方势阱

模型描述与结论

一维有限深对称方势阱表示为

V(x)={0,x<a2V0,x>a2 V(x) = \begin{cases} 0, & |x|<\frac{a}{2} \ V_0, & |x|>\frac{a}{2} \end{cases}

一维有限深对称方势阱

在该势阱中的质量为 mm 的粒子,讨论其处于束缚态 (0<E<V0)(0<E<V_0) 的情况。

k=2mE ,β=2m(V0E) k = \frac{\sqrt{2mE}}{\hbar}\ , \kern 1em \beta = \frac{\sqrt{2m(V_0-E)}}{\hbar}

引入无量纲参数

ξ=ka2 ,η=βa2 \xi = \frac{ka}{2}\ , \kern 1em \eta = \frac{\beta a}{2}

这两个无量纲参数满足一定的方程组,使得其取值是离散的,对应的能量本征态为

En=22ma2ξn2 E_n = \frac{2\hbar^2}{ma^2} \xi^2_n

对于偶宇称态

波函数形式为(可利用波函数连续性与归一化进一步求出 AACC

ψ(x)={Ceβx,x<a2Acos(kx),a2<x<a2Ceβx,x>a2 \psi(x) = \begin{cases} C\mathrm{e}^{\beta x}, & x<-\frac{a}{2} \ A\cos(kx), & -\frac{a}{2}<x<\frac{a}{2} \ C\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

无量纲参数满足方程组

{ξtanξ=ηξ2+η2=mV0a222 \left{\begin{matrix} \xi\tan\xi = \eta \ \xi^2+\eta^2 = \frac{mV_0a^2}{2\hbar^2} \end{matrix}\right.

对于奇宇称态

波函数形式为(可利用波函数连续性与归一化进一步求出 BBCC

ψ(x)={Ceβx,x<a2Bsin(kx),a2<x<a2Ceβx,x>a2 \psi(x) = \begin{cases} C\mathrm{e}^{\beta x}, & x<-\frac{a}{2} \ B\sin(kx), & -\frac{a}{2}<x<\frac{a}{2} \ -C\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

无量纲参数满足方程组

{ξcotξ=ηξ2+η2=mV0a222 \left{\begin{matrix} -\xi\cot\xi = \eta \ \xi^2+\eta^2 = \frac{mV_0a^2}{2\hbar^2} \end{matrix}\right.

模型求解

先考虑势阱外的情况,能量本征方程为

d2dx2ψ(x)β2ψ(x)=0 \frac{\mathrm{d}^2}{\mathrm{d}x^2} \psi(x) - \beta^2 \psi(x) = 0

其中 β=2m(V0E)\beta = \frac{\sqrt{2m(V_0-E)}}{\hbar} ,解得

ψ(x)=Ceβx+Deβx \psi(x) = C\mathrm{e}^{\beta x} + D\mathrm{e}^{-\beta x}

考虑束缚态边界条件,即在 xx\to\infty 处,要求 ψ(x)0\psi(x)\to0 ,则

ψ(x)={Ceβx,x<a2Deβx,x>a2 \psi(x) = \begin{cases} C\mathrm{e}^{\beta x}, & x<-\frac{a}{2} \ D\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

再考虑势阱内的情况,能量本征方程为

d2dx2ψ(x)+k2ψ(x)=0 \frac{\mathrm{d}^2}{\mathrm{d}x^2} \psi(x) + k^2 \psi(x) = 0

其中 k=2mEk = \frac{\sqrt{2mE}}{\hbar} ,解得

ψ(x)=Acos(kx)+Bsin(kx) \psi(x) = A\cos(kx) + B\sin(kx)

考虑到势阱具有空间反射不变性 V(x)=V(x)V(-x)=V(x) ,由定理3推论可知,束缚态能量本征函数(由定理7知其不简并)必具有确定的宇称,因此只能单独取 cos(kx)\cos(kx)sin(kx)\sin(kx) 形式,以下分别讨论。

对于偶宇称态

在势阱内 (x<a2)(|x|<\frac{a}{2})B=0B=0

ψ(x)=Acos(kx) \psi(x) = A\cos(kx) \kern 2em

在势阱外 (x>a2)(|x|>\frac{a}{2})C=DC=D ,

ψ(x)={Ceβx,x<a2Ceβx,x>a2 \psi(x) = \begin{cases} C\mathrm{e}^{\beta x}, & x<-\frac{a}{2} \ C\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

由定理5可知,波函数 ψ(x)\psi(x) 及导数 ψ(x)\psi'(x)x=a2|x|=\frac{a}{2} 处是连续的,由于波函数具有偶宇称,在 a2-\frac{a}{2}a2\frac{a}{2} 处的情况实际上是等效的,这里只用分析 x=a2x=\frac{a}{2} 的情况。

ψ(x)={Acos(kx),a2<x<a2Ceβx,x>a2ψ(x)={Aksin(kx),a2<x<a2Cβeβx,x>a2 \psi(x) = \begin{cases} A\cos(kx), & -\frac{a}{2}<x<\frac{a}{2} \ C\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases} \ \kern 1em \ \psi'(x) = \begin{cases} -Ak\sin(kx), & -\frac{a}{2}<x<\frac{a}{2} \ -C\beta\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

Acos(ka2)=Ceβa2Aksin(ka2)=Cβeβa2 A\cos(k\frac{a}{2}) = C\mathrm{e}^{-\beta\frac{a}{2}} \ \kern 1em \ -Ak\sin(k\frac{a}{2}) = -C\beta\mathrm{e}^{-\beta\frac{a}{2}}

两式相除,可消去 A,CA,C ,得

ktan(ka2)=β k\tan(k\frac{a}{2}) = \beta

注:也可以直接考虑 (lnψ)(\ln\psi)' 的连续性,从而直接消去 A,CA,C 这两个常数。

引入无量纲参数

ξ=ka2 ,η=βa2 \xi = \frac{ka}{2}\ , \kern 1em \eta = \frac{\beta a}{2}

可得

ξtanξ=η \xi\tan\xi = \eta

同时, ξ\xiη\eta 还满足

ξ2+η2=a24(k2+β2)=a24[2mE2+2m(V0E)2]=mV0a222 \xi^2+\eta^2 = \frac{a^2}{4}(k^2+\beta^2) = \frac{a^2}{4}\left[\frac{2mE}{\hbar^2} + \frac{2m(V_0-E)}{\hbar^2}\right] = \frac{mV_0a^2}{2\hbar^2}

整理即得 ξ\xiη\eta 满足方程组

{ξtanξ=ηξ2+η2=mV0a222 \left{\begin{matrix} \xi\tan\xi = \eta \ \xi^2+\eta^2 = \frac{mV_0a^2}{2\hbar^2} \end{matrix}\right.

对于奇宇称态

在势阱内 (x<a2)(|x|<\frac{a}{2})A=0A=0

ψ(x)=Bsin(kx) \psi(x) = B\sin(kx) \kern 2em

在势阱外 (x>a2)(|x|>\frac{a}{2})C=DC=-D ,

ψ(x)={Ceβx,x<a2Ceβx,x>a2 \psi(x) = \begin{cases} C\mathrm{e}^{\beta x}, & x<-\frac{a}{2} \ -C\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

波函数 ψ(x)\psi(x) 及导数 ψ(x)\psi'(x)x=a2x=\frac{a}{2} 处是连续的,

ψ(x)={Bsin(kx),a2<x<a2Ceβx,x>a2ψ(x)={Bkcos(kx),a2<x<a2Cβeβx,x>a2 \psi(x) = \begin{cases} B\sin(kx), & -\frac{a}{2}<x<\frac{a}{2} \ -C\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases} \ \kern 1em \ \psi'(x) = \begin{cases} Bk\cos(kx), & -\frac{a}{2}<x<\frac{a}{2} \ C\beta\mathrm{e}^{-\beta x}, & x>\frac{a}{2} \end{cases}

Bsin(ka2)=Ceβa2Bkcos(ka2)=Cβeβa2 B\sin(k\frac{a}{2}) = -C\mathrm{e}^{-\beta\frac{a}{2}} \ \kern 1em \ Bk\cos(k\frac{a}{2}) = C\beta\mathrm{e}^{-\beta\frac{a}{2}}

两式相除,可消去 B,CB,C ,得

kcot(ka2)=β -k\cot(k\frac{a}{2}) = \beta

引入无量纲参数

ξ=ka2 ,η=βa2 \xi = \frac{ka}{2}\ , \kern 1em \eta = \frac{\beta a}{2}

可得

ξtanξ=η \xi\tan\xi = \eta

同时, ξ\xiη\eta 还满足

ξ2+η2=mV0a222 \xi^2+\eta^2 =\frac{mV_0a^2}{2\hbar^2}

整理即得 ξ\xiη\eta 满足方程组

{ξcotξ=ηξ2+η2=mV0a222 \left{\begin{matrix} -\xi\cot\xi = \eta \ \xi^2+\eta^2 = \frac{mV_0a^2}{2\hbar^2} \end{matrix}\right.

讨论
无量纲参数满足的方程组的解

对于 ξ\xiη\eta 满足的方程组,可以采用图解法近似求解,无论是奇宇称态还是偶宇称态,方程组中的第二个方程的图象都是圆弧,半径为 mV0a222\sqrt{\frac{mV_0a^2}{2\hbar^2}} 。对于偶宇称态,图为:

有限深对称方势阱偶宇称态

对于奇宇称态,图为:

有限深对称方势阱奇宇称态

注:实际上对于确定的 m,a,V0m,a,V_0 ,图中的圆弧应该只有一条。

在一维有限深对称方势阱问题中,无论势阱多浅或多窄(即无论 V0a2V_0a^2 的值多小),偶宇称态的方程组都至少有一个根,这表明至少存在一个束缚态(即基态),其宇称为偶。

而对于奇宇称态的方程组,只有当

ξ2+η2=mV0a222(π2)2 \xi^2+\eta^2 = \frac{mV_0a^2}{2\hbar^2} \ge \left(\frac{\pi}{2}\right)^2

V0a2π222m V_0a^2 \ge \frac{\pi^2\hbar^2}{2m}

方程组才会有解,即才可能出现最低的奇宇称能级。

随着 V0a2V_0a^2 的增大,方程组的解的个数会逐渐增多,出现更高的激发态能级,宇称奇偶相间。由图可得,圆弧的半径 mV0a222\sqrt{\frac{mV_0a^2}{2\hbar^2}} 每增大 π2\frac{\pi}{2} ,两图中交点的总个数会增加一个,由此可以推得束缚态能级总数为

N=1+aπ2mV0 N = 1+\left\lfloor\frac{a}{\hbar\pi}\sqrt{2mV_0}\right\rfloor

能量本征值

ξ=ka2=a22mEEn=22ma2ξn2 \xi = \frac{ka}{2} = \frac{a}{2}\frac{\sqrt{2mE}}{\hbar} \Longrightarrow E_n = \frac{2\hbar^2}{ma^2} \xi_n^2

由上图可得

0<ξ1<π2<ξ2<π<<π2(n1)<ξn<π2n< 0<\xi_1<\frac{\pi}{2}<\xi_2<\pi<\cdots<\frac{\pi}{2}(n-1)<\xi_n<\frac{\pi}{2}n<\cdots

故有限深方势阱每个能级都比无限深方势阱的相应能级低一些:

En=22ma2ξn2<π222ma2n2 E_n = \frac{2\hbar^2}{ma^2} \xi_n^2 < \frac{\pi^2\hbar^2}{2ma^2} n^2

V0+V_0\to+\infty 时,有 ξnπ2n\xi_n\to\frac{\pi}{2}n ,则 Enπ222ma2n2E_n\to\frac{\pi^2\hbar^2}{2ma^2}n^2 ,即趋向于无限深方势阱的能级。

一维谐振子

模型描述与结论

取谐振子的平衡位置为坐标原点,并选原点为势能的零点,则以为谐振子的势能可以表示为

V(x)=12mω2x2 V(x) = \frac12m\omega^2x^2

其中 mm 为谐振子的质量, ω\omega 为经典谐振子的自然频率。理想的谐振子势是一个无限深势阱,只存在束缚态,谐振子的能量本征值

En=(n+12)ω(n=0,1,2,) E_n = (n+\frac12)\hbar\omega \kern 2em (n=0,1,2,\cdots)

正交归一化的能量本征函数

ψn(x)=Aneα2x22Hn(αx) \psi_n(x) = A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x)

其中 α=mω\alpha=\sqrt{\frac{m\omega}{\hbar}}Hn(x)\mathrm{H}_n(x) 为Hermite多项式,归一化系数为

An=απ2nn! A_n = \sqrt{\frac{\alpha}{\sqrt{\pi}2^n\cdot n!}}

模型求解

一维谐振子的能量本征方程为

[22md2dx2+12mω2x2]ψ(x)=Eψ(x) \left[-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}+\frac12m\omega^2x^2\right] \psi(x) = E \psi(x)

α=mω\alpha=\sqrt{\frac{m\omega}{\hbar}} ,并引进无量纲参量

ξ=αx,λ=E12ω \xi = \alpha x , \kern 1em \lambda = \frac{E}{\frac12\hbar\omega}

则方程可整理为

d2dξ2ψ+(λξ2)ψ=0 \frac{\mathrm{d}^2}{\mathrm{d}\xi^2} \psi + (\lambda-\xi^2) \psi = 0

设解的形式为

ψ=eξ22u(ξ) \psi = \mathrm{e}^{-\frac{\xi^2}{2}} u(\xi)

之所以这么设,可以按如下方式考虑:当 ξ\xi\to\infty 时,方程近似表示为 d2dξ2ψξ2ψ=0\frac{\mathrm{d}^2}{\mathrm{d}\xi^2} \psi -\xi^2 \psi = 0 ,当 ψ=e±ξ22\psi = \mathrm{e}^{\pm\frac{\xi^2}{2}} 时, ψ=±ξe±ξ22\psi' = \pm\xi \mathrm{e}^{\pm\frac{\xi^2}{2}}ψ=(ξ2±1)e±ξ22ξ2e±ξ22=ξ2ψ\psi'' = (\xi^2\pm1) \mathrm{e}^{\pm\frac{\xi^2}{2}} \approx \xi^2 \mathrm{e}^{\pm\frac{\xi^2}{2}} = \xi^2\psi ,故方程的近似解为 ψe±ξ22\psi \sim \mathrm{e}^{\pm\frac{\xi^2}{2}} ,而根据束缚态边界条件,即 ξ\xi\to\inftyψ0\psi\to0 ,应舍去 ψeξ22\psi \sim \mathrm{e}^{\frac{\xi^2}{2}}

将上述解的形式代入原方程,可得到 u(ξ)u(\xi) 满足的方程

d2dξ2u+2ξddξu+(λ1)u=0 \frac{\mathrm{d}^2}{\mathrm{d}\xi^2} u + 2\xi \frac{\mathrm{d}}{\mathrm{d}\xi} u + (\lambda-1) u = 0

此即Hermite方程,可以通过级数解法求解:在 ξ=0\xi=0 附近,用幂级数展开

u(ξ)=k=0+ckξk u(\xi) = \sum_{k=0}^{+\infty} c_k \xi^k

代入Hermite方程,比较同幂项的系数,可得

ck+2=2k(λ1)(k+2)(k+1)ck(k=0,1,2,) c_{k+2} = \frac{2k-(\lambda-1)}{(k+2)(k+1)} c_k \kern 2em (k=0,1,2,\cdots)

故所有的偶次项系数都可以用 c0c_0 来表示,所有的奇次项系数都可以用 c1c_1 来表示,把 c0c_0c1c_1 作为两个任意常数,就可以得到Hermite方程两个线性无关的解,即级数的偶次项部分与奇次项部分

u1(ξ)=m=0+c2mξ2m=c0+c2ξ2+c4ξ4+u2(ξ)=m=0+c2m+1ξ2m+1=c1ξ+c3ξ3+c5ξ5+ u_1(\xi) = \sum_{m=0}^{+\infty} c_{2m} \xi^{2m} = c_0 + c_2\xi^2 + c_4\xi^4 + \cdots \ \kern 1em \ u_2(\xi) = \sum_{m=0}^{+\infty} c_{2m+1} \xi^{2m+1} = c_1\xi + c_3\xi^3 + c_5\xi^5 + \cdots

考虑当 ξ\xi\to\infty 时的情况,当 k+k\to+\infty 时,

ck+2ck=2k(λ1)(k+2)(k+1)2k \frac{c_{k+2}}{c_k} = \frac{2k-(\lambda-1)}{(k+2)(k+1)} \to \frac{2}{k}

对于偶数的情况,即 k=2mk=2m ,有 c2m+2/c2m1/mc_{2m+2}/c_{2m} \sim 1/m ,这与 eξ2\mathrm{e}^{\xi^2} 的Taylor展开

eξ2=m=0+ξ2mm! \mathrm{e}^{\xi^2} = \sum_{m=0}^{+\infty} \frac{\xi^{2m}}{m!}

相邻两项的系数比相同,因此,

u1(ξ)eξ2 u_1(\xi) \sim \mathrm{e}^{\xi^2}

同理可得

u2(ξ)ξeξ2 u_2(\xi) \sim \xi\mathrm{e}^{\xi^2}

代回到波函数可得

ψ1=eξ22u1(ξ)eξ22ψ2=eξ22u2(ξ)ξeξ22 \psi_1 = \mathrm{e}^{-\frac{\xi^2}{2}} u_1(\xi) \sim \mathrm{e}^{\frac{\xi^2}{2}} \ \kern 1em \ \psi_2 = \mathrm{e}^{-\frac{\xi^2}{2}} u_2(\xi) \sim \xi\mathrm{e}^{\frac{\xi^2}{2}}

这不满足束缚态的边界条件(当 ξ\xi\to\inftyψ0\psi\to0 ),故 u1u_1u2u_2 两个无穷级数解中,必须至少有一个中断为多项式,也就是要找到合适的 λ\lambda ,使得存在 kNk\in\mathbb{N} 满足 2k(λ1)(k+2)(k+1)=0\frac{2k-(\lambda-1)}{(k+2)(k+1)} = 0 ,故当

λ1=2n(n=0,1,2,) \lambda-1 = 2n \kern 2em (n=0,1,2,\cdots)

时,级数将中断一个多项式( cn+2=cn+4=cn+6==0c_{n+2} = c_{n+4} = c_{n+6} = \cdots = 0 )。当 nn 为偶时, u1u_1 中断为Hermite多项式 Hn(ξ)\mathrm{H}_n(\xi)u2u_2 仍为无穷级数;当 nn 为奇时, u2u_2 中断为Hermite多项式 Hn(ξ)\mathrm{H}_n(\xi)u1u_1 仍为无穷级数。其中Hermite多项式表示为

Hn(ξ)=(1)neξ2dndξneξ2=(2ξ)n+n(n1)(2ξ)n2++(1)n2n!n2!(2ξ)n2n2 \mathrm{H}_n(\xi) = (-1)^n \mathrm{e}^{\xi^2} \frac{\mathrm{d}^n }{\mathrm{d}\xi^n} \mathrm{e}^{-\xi^2} \ = (2\xi)^n + n(n-1)(2\xi)^{n-2} + \cdots + (-1)^{\left\lfloor \frac{n}{2} \right\rfloor} \frac{n!}{\left\lfloor \frac{n}{2} \right\rfloor!} (2\xi)^{n-2\left\lfloor \frac{n}{2} \right\rfloor}

例如

H0(ξ)=1H1(ξ)=2ξH2(ξ)=4ξ22 \mathrm{H}_0(\xi) = 1 \ \mathrm{H}_1(\xi) = 2\xi \ \mathrm{H}_2(\xi) = 4\xi^2 -2

Hermite多项式的带权正交归一性表示为

+Hm(ξ)Hn(ξ)eξ2dξ=π2nn!δmn \int_{-\infty}^{+\infty} \mathrm{H}_m(\xi) \mathrm{H}n(\xi) \mathrm{e}^{-\xi^2} \mathrm{d}\xi = \sqrt{\pi} 2^n \cdot n! \delta{mn}

根据 λ\lambda 满足的离散化条件,可以求出一维谐振子的能量本征值

λ=E12ω=2n+1En=(n+12)ω(n=0,1,2,) \lambda = \frac{E}{\frac12\hbar\omega} = 2n+1 \Longrightarrow E_n = (n+\frac12)\hbar\omega \kern 2em (n=0,1,2,\cdots)

借助Hermite多项式,并把 ξ=αx\xi=\alpha x 代入,可以表示出一维谐振子的能量本征函数

ψneξ22Hn(ξ)ψn(x)=Aneα2x22Hn(αx) \psi_n \propto \mathrm{e}^{-\frac{\xi^2}{2}} \mathrm{H}_n(\xi) \Longrightarrow \psi_n(x) = A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x)

根据Hermite多项式的带权正交归一性,

(ψm,ψn)=AmAn+Hm(ξ)Hn(ξ)eξ2dξ=AmAnπ2nn!δmn (\psi_m,\psi_n) = A_mA_n \int_{-\infty}^{+\infty} \mathrm{H}_m(\xi) \mathrm{H}n(\xi) \mathrm{e}^{-\xi^2} \mathrm{d}\xi = A_mA_n\sqrt{\pi} 2^n \cdot n! \delta{mn}

可得归一化系数

An=απ2nn! A_n = \sqrt{\frac{\alpha}{\sqrt{\pi}2^n\cdot n!}}

这样波函数就满足了正交归一化条件

(ψm,ψn)=δmn (\psi_m,\psi_n) = \delta_{mn}

讨论
能量本征值

En=(n+12)ω(n=0,1,2,) E_n = (n+\frac12)\hbar\omega \kern 2em (n=0,1,2,\cdots)

一维谐振子的能量是均匀分布的,相邻的两条能级间距为 En+1En=ωE_{n+1} - E_n = \hbar\omega

波函数与能级宇称

最低的三条能级上的谐振子波函数如下:

ψ0(x)=απ14eα2x22ψ1(x)=2απ14αx eα2x22ψ2(x)=1π14α2(2α2x21)eα2x22 \psi_0(x) = \frac{\sqrt{\alpha}}{\pi^{\frac{1}{4}}} \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \ \kern 1em \ \psi_1(x) = \frac{\sqrt{2\alpha}}{\pi^{\frac{1}{4}}} \alpha x\ \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \ \kern 1em \ \psi_2(x) = \frac{1}{\pi^{\frac{1}{4}}} \sqrt{\frac{\alpha}{2}} (2\alpha^2x^2-1) \mathrm{e}^{-\frac{\alpha^2x^2}{2}}

2.2一维谐振子较低几条能级的波函数

2.2一维谐振子较低几条能级的概率位置密度

其中 ψn(x)\psi_n(x)nn 个节点。

由于一维谐振子势具有空间反射不变性( V(x)=V(x)V(-x)=V(x) ),根据定理3推论, ψn(x)\psi_n(x) 必有确定的宇称,事实上,可以证明

ψn(x)=(1)nψn(x) \psi_n(-x) = (-1)^n \psi_n(x)

能级的宇称偶奇相间,基态是偶宇称。

基态能量与概率分布

一维谐振子基态能量为

E0=12ω E_0 = \frac12 \hbar\omega

其并不为零(可以用不确定性关系解释),称为零点能

处于基态的谐振子在空间的概率分布为

ψ0(x)2=απeα2x2 |\psi_0(x)|^2 = \frac{\alpha}{\sqrt{\pi}} \mathrm{e}^{-\alpha^2x^2}

这是一个Gauss型分布,在原点 (x=0)(x=0) 处找到粒子的概率最大。

补充(二维、三维谐振子)

对于一维谐振子,其Hamilton算符

H^=22m2x2+12mω2x2 \hat{H} = -\frac{\hbar^2}{2m} \frac{\partial^2}{\partial x^2} + \frac12m\omega^2x^2

二维谐振子

二维谐振子的势能可以表示为

V(x,y)=12mω2r2=12mω2(x2+y2) V(x,y) = \frac12m\omega^2r^2 = \frac12m\omega^2(x^2+y^2)

其Hamilton算符可以表示为

H^=22m(2x2+2y2)+12mω2(x2+y2)=H^x+H^y \hat{H} = -\frac{\hbar^2}{2m} \left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right) + \frac12m\omega^2(x^2+y^2) = \hat{H}_x + \hat{H}_y

对于二维谐振子的能量本征方程

H^ψ(x,y)=Eψ(x,y) \hat{H} \psi(x,y) = E \psi(x,y)

由于 x,yx,y 相独立,可以使用分离变量法求解,令 ψ(x,y)=ψx(x)ψy(y)\psi(x,y)=\psi_x(x)\psi_y(y) ,则能量本征方程可表示为

(H^x+H^y)ψxψy=EψxψyψyH^xψx+ψxH^yψy=EψxψyH^xψxψx+H^yψyψy=E (\hat{H}_x + \hat{H}_y)\psi_x\psi_y = E\psi_x\psi_y \ \Downarrow \ \psi_y\hat{H}_x\psi_x + \psi_x\hat{H}_y\psi_y = E\psi_x\psi_y \ \Downarrow \ \frac{\hat{H}_x\psi_x}{\psi_x} + \frac{\hat{H}_y\psi_y}{\psi_y} = E

这样,能量本征方程就可以分离为 x,yx,y 两个方向上的方程:

H^xψx=Exψx ,H^yψy=Eyψy \hat{H}_x\psi_x = E_x\psi_x\ , \kern 1em \hat{H}_y\psi_y = E_y\psi_y

则二维谐振子的能量本征函数

ψnxny(x,y)=ψnx(x)ψny(y)(nx,ny=0,1,2,) \psi_{n_xn_y}(x,y) = \psi_{n_x}(x) \psi_{n_y}(y) \kern 2em (n_x,n_y=0,1,2,\cdots)

其中 ψnx,ψny\psi_{n_x},\psi_{n_y} 与一维谐振子的 ψn\psi_n 函数相同。

二维谐振子的能量本征值

Enxny=Enx+Eny=(12+nx)ω+(12+ny)ω=(1+nx+ny)ω(nx,ny=0,1,2,) E_{n_xn_y} = E_{n_x} + E_{n_y} = (\frac12+n_x)\hbar\omega + (\frac12+n_y)\hbar\omega = (1+n_x+n_y)\hbar\omega \kern 2em (n_x,n_y=0,1,2,\cdots)

其中 Enx,EnyE_{n_x},E_{n_y} 与一维谐振子的 EnE_n 表达式相同,记 N=nx+nyN=n_x+n_y ,则

Enxny=(1+N)ω(N=0,1,2,) E_{n_xn_y} = (1+N)\hbar\omega \kern 2em (N=0,1,2,\cdots)

对于给定的 NN(nx,ny)(n_x,n_y) 的可能取值共有 N+1N+1 种(即 (0,N),(1,N1),,(N,0)(0,N),(1,N-1),\cdots,(N,0) ),故能级简并度

fN=N+1 f_N = N+1

三维谐振子

三维谐振子的结论与二维谐振子类似,能量本征函数

ψnxnyny(x,y,z)=ψnx(x)ψny(y)ψnz(z)(nx,ny,nz=0,1,2,) \psi_{n_xn_yn_y}(x,y,z) = \psi_{n_x}(x) \psi_{n_y}(y) \psi_{n_z}(z) \kern 2em (n_x,n_y,n_z=0,1,2,\cdots)

能量本征值

Enxnyny=Enx+Eny+Enz=(12+nx)ω+(12+ny)ω+(12+nz)ω=(32+nx+ny+nz)ω(nx,ny,ny=0,1,2,) E_{n_xn_yn_y} = E_{n_x} + E_{n_y} + E_{n_z} = (\frac12+n_x)\hbar\omega + (\frac12+n_y)\hbar\omega + (\frac12+n_z)\hbar\omega \ = (\frac32+n_x+n_y+n_z)\hbar\omega\kern 2em (n_x,n_y,n_y=0,1,2,\cdots)

N=nx+ny+nzN=n_x+n_y+n_z ,则

Enxnynz=(32+N)ω(N=0,1,2,) E_{n_xn_yn_z} = (\frac32+N)\hbar\omega \kern 2em (N=0,1,2,\cdots)

对于给定的 NN(nx,ny,nz)(n_x,n_y,n_z) 的可能取值共有

nx=0N(N+1nx)=k=1N+1k=12(N+1)(N+2) \sum_{n_x=0}^{N} (N+1-n_x) = \sum_{k=1}^{N+1} k = \frac12(N+1)(N+2)

故能级简并度

fN=12(N+1)(N+2) f_N = \frac12(N+1)(N+2)

δ\delta 势阱

δ\delta 势阱表示为

V(x)=γδ(x)(γ>0) V(x) = -\gamma \delta(x) \kern 2em (\gamma > 0)

质量为 mm 的粒子在 δ\delta 势阱中运动:在 x0x\ne0 处有 V(x)=0V(x)=0 ,所以 E>0E>0 为游离态, EE 可以取一切正实数值,是连续变化的;而 E<0E<0 时,则可能存在束缚能量本征态, EE 只能取离散值。以下讨论束缚态,即 E<0E<0 的情况。

能量本征方程为

ψ(x)=2m2[E+γδ(x)]ψ(x) \psi''(x) = -\frac{2m}{\hbar^2} \left[E+\gamma\delta(x)\right] \psi(x)

左右两边同时积分可以得到 δ\delta 势阱中 ψ\psi'跃变条件

limε0+εεψ(x)dx=limε0+εε2m2[E+γδ(x)]ψ(x)dxψ(0+)ψ(0)=2mγ2ψ(0) \lim_{\varepsilon\to0^+} \int_{-\varepsilon}^{\varepsilon} \psi''(x) \mathrm{d}x = \lim_{\varepsilon\to0^+} \int_{-\varepsilon}^{\varepsilon} -\frac{2m}{\hbar^2} \left[E+\gamma\delta(x)\right] \psi(x) \mathrm{d}x \ \Downarrow \ \psi'(0^+) - \psi'(0^-) = -\frac{2m\gamma}{\hbar^2} \psi(0)

β=2mE>0 \beta = \frac{\sqrt{-2mE}}{\hbar} > 0

则在 x0x\ne0 的区域,能量本征方程可化为

ψ(x)β2ψ(x)=0 \psi''(x) - \beta^2\psi(x) = 0

解得

ψ(x)=Aeβx+Beβx \psi(x) = A\mathrm{e}^{\beta x} + B\mathrm{e}^{-\beta x}

考虑到束缚态边界条件,即在 xx\to\infty 处,要求 ψ(x)0\psi(x)\to0 ,则

ψ(x)={Aeβx,x<0Beβx,x>0 \psi(x) = \begin{cases} A\mathrm{e}^{\beta x}, & x<0 \ B\mathrm{e}^{-\beta x}, & x>0 \end{cases}

考虑到势阱具有空间反射不变性 V(x)=V(x)V(-x)=V(x) ,由定理3推论可知,束缚态能量本征函数(由定理7知其不简并)必具有确定的宇称,以下分别讨论:

对于偶宇称态

波函数应表示为

ψ(x)={Aeβx,x<0Aeβx,x>0 \psi(x) = \begin{cases} A\mathrm{e}^{\beta x}, & x<0 \ A\mathrm{e}^{-\beta x}, & x>0 \end{cases}

按照 ψ\psi' 跃变条件

ψ(0+)ψ(0)=2mγ2ψ(0)AβAβ=2mγ2A \psi'(0^+) - \psi'(0^-) = -\frac{2m\gamma}{\hbar^2} \psi(0) \ \Downarrow \ -A\beta - A\beta = -\frac{2m\gamma}{\hbar^2}A

可得

β=mγ2 \beta = \frac{m\gamma}{\hbar^2}

则可得出粒子的能量本征值

β=mγ2=2mEE=mγ222 \beta = \frac{m\gamma}{\hbar^2} = \frac{\sqrt{-2mE}}{\hbar} \Longrightarrow E = -\frac{m\gamma^2}{2\hbar^2}

由归一化条件可得

(ψ,ψ)=20+A2e2βxdx=A2β=1 (\psi,\psi) = 2\int_{0}^{+\infty} |A|^2\mathrm{e}^{-2\beta x} \mathrm{d}x = \frac{|A|^2}{\beta} = 1

δ\delta 势的特征长度

L=1β=2mγ L = \frac{1}{\beta} =\frac{\hbar^2}{m\gamma}

A=β=1L |A| = \sqrt{\beta} = \frac{1}{\sqrt{L}}

这样归一化的束缚能量本征态波函数可表示为

ψ(x)=1LexL \psi(x) = \frac{1}{\sqrt{L}} \mathrm{e}^{-\frac{|x|}{L}}

势阱归一化的束缚能量本征态波函数

对于奇宇称态

波函数应表示为

ψ(x)={Aeβx,x<0Aeβx,x>0 \psi(x) = \begin{cases} A\mathrm{e}^{\beta x}, & x<0 \ -A\mathrm{e}^{-\beta x}, & x>0 \end{cases}

由波函数在 x=0x=0 点连续,可以得到

ψ(0)=ψ(0+)A=AA=0 \psi(0^-) = \psi(0^+) \Longrightarrow A = -A \Longrightarrow A = 0

所以不可能存在奇宇称束缚能量本征态。

从物理上考虑,奇宇称波函数在 x=0x=0 点必为零,而 δ\delta 势又恰好只在 x=0x=0 点其作用,所以 δ\delta 势阱对奇宇称态没有影响,因而不可能形成束缚态。

2.3 散射问题:量子隧穿(隧道)效应

方势垒的反射与透射

模型描述与结论

设具有一定能量 EE 的质量为 mm 的粒子沿 xx 轴正方向射向方势垒

V(x)={V0 ,0<x<a0 ,x<0,x>a V(x) = \begin{cases} V_0\ , & 0<x<a \ 0\ , & x<0,x>a \end{cases}

无论粒子能量 E>V0E>V_0 还是 E<V0E<V_0 ,都有一定概率穿透势垒,也有一定概率被反射回去。主要考虑 0<E<V00<E<V_0 的情况,令

k=2mE,κ=2m(V0E) k = \frac{\sqrt{2mE}}{\hbar}, \kern 1em \kappa = \frac{\sqrt{2m(V_0-E)}}{\hbar}

则波函数为

ψ(x)={eikx+Reikx ,x<0Aeκx+Beκx ,0<x<aSeikx ,x>a \psi(x) = \begin{cases} \mathrm{e}^{\mathrm{i}kx}+R\mathrm{e}^{-\mathrm{i}kx}\ , & x<0 \ A\mathrm{e}^{\kappa x}+B\mathrm{e}^{-\kappa x}\ , & 0<x<a \ S\mathrm{e}^{\mathrm{i}kx}\ , & x>a \end{cases}

其中 ReikxR\mathrm{e}^{-\mathrm{i}kx} 为反射波, SeikxS\mathrm{e}^{\mathrm{i}kx} 为透射波

隧道效应

透射系数

T=S2=4k2κ2(k2+κ2)2sinh2(κa)+4k2κ2=[1+14EV0(1EV0)sinh2(κa)]1 T = |S|^2 = \frac{4k^2\kappa^2}{(k^2+\kappa^2)^2\sinh^2(\kappa a)+4k^2\kappa^2} = \left[1+\frac{1}{\frac{4E}{V_0}(1-\frac{E}{V_0})}\sinh^2(\kappa a)\right]^{-1}

反射系数

R2=(k2+κ2)2sinh2(κa)(k2+κ2)2sinh2(κa)+4k2κ2 |R|^2 = \frac{(k^2+\kappa^2)^2\sinh^2(\kappa a)}{(k^2+\kappa^2)^2\sinh^2(\kappa a)+4k^2\kappa^2}

模型求解

在势垒外 (x<0,x>a)(x<0,x>a) ,能量本征方程表示为

ψ(x)+2mE2ψ(x)=0 \psi''(x) + \frac{2mE}{\hbar^2}\psi(x) = 0

k=2mEk=\frac{\sqrt{2mE}}{\hbar} ,该方程的两个线性无关解可取为 ψ(x)e±ikx\psi(x) \sim \mathrm{e}^{\pm\mathrm{i}kx} 。粒子是从左入射,由于势垒的存在,在 x<ax<a 的区域中,既有入射波 eikx\mathrm{e}^{\mathrm{i}kx} ,也有反射波 eikx\mathrm{e}^{-\mathrm{i}kx} ;而在 x>ax>a 的区域中,则只有透射波 eikx\mathrm{e}^{\mathrm{i}kx} ,所以

ψ(x)={eikx+Reikx ,x<0Seikx ,x>a \psi(x) = \begin{cases} \mathrm{e}^{\mathrm{i}kx}+R\mathrm{e}^{-\mathrm{i}kx}\ , & x<0 \ S\mathrm{e}^{\mathrm{i}kx}\ , & x>a \end{cases}

这里把入射波的波幅任意地取为 11 ,只是为了方便求解,由于还没有归一化,只要相对比例一定,对透射和反射系数都没有影响。

在势垒内部 (0<x<a)(0<x<a) ,能量本征方程表示为

ψ(x)2m(V0E)2ψ(x)=0 \psi''(x) - \frac{2m(V_0-E)}{\hbar^2}\psi(x) = 0

κ=2m(V0E)\kappa=\frac{\sqrt{2m(V_0-E)}}{\hbar} ,解得

ψ(x)=Aeκx+Beκx(0<x<a) \psi(x) = A\mathrm{e}^{\kappa x}+B\mathrm{e}^{-\kappa x} \kern 2em (0<x<a)

根据 ψ\psiψ\psi' 分别在 x=0x=0x=ax=a 处连续,可以得到如下关于 R,S,A,BR,S,A,B 的方程组

{1+R=A+Bik(1R)=κ(AB)Seika=Aeκa+BeκaikSeikx=κ(AeκaBeκa) \left{\begin{matrix} 1+R = A+B \ \mathrm{i}k(1-R) = \kappa(A-B) \ S\mathrm{e}^{\mathrm{i}ka} = A\mathrm{e}^{\kappa a}+B\mathrm{e}^{-\kappa a} \ \mathrm{i}kS\mathrm{e}^{\mathrm{i}kx} = \kappa(A\mathrm{e}^{\kappa a}-B\mathrm{e}^{-\kappa a}) \end{matrix}\right.

为了求解该方程组,可由前两个方程用 RR 表示 A,BA,B ,再由后两个方程用 SS 表示 A,BA,B ,两种表示对比可得到关于 S,RS,R 的方程组,进一步求出 SSRR ,回代得到 A,BA,B ;或者使用线性代数的知识求解也可。完整的解较为复杂,这里不再展示。

入射的粒子流密度为

ji=i2m(ψiddxψiψiddxψi)=i2m(eikxddxeikxeikxddxeikx)=km=v j_i = -\frac{\mathrm{i}\hbar}{2m} (\psi_i^* \frac{\mathrm{d}}{\mathrm{d}x} \psi_i - \psi_i \frac{\mathrm{d}}{\mathrm{d}x} \psi_i^*) = -\frac{\mathrm{i}\hbar}{2m} (\mathrm{e}^{-\mathrm{i}kx} \frac{\mathrm{d}}{\mathrm{d}x} \mathrm{e}^{\mathrm{i}kx} - \mathrm{e}^{\mathrm{i}kx} \frac{\mathrm{d}}{\mathrm{d}x} \mathrm{e}^{-\mathrm{i}kx}) = \frac{\hbar k}{m} = v

类似的,可以计算出反射流密度 jrj_r 和透射流密度 jtj_t 分别为

jr=R2v,jt=S2v j_r = -|R|^2v , \kern 1em j_t = |S|^2v

所以

反射系数=jrji=R2透射系数=jtji=S2 反射系数 = \frac{|j_r|}{|j_i|} = |R|^2 \ \kern 1em \ 透射系数 = \frac{|j_t|}{|j_i|} = |S|^2

代入求解方程组得到的 R,SR,S ,得

透射系数

T=S2=4k2κ2(k2+κ2)2sinh2(κa)+4k2κ2 T = |S|^2 = \frac{4k^2\kappa^2}{(k^2+\kappa^2)^2\sinh^2(\kappa a)+4k^2\kappa^2}

反射系数

R2=(k2+κ2)2sinh2(κa)(k2+κ2)2sinh2(κa)+4k2κ2 |R|^2 = \frac{(k^2+\kappa^2)^2\sinh^2(\kappa a)}{(k^2+\kappa^2)^2\sinh^2(\kappa a)+4k^2\kappa^2}

讨论
概率守恒

R2+S2=1 |R|^2 + |S|^2 = 1

隧道效应

通过整理,透射系数还可以表示为

T=[1+14EV0(1EV0)sinh2(κa)]1 T = \left[1+\frac{1}{\frac{4E}{V_0}(1-\frac{E}{V_0})}\sinh^2(\kappa a)\right]^{-1}

0<E<V00<E<V_0 时,透射系数 T0T\ne0 ,这种粒子能穿透比它动能更高的势垒的现象,称为量子隧穿效应(或称隧道效应(tunnel effect)、势垒贯穿),它是粒子具有波动性的表现。当然,这种现象一般概率较低,只有在一定的条件下才比较显著。

近似公式

κa1\kappa a\gg1 ,则 sinh(κa)=12(eκaeκa)12eκa1\sinh(\kappa a) = \frac12(\mathrm{e}^{\kappa a}-\mathrm{e}^{-\kappa a}) \approx \frac12\mathrm{e}^{\kappa a} \gg 1 ,则透射系数可近似表示为

T4k2κ2(k2+κ2)2sinh2(κa)4k2κ2(k2+κ2)2(12eκa)2=16k2κ2(k2+κ2)2e2κa=16E(V0E)V02e2a2m(V0E) T \approx \frac{4k^2\kappa^2}{(k^2+\kappa^2)^2\sinh^2(\kappa a)} \approx \frac{4k^2\kappa^2}{(k^2+\kappa^2)^2(\frac12\mathrm{e}^{\kappa a})^2} = \frac{16k^2\kappa^2}{(k^2+\kappa^2)^2}\mathrm{e}^{-2\kappa a} = \frac{16E(V_0-E)}{V_0^2}\mathrm{e}^{-\frac{2a}{\hbar}\sqrt{2m(V_0-E)}}

若记

T0=16EV0(1EV0) T_0 = 16 \frac{E}{V_0} \left(1-\frac{E}{V_0}\right)

TT0exp(2a2m(V0E)) T \approx T_0 \exp\left(-\frac{2a}{\hbar}\sqrt{2m(V_0-E)}\right)

可以看出 TT 灵敏地依赖于粒子的质量 mm 、势垒宽度 aa 以及 (V0E)(V_0-E)

对于一般形状的势垒,可以将其视为许多方势垒相邻排布,若透射系数 T1T\ll1 ,则对于在 axba\le x\le b 之间的势垒,有 WKB准经典近似公式

TT0exp{2ab2m[V(x)E] dx} T \approx T_0 \exp\left{ -\frac{2}{\hbar} \int_a^b \sqrt{2m[V(x)-E]}\ \mathrm{d}x \right}

方势阱的反射、透射与共振

首先考虑方势垒中 E>V0E>V_0 的情况,令

k=2m(EV0) k' = \frac{\sqrt{2m(E-V_0)}}{\hbar}

只需要将 κik\kappa \longrightarrow \mathrm{i}k' ,可得透射系数

T=4k2k2(k2k2)2sin2(ka)+4k2k2=[1+14(kkkk)2sin2(ka)]1 T = \frac{4k^2k'^2}{(k^2-k'^2)^2\sin^2(k'a)+4k^2k'^2} = \left[1+\frac14 \left(\frac{k}{k'}-\frac{k'}{k}\right)^2 \sin^2(k'a)\right]^{-1}

ka=nπk'a=n\pi 时, sin(ka)=0\sin(k'a)=0 ,故 T=1T=1 ,称为共振透射

而对于方势阱的透射,上述理论仍然适用,只需要把 V0V0V_0 \longrightarrow -V_0 ,则相应的

k=2m(E+V0)2mE=k k' = \frac{\sqrt{2m(E+V_0)}}{\hbar} \ge \frac{\sqrt{2mE}}{\hbar} = k

此时透射系数

T=[1+14(kkkk)2sin2(ka)]1=[1+sin2(ka)4EV0(1+EV0)]1 T = \left[1+\frac14 \left(\frac{k}{k'}-\frac{k'}{k}\right)^2 \sin^2(k'a)\right]^{-1} = \left[1 + \frac{\sin^2(k'a)}{4\frac{E}{V_0}\left(1+\frac{E}{V_0}\right)}\right]^{-1}

可以看出,若 V0=0V_0=0 ,则 T=1T=1 ;若 V00V_0\ne0 ,则一般情况下 T<1,R20T<1,|R|^2\ne0 ,即粒子有一定概率被势阱弹回。

对于给定势阱,透射系数 TT 完全依赖于入射粒子的能量 EE ,透射系数 T(E)T(E)EE 的变化如图所示

方势阱T(E)

如果 EV0E \ll V_0 ,则一般来说 TT 值很小,除非入射粒子的能量 EE 合适,使 sin(ka)=0\sin(k'a)=0 ,此时 T=1T=1 (反射系数 R2=0|R|^2=0 ),这种现象被称为共振透射,它出现的条件是

ka=nπ(n=1,2,3,) k'a = n\pi \kern 2em (n=1,2,3,\cdots)

可以得到共振能级 EnE_n 的表达式为

k=2m(E+V0)=nπaEn=V0+n2π222ma2 k' = \frac{\sqrt{2m(E+V_0)}}{\hbar} = \frac{n\pi}{a} \Longrightarrow E_n = -V_0 + \frac{n^2\pi^2\hbar^2}{2ma^2}

与此相对,反射最强的条件是

ka=(n+12)π(n=0,1,2,) k'a = (n+\frac12)\pi \kern 2em (n=0,1,2,\cdots)

δ\delta 势的穿透

模型建立与求解

设具有一定能量 EE 的质量为 mm 的粒子沿 xx 轴正方向射向 δ\delta 势垒

V(x)=γδ(x)(γ>0) V(x) = \gamma\delta(x) \kern 2em (\gamma>0)

能量本征方程为

ψ(x)=2m2[Eγδ(x)]ψ(x) \psi''(x) = -\frac{2m}{\hbar^2} \left[E-\gamma\delta(x)\right] \psi(x)

左右两边同时积分可以得到 δ\delta 势阱中 ψ\psi'跃变条件

limε0+εεψ(x)dx=limε0+εε2m2[Eγδ(x)]ψ(x)dxψ(0+)ψ(0)=2mγ2ψ(0) \lim_{\varepsilon\to0^+} \int_{-\varepsilon}^{\varepsilon} \psi''(x) \mathrm{d}x = \lim_{\varepsilon\to0^+} \int_{-\varepsilon}^{\varepsilon} -\frac{2m}{\hbar^2} \left[E-\gamma\delta(x)\right] \psi(x) \mathrm{d}x \ \Downarrow \ \psi'(0^+) - \psi'(0^-) = \frac{2m\gamma}{\hbar^2} \psi(0)

k=2mE>0 k = \frac{\sqrt{2mE}}{\hbar} > 0

则在 x0x\ne0 的区域,能量本征方程可化为

ψ(x)+k2ψ(x)=0 \psi''(x) + k^2\psi(x) = 0

该方程的两个线性无关解可取为 ψ(x)e±ikx\psi(x) \sim \mathrm{e}^{\pm\mathrm{i}kx} 。粒子是从左入射,由于势垒的存在,在 x<0x<0 的区域中,既有入射波 eikx\mathrm{e}^{\mathrm{i}kx} ,也有反射波 eikx\mathrm{e}^{-\mathrm{i}kx} ;而在 x>0x>0 的区域中,则只有透射波 eikx\mathrm{e}^{\mathrm{i}kx} ,所以解得

ψ(x)={eikx+Reikx ,x<0Seikx ,x>0 \psi(x) = \begin{cases} \mathrm{e}^{\mathrm{i}kx}+R\mathrm{e}^{-\mathrm{i}kx}\ , & x<0 \ S\mathrm{e}^{\mathrm{i}kx}\ , & x>0 \end{cases}

根据在 x=0x=0ψ\psi 连续与 ψ\psi' 的跃变条件,可以得到如下关于 R,SR,S 的方程组

{1+R=SikSik(1R)=2mγ2S \left{\begin{matrix} 1+R = S \ \mathrm{i}kS - \mathrm{i}k(1-R) = \frac{2m\gamma}{\hbar^2}S \end{matrix}\right.

解得

S=11+imγ2kR=imγ2k1+imγ2k S = \frac{1}{1+\frac{\mathrm{i}m\gamma}{\hbar^2k}} \ \kern 1em \ R = \frac{-\frac{\mathrm{i}m\gamma}{\hbar^2k}}{1+\frac{\mathrm{i}m\gamma}{\hbar^2k}}

透射系数

T=S2=11+m2γ24k2=11+mγ222E T = |S|^2 = \frac{1}{1+\frac{m^2\gamma^2}{\hbar^4k^2}} = \frac{1}{1+\frac{m\gamma^2}{2\hbar^2E}}

反射系数

R2=m2γ24k21+m2γ24k2=mγ222E1+mγ222E |R|^2 = \frac{\frac{m^2\gamma^2}{\hbar^4k^2}}{1+\frac{m^2\gamma^2}{\hbar^4k^2}} = \frac{\frac{m\gamma^2}{2\hbar^2E}}{1+\frac{m\gamma^2}{2\hbar^2E}}

讨论
概率守恒

R2+S2=1 |R|^2 + |S|^2 = 1

δ\delta 势阱的穿透

如果把 δ\delta 势垒换为 δ\delta 势阱( γγ\gamma\longrightarrow-\gamma ),透射系数与反射系数的值均不变。

特征长度与特征能量

δ\delta 势的特征长度 L=2mγL=\frac{\hbar^2}{m\gamma} ,特征能量为 mγ22\frac{m\gamma^2}{\hbar^2}

透射波的波幅 SS 只依赖于 mγ2k=1k/2mγ\frac{m\gamma}{\hbar^2k} = \frac{1}{k} / \frac{\hbar^2}{m\gamma} ,即入射粒子波长与 δ\delta 势特征长度之比;而透射系数 TT 只依赖于 mγ22E=mγ22/E\frac{m\gamma^2}{\hbar^2E} = \frac{m\gamma^2}{\hbar^2} / E ,即特征能量与入射粒子能量之比。当 Emγ22E \gg \frac{m\gamma^2}{\hbar^2} 时, T1T\approx1 ,即高能极限下粒子将完全穿透 δ\delta 势垒。

第3章 力学量用算符表达

3.1 算符的运算

算符基本介绍

算符是作用于波函数把它变成另一个函数的运算符号,代表力学量 AA 的算符记做 A^\hat{A} 。量子力学中任一可观测力学量 AA 可以用线性Hermite算符 A^\hat{A} 来表示,这些算符作用于态的波函数。

为强调算符的特点,常常在算符的符号上加一个 “^\hat{\kern 1em}” 号,但在不会引起误解的地方,也常把 “^\hat{\kern 1em}” 略去。

在全部量子理论中,时间一直保持为连续变化的参量,不存在相应的“时间算符”

此小节之后搬运了第一章的内容,主要是从平均值的角度引出算符。

平均值假定

粒子处于波函数 ψ(r)\psi(\vec{r}) 所描述的状态下,力学量(又叫可观测量)都有确定的概率分布,因而有确定的平均值(又叫期待值)。在任意状态 ψ\psi 上,对力学量 AA 进行足够多次的测量,所得结果的平均值为

Aˉ=(ψ,A^ψ)(ψ,ψ) \bar{A} = \frac{(\psi,\hat{A}\psi)}{(\psi,\psi)}

其中 A^\hat{A} 是力学量 AA 对应的算符,若波函数已归一化,则

Aˉ=(ψ,A^ψ) \bar{A} = (\psi,\hat{A}\psi)

坐标算符与动量算符

在波函数 ψ\psi 已归一化的条件下,位置 xx 的平均值为

xˉ=+ψ(r)2xd3r=+ψ(r)xψ(r)d3r \bar{x} = \int_{-\infty}^{+\infty} \left|\psi(\vec{r})\right|^2, x, \mathrm{d}^3r= \int_{-\infty}^{+\infty} \psi^*(\vec{r}), x, \psi(\vec{r}), \mathrm{d}^3r

可以得到坐标表象下的坐标算符为

x^=x \hat{x} = x

同理

y^=y,z^=z,r^=r \hat{y} = y,\kern 12pt\hat{z} = z,\kern 12pt\hat{\vec{r}} = \vec{r}

如果状态用动量表象波函数 φ(p)\varphi(\vec{p}) 来表示,则粒子动量的平均值为

pˉ=+φ(p)2pd3p=+φ(p)pφ(p)d3p \bar{\vec{p}} = \int_{-\infty}^{+\infty} \left|\varphi(\vec{p})\right|^2, \vec{p}, \mathrm{d}^3p= \int_{-\infty}^{+\infty} \varphi^*(\vec{p}), \vec{p}, \varphi(\vec{p}), \mathrm{d}^3p

可以得到动量表象下的动量算符为

p^=p,p^x=px,p^y=py,p^z=pz \hat{\vec{p}} = \vec{p},\kern 12pt\hat{p}_x = p_x,\kern 12pt\hat{p}_y = p_y,\kern 12pt\hat{p}_z = p_z

通过表象的转换,可以推得坐标表象下的动量算符为

p^=i,p^x=ix,p^y=iy,p^z=iz \hat{\vec{p}} = -\mathrm{i}\hbar\nabla,\kern 12pt\hat{p}_x = -\mathrm{i}\hbar \frac{\partial}{\partial x},\kern 12pt\hat{p}_y = -\mathrm{i}\hbar \frac{\partial}{\partial y},\kern 12pt\hat{p}_z = -\mathrm{i}\hbar \frac{\partial}{\partial z}

动量表象下的坐标算符为

r^=ip,x^=ipx,y^=ipy,z^=ipz \hat{\vec{r}} = \mathrm{i}\hbar \frac{\partial}{\partial \vec{p}},\kern 12pt\hat{x} = \mathrm{i}\hbar \frac{\partial}{\partial p_x},\kern 12pt\hat{y} = \mathrm{i}\hbar \frac{\partial}{\partial p_y},\kern 12pt\hat{z} = \mathrm{i}\hbar \frac{\partial}{\partial p_z}

注:梯度、散度、旋度的介绍以及其在各种坐标系下的表示,可参考魏斌老师的课件。这里给出柱坐标与球坐标下的梯度算符,

柱坐标 (r,ϕ,z)(r,\phi,z)

f=frer+1rfθeθ+fzez \nabla f = \frac{\partial f}{\partial r}\vec{e}r + \frac{1}{r}\frac{\partial f}{\partial \theta}\vec{e}\theta + \frac{\partial f}{\partial z}\vec{e}_z

球坐标 (r,θ,φ)(r,\theta,\varphi)

f=frer+1rfθeθ+1rsinθfφeφ \nabla f = \frac{\partial f}{\partial r}\vec{e}r + \frac{1}{r}\frac{\partial f}{\partial \theta}\vec{e}\theta + \frac{1}{r\sin\theta}\frac{\partial f}{\partial \varphi}\vec{e}_\varphi

力学量算符

对于有经典对应的力学量,例如动能、势能和轨道角动量,由经典力学中的函数形式假定量子力学中的算符形式,可以由坐标算符与动量算符通过运算得到,即

A=A(r,p)A^=A(r^,p^) A = A(\vec{r},\vec{p})\Longrightarrow\hat{A} = A(\hat{\vec{r}},\hat{\vec{p}})

如一维谐振子的能量算符

H=(px)22m+12kx2H^=(p^x)22m+12kx^2 H = \frac{(p_x)^2}{2m} + \frac12kx^2\Longrightarrow\hat{H} = \frac{(\hat{p}_x)^2}{2m} + \frac12k\hat{x}^2

如粒子的轨道角动量算符

L=r×pL^=r^×p^=ijkx^y^z^p^xp^yp^z \vec{L} = \vec{r} \times \vec{p}\Longrightarrow\hat{\vec{L}} = \hat{\vec{r}} \times \hat{\vec{p}}=\begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \ \hat{x} & \hat{y} & \hat{z} \ \hat{p}_x & \hat{p}_y & \hat{p}_z\end{vmatrix}

L^x=y^p^zz^p^yL^y=z^p^xx^p^zL^z=x^p^yy^p^x \hat{L}_x = \hat{y}\hat{p}_z - \hat{z}\hat{p}_y\\kern 12pt\\hat{L}_y = \hat{z}\hat{p}_x - \hat{x}\hat{p}_z\\kern 12pt\\hat{L}_z = \hat{x}\hat{p}_y - \hat{y}\hat{p}_x

在坐标表象下,上述算符的表达式为

H^=22md2dx2+12kx2 \hat{H} = -\frac{\hbar^2}{2m} \frac{\mathrm{d}^2}{\mathrm{d}x^2} + \frac12kx^2

L^=r×(i)=ijkxyzixiyiz \hat{\vec{L}} = \vec{r} \times (-\mathrm{i}\hbar\nabla)=\begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \ x & y & z \ -\mathrm{i}\hbar \frac{\partial}{\partial x} & -\mathrm{i}\hbar \frac{\partial}{\partial y} & -\mathrm{i}\hbar \frac{\partial}{\partial z}\end{vmatrix}

L^x=i(yzzy)L^y=i(zxxz)L^z=i(xyyx) \hat{L}_x = -\mathrm{i}\hbar (y\frac{\partial}{\partial z}-z\frac{\partial}{\partial y})\\kern 12pt\\hat{L}_y = -\mathrm{i}\hbar (z\frac{\partial}{\partial x}-x\frac{\partial}{\partial z})\\kern 12pt\\hat{L}_z = -\mathrm{i}\hbar (x\frac{\partial}{\partial y}-y\frac{\partial}{\partial x})

若使用球坐标系,角动量算符表示为

L^x=i(sinφz+cotθcosφy)L^y=i(cosφθ+cotθsinφφ)L^z=iφL^2=2[1sinθθ(sinθθ)+1sin2θ2φ2] \hat{L}_x = \mathrm{i}\hbar \left( \sin\varphi\frac{\partial}{\partial z}+\cot\theta\cos\varphi\frac{\partial}{\partial y} \right) \ \kern 1em \ \hat{L}_y = \mathrm{i}\hbar \left( -\cos\varphi\frac{\partial}{\partial\theta}+\cot\theta\sin\varphi\frac{\partial}{\partial\varphi} \right) \ \kern 1em \ \hat{L}_z = -\mathrm{i}\hbar \frac{\partial}{\partial\varphi} \ \kern 1em \ \hat{L}^2 = -\hbar^2 \left[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\varphi^2} \right]

对于已归一化的波函数,力学量 AA 在坐标表象与动量表象下的平均值表达式分别为

Aˉ=+ψ(r)A(r,i)ψ(r)d3r \bar{A} = \int_{-\infty}^{+\infty} \psi^*(\vec{r}), A(\vec{r},-\mathrm{i}\hbar\nabla), \psi(\vec{r}), \mathrm{d^3}r

Aˉ=+φ(p)A(ip,p)φ(p)d3p \bar{A} = \int_{-\infty}^{+\infty} \varphi^*(\vec{p}), A(\mathrm{i}\hbar \frac{\partial}{\partial \vec{p}},\vec{p}), \varphi(\vec{p}), \mathrm{d^3}p

算符的运算规则与一些算符类型

对于本小节的理解,可以参照对线性代数的矩阵的理解。

类型:线性算符

对于任意复数 c1,c2c_1,c_2 ,任意波函数 ψ1,ψ2\psi_1,\psi_2 ,满足下列运算规则的算符 A^\hat{A} 称为线性算符

A^(c1ψ1+c2ψ2)=c1A^ψ1+c2A^ψ2 \hat{A} (c_1\psi_1+c_2\psi_2) = c_1\hat{A}\psi_1 + c_2\hat{A}\psi_2

刻画可观测量的算符都是线性算符。

并非所有算符都是线性算符,如取复共轭算符,在一般情况下 (c1ψ1+c2ψ2)=c1ψ1+c2ψ2c1ψ1+c2ψ2(c_1\psi_1+c_2\psi_2)^* = c_1^\psi_1^ + c_2^\psi_2^ \ne c_1\psi_1^* + c_2\psi_2^*

运算:算符相等

若两个算符 A^,B^\hat{A},\hat{B}任意波函数 ψ\psi 的运算所得结果都相同,即

A^ψ=B^ψ \hat{A} \psi = \hat{B} \psi

则称这两个算符相等,记为 A^=B^\hat{A}=\hat{B}

算符相等的定义给出了计算或化简算符表达式的方法,即将算符表达式作用于波函数上之后再进行计算;若不作用于波函数上直接计算,有可能会计算错误。同时也应注意只有对任意波函数都成立,才可说明算符相等,只有个别波函数成立则无法说明。

类型:单位算符

单位算符 I^\hat{I} ,是指保持任意波函数 ψ\psi 不变的运算,即

I^ψ=ψ \hat{I} \psi = \psi

运算:算符之和

算符 A^\hat{A}B^\hat{B} 之和,记为 A^+B^\hat{A}+\hat{B} ,定义如下:对于任意波函数 ψ\psi ,有

(A^+B^)ψ=A^ψ+B^ψ (\hat{A}+\hat{B}) \psi = \hat{A} \psi + \hat{B} \psi

算符之和满足交换律结合律

A^+B^=B^+A^A^+(B^+C^)=(A^+B^)+C^ \hat{A} + \hat{B} = \hat{B} + \hat{A} \ \hat{A} + (\hat{B} + \hat{C}) = (\hat{A} + \hat{B}) + \hat{C}

两个线性算符之和仍为线性算符。

运算:算符之积

算符 A^\hat{A}B^\hat{B} 之积,记为 A^B^\hat{A}\hat{B} ,定义如下:对于任意波函数 ψ\psi ,有

(A^B^)ψ=A^(B^ψ) (\hat{A}\hat{B}) \psi = \hat{A} (\hat{B}\psi)

A^B^\hat{A}\hat{B}ψ\psi 的运算结果,等于先用 B^\hat{B}ψ\psi 运算,再用 A^\hat{A}(B^ψ)(\hat{B}\psi) 运算得到的结果。

一般算符之积不满足交换律,即

A^B^B^A^ \hat{A} \hat{B} \ne \hat{B} \hat{A}

任意算符与单位算符之间可交换,即 A^I^=I^A^\hat{A}\hat{I}=\hat{I}\hat{A}

运算:算符的幂运算

A^n=A^A^A^nA^ \hat{A}^n = \underbrace{\hat{A}\hat{A}\cdots\hat{A}}_{n个\hat{A}}

A^mA^n=A^m+n \hat{A}^m \hat{A}^n = \hat{A}^{m+n}

运算:逆算符

已知算符 A^\hat{A} 与波函数 φ\varphi ,若由 A^ψ=φ\hat{A}\psi=\varphi 可以唯一地解出波函数 ψ\psi ,则可以定义算符 A^\hat{A}逆算符 A^1\hat{A}^{-1}

A^1φ=ψ \hat{A}^{-1} \varphi = \psi

并非所有的算符都有逆算符,如投影算符就不存在逆,因为 A^ψ=φ\hat{A}\psi=\varphi 的解不唯一。

A^,B^\hat{A},\hat{B} 的逆 A^1,B^1\hat{A}^{-1},\hat{B}^{-1} 存在,则

A^A^1=A^1A^=I^(A^B^)1=B^1A^1 \hat{A} \hat{A}^{-1} = \hat{A}^{-1} \hat{A} = \hat{I}\ (\hat{A}\hat{B})^{-1} = \hat{B}^{-1} \hat{A}^{-1} \

运算:转置算符

这一部分在课上没有涉及。

对于任意的波函数 ψ,φ\psi,\varphi ,算符 A^\hat{A}转置算符 A^~\widetilde{\hat{A}} 定义为

(ψ,A^~φ)=(φ,A^ψ) (\psi,\widetilde{\hat{A}}\varphi) = (\varphi^,\hat{A}\psi^)

用积分表达为

dτ ψA^~φ=dτ φA^ψ \int \mathrm{d}\tau\ \psi^* \widetilde{\hat{A}} \varphi = \int \mathrm{d}\tau\ \varphi \hat{A} \psi^*

即在积分式中,用 A^~\widetilde{\hat{A}} 作用于 φ\varphi 相当于用 A^\hat{A} 作用于 ψ\psi^*

可以证明,

A^B^=B^~ \widetilde{\hat{A}\hat{B}} = \widetilde{\hat{B}} \widetilde{}

运算:复共轭算符

算符 A^\hat{A}复共轭算符 A^\hat{A}^* 定义为

A^ψ=(A^ψ) \hat{A}^* \psi = (\hat{A}\psi^)^

通常算符 A^\hat{A} 的复共轭可通过把 A^\hat{A} 的表达式中所有量换成其复共轭得到,且算符 A^\hat{A}^* 的表达式与表象有关,如在坐标表象中

p^=(i)=i=p^ \hat{\vec{p}}^* = (-\mathrm{i} \hbar \nabla)^* = \mathrm{i} \hbar \nabla = -\hat{\vec{p}}

在动量表象中

p^=p^ \hat{\vec{p}}^* = \hat{\vec{p}}

运算:厄米共轭算符

算符 A^\hat{A}厄米共轭算符 A^+\hat{A}^+ 定义为

(ψ,A^+φ)=(A^ψ,φ) (\psi , \hat{A}^+ \varphi) = (\hat{A} \psi , \varphi)

用积分表达为

dτ ψ(A^+φ)=dτ (A^ψ)φ \int \mathrm{d}\tau\ \psi^* (\hat{A}^+ \varphi) = \int \mathrm{d}\tau\ (\hat{A} \psi)^* \varphi

厄米共轭算符有如下性质:

(A^+)+=A^(A^+B^)+=A^++B^+(A^B^)+=B^+A^+ (\hat{A}^+)^+ = \hat{A} \ (\hat{A} + \hat{B})^+ = \hat{A}^+ + \hat{B}^+ \ (\hat{A}\hat{B})^+ = \hat{B}^+ \hat{A}^+

类型:厄米算符

对于任意波函数 ψ,φ\psi,\varphi ,若算符 A^\hat{A} 满足

(ψ,A^φ)=(A^ψ,φ) (\psi , \hat{A} \varphi) = (\hat{A}\psi , \varphi)

A^\hat{A}厄米算符,也称自共轭算符

上述条件用积分表达为

dτ ψ(A^φ)=dτ (A^ψ)φ \int \mathrm{d}\tau\ \psi^* (\hat{A} \varphi) = \int \mathrm{d}\tau\ (\hat{A} \psi)^* \varphi

与厄米共轭算符的定义相对比,可知若 A^+=A^\hat{A}^+ = \hat{A} ,则 A^\hat{A}厄米算符

两个厄米算符 A^,B^\hat{A},\hat{B} 之和仍为厄米算符,即 (A^+B^)+=A^++B^+=A^+B^(\hat{A}+\hat{B})^+ = \hat{A}^+ + \hat{B}^+ = \hat{A} + \hat{B} ;当且仅当两个厄米算符 A^,B^\hat{A},\hat{B} 可交换时,其积为厄米算符,这是因为 (A^B^)+=B^+A^+=B^A^(\hat{A}\hat{B})^+ = \hat{B}^+ \hat{A}^+ = \hat{B}\hat{A} ,当且仅当 A^,B^\hat{A},\hat{B} 可交换时, (A^B^)+=A^B^(\hat{A}\hat{B})^+ = \hat{A}\hat{B}

类型:幺正算符

若算符 A^\hat{A} 满足

A^+=A^1 \hat{A}^{+} = \hat{A}^{-1}

则称 A^\hat{A} 为幺正算符。

例:空间反演算符 P^\hat{P} 表达体系的宇称,是厄米算符, P^+=P^\hat{P}^{+} = \hat{P} ;根据 P^P^=I^P^=P^1P^+=P^1\hat{P}\hat{P} = \hat{I} \Longrightarrow \hat{P} = \hat{P}^{-1} \Longrightarrow \hat{P}^+ = \hat{P}^{-1} 可知 P^\hat{P} 也是幺正算符。

运算:算符的函数

给定一函数 F(x)F(x) ,其各阶导数均存在,幂级数展开收敛,

F(x)=n=0+F(n)(0)n!xn F(x) = \sum_{n=0}^{+\infty} \frac{F^{(n)}(0)}{n!} x^n

可以定义算符 A^\hat{A} 的函数 F(A^)F(\hat{A})

F(A^)=n=0+F(n)(0)n!A^n F(\hat{A}) = \sum_{n=0}^{+\infty} \frac{F^{(n)}(0)}{n!} \hat{A}^n

两个(或多个)算符的函数也可以类似定义,如

F(A^,B^)=m,n=0+1m!n!m+nFxmyn(0,0) A^mB^n F(\hat{A},\hat{B}) = \sum_{m,n=0}^{+\infty} \frac{1}{m!n!} \frac{\partial^{m+n}F}{\partial x^m\partial y^n} (0,0)\ \hat{A}^m\hat{B}^n

算符的对易式

对易式的定义

定义算符 A^,B^\hat{A},\hat{B}对易式

[A^,B^]=A^B^B^A^ [\hat{A},\hat{B}] = \hat{A}\hat{B} - \hat{B}\hat{A}

[A^,B^]=0[\hat{A},\hat{B}] = 0 ,即 A^B^=B^A^\hat{A}\hat{B} = \hat{B}\hat{A} ,则称 A^\hat{A}B^\hat{B} 对易(可交换)。显然一个算符与它本身对易,即 [A^,A^]=0[\hat{A},\hat{A}] = 0

若要计算一个对易式 [A^,B^][\hat{A},\hat{B}] ,可以使用作用法,即 [A^,B^]ψ=C^ψ[A^,B^]=C^[\hat{A},\hat{B}]\psi = \hat{C}\psi \Longrightarrow [\hat{A},\hat{B}] = \hat{C} ;也可以由对易式的运算规则直接计算。

对易式的运算规则

 [A^,B^]=[B^,A^] [A^,B^+C^]=[A^,B^]+[A^,C^] [A^,B^C^]=B^[A^,C^]+[A^,B^]C^ [A^B^,C^]=A^[B^,C^]+[A^,C^]B^ [A^,[B^,C^]]+[B^,[C^,A^]]+[C^,[A^,B^]]=0 \ [\hat{A},\hat{B}] = -[\hat{B},\hat{A}] \ \ [\hat{A},\hat{B}+\hat{C}] = [\hat{A},\hat{B}] + [\hat{A},\hat{C}] \ \ [\hat{A},\hat{B}\hat{C}] = \hat{B}[\hat{A},\hat{C}] + [\hat{A},\hat{B}]\hat{C} \ \ [\hat{A}\hat{B},\hat{C}] = \hat{A}[\hat{B},\hat{C}] + [\hat{A},\hat{C}]\hat{B} \ \ [\hat{A},[\hat{B},\hat{C}]] + [\hat{B},[\hat{C},\hat{A}]] + [\hat{C},[\hat{A},\hat{B}]] = 0

其中最后一个式子称为Jaccobi恒等式。

最基本的对易关系:坐标与动量的对易关系

坐标的各个分量算符之间对易,动量的各个分量算符之间对易。

 [x^,y^]=[x^,z^]=[y^,z^]=0 [p^x,p^y]=[p^x,p^z]=[p^y,p^z]=0 \ [\hat{x},\hat{y}] = [\hat{x},\hat{z}] = [\hat{y},\hat{z}] = 0 \ \kern 1em \ \ [\hat{p}_x,\hat{p}_y] = [\hat{p}_x,\hat{p}_z] = [\hat{p}_y,\hat{p}_z] = 0

坐标算符与动量算符之间的对易关系为

[x^α,p^β]=iδαβ(α,β=x,y,z1,2,3) [\hat{x}{\alpha} , \hat{p}{\beta}] = \mathrm{i} \hbar \delta_{\alpha\beta} \kern 2em (\alpha,\beta = x,y,z 或 1,2,3)

分开来看,同分量的坐标算符与动量算符不对易

[x^,p^x]=[y^,p^y]=[z^,p^z]=i [\hat{x} , \hat{p}{x}] = [\hat{y} , \hat{p}{y}] = [\hat{z} , \hat{p}_{z}] = \mathrm{i} \hbar

不同分量的坐标算符与动量算符对易

[x^,p^y]=[x^,p^z]=[y^,p^x]=[y^,p^z]=[z^,p^x]=[z^,p^y]=0 [\hat{x} , \hat{p}{y}] = [\hat{x} , \hat{p}{z}] = [\hat{y} , \hat{p}{x}] = [\hat{y} , \hat{p}{z}] = [\hat{z} , \hat{p}{x}] = [\hat{z} , \hat{p}{y}] = 0

角动量的对易关系

角动量算符三个分量之间的对易关系为

[L^α,L^β]=εαβγiL^γ [\hat{L}{\alpha} , \hat{L}{\beta}] = \varepsilon_{\alpha\beta\gamma} \mathrm{i} \hbar \hat{L}_{\gamma}

式中 εαβγ\varepsilon_{\alpha\beta\gamma} 称为 Levi-Civita 符号,是一个三阶反对称张量,定义如下

{εαβγ=εβαγ=εαγβε123=1 \begin{cases} \varepsilon_{\alpha\beta\gamma} = - \varepsilon_{\beta\alpha\gamma} = - \varepsilon_{\alpha\gamma\beta} \ \varepsilon_{123} = 1 \end{cases}

其中 α,β,γ=1,2,3x,y,z\alpha,\beta,\gamma=1,2,3或x,y,z ,上式中第一个式子的含义是任何两个指标交换时 εαβγ\varepsilon_{\alpha\beta\gamma} 改变正负号,由此可得任何两个指标相同时 εαβγ=0\varepsilon_{\alpha\beta\gamma}=0

上述的三阶反对称张量可认为表示如下

[[000001010],[001000100], [010100000]] \begin{bmatrix} \begin{bmatrix} 0 & 0 & 0 \ 0 & 0 & 1 \ 0 & -1 & 0 \end{bmatrix} , \begin{bmatrix} 0 & 0 & -1 \ 0 & 0 & 0 \ 1 & 0 & 0 \end{bmatrix} ,\ \begin{bmatrix} 0 & 1 & 0 \ -1 & 0 & 0 \ 0 & 0 & 0 \end{bmatrix} \end{bmatrix}

分开来看,不同分量的角动量算符之间不对易

[L^x,L^y]=iL^z,[L^y,L^z]=iL^x,[L^z,L^x]=iL^y [\hat{L}{x} , \hat{L}{y}] = \mathrm{i} \hbar \hat{L}{z} ,\kern 1em [\hat{L}{y} , \hat{L}{z}] = \mathrm{i} \hbar \hat{L}{x} ,\kern 1em [\hat{L}{z} , \hat{L}{x}] = \mathrm{i} \hbar \hat{L}_{y}

可以根据该对易关系拓展角动量算符的定义:若一个矢量算符的三个分量满足上述对易关系,则这个算符就是角动量算符。

角动量算符的三个分量都和角动量的平方对易

[L^x,L^2]=[L^y,L^2]=[L^z,L^2]=0 [\hat{L}{x} , \hat{L}^2] = [\hat{L}{y} , \hat{L}^2] = [\hat{L}_{z} , \hat{L}^2] = 0

角动量算符与坐标、动量算符之间满足类似于角动量的三个分量之间的对易关系

 [L^α,x^β]=εαβγix^γ [L^α,p^β]=εαβγip^γ \ [\hat{L}{\alpha} , \hat{x}{\beta}] = \varepsilon_{\alpha\beta\gamma} \mathrm{i} \hbar \hat{x}{\gamma} \ \kern 1em \ \ [\hat{L}{\alpha} , \hat{p}{\beta}] = \varepsilon{\alpha\beta\gamma} \mathrm{i} \hbar \hat{p}_{\gamma}

分开来看,同分量的角动量算符与坐标、动量算符对易不同分量的角动量算符与坐标、动量算符不对易

[l^x,x^]=0[l^x,y^]=iz[l^x,z^]=iy[l^y,x^]=iz[l^y,y^]=0[l^y,z^]=ix[l^z,x^]=iy[l^z,y^]=ix[l^z,z^]=0[l^x,p^x]=0[l^x,p^y]=ipz[l^x,p^z]=ipy[l^y,p^x]=ipz[l^y,p^y]=0[l^y,p^z]=ipx[l^z,p^x]=ipy[l^z,p^y]=ipx[l^z,p^z]=0 \begin{matrix} [\hat{l}_x,\hat{x}] = 0 & [\hat{l}_x,\hat{y}] = \mathrm{i}\hbar z & [\hat{l}_x,\hat{z}] = -\mathrm{i}\hbar y \ [\hat{l}_y,\hat{x}] = -\mathrm{i}\hbar z & [\hat{l}_y,\hat{y}] = 0 & [\hat{l}_y,\hat{z}] = \mathrm{i}\hbar x \ [\hat{l}_z,\hat{x}] = \mathrm{i}\hbar y & [\hat{l}_z,\hat{y}] = -\mathrm{i}\hbar x & [\hat{l}_z,\hat{z}] = 0 \end{matrix} \ \kern 1em \ \begin{matrix} [\hat{l}_x,\hat{p}_x] = 0 & [\hat{l}_x,\hat{p}_y] = \mathrm{i}\hbar p_z & [\hat{l}_x,\hat{p}_z] = -\mathrm{i}\hbar p_y \ [\hat{l}_y,\hat{p}_x] = -\mathrm{i}\hbar p_z & [\hat{l}_y,\hat{p}_y] = 0 & [\hat{l}_y,\hat{p}_z] = \mathrm{i}\hbar p_x \ [\hat{l}_z,\hat{p}_x] = \mathrm{i}\hbar p_y & [\hat{l}_z,\hat{p}_y] = -\mathrm{i}\hbar p_x & [\hat{l}_z,\hat{p}_z] = 0 \end{matrix}

3.2 厄米算符的本征值与本征函数

算符的本征方程

算符的本征方程及其意义

对于算符 A^\hat{A} ,有如下本征方程

A^ψλ=λψλ \hat{A} \psi_{\lambda} = \lambda \psi_{\lambda}

算符 A^\hat{A}本征值{λ}{\lambda} 就是力学量 AA 的测量值集,算符 A^\hat{A}本征函数 ψλ\psi_{\lambda} 代表力学量 AA 在本征值 λ\lambda 下的状态。

本征值的类型

本征值可以是分立谱(discrete spectra)、连续谱(continuous spectra)和混合谱

分立谱的本征值 AnA_n 是离散的,本征方程表示为

A^ψn=Anψn ,(n=1,2,3,) \hat{A} \psi_n = A_n \psi_n\ , \kern 2em (n=1,2,3,\cdots)

连续谱的本征值 λ\lambda 是连续的,本征方程表示为

A^ψλ=λψλ \hat{A} \psi_{\lambda} = \lambda \psi_{\lambda}

混合谱是由一部分分立谱与一部分离散谱组成的,如氢原子的能级,在电离前(束缚态)是离散谱,电离后(游离态)是连续谱。

厄米算符的本征值与平均值

厄米算符的本征值

厄米算符的本征值必为实数

证明

对于厄米算符 A^\hat{A} 的本征值 λ\lambda 与本征函数 ψ\psi ,有

A^ψ=λψ \hat{A} \psi = \lambda \psi

由厄米算符定义可得

dτ ψ(A^ψ)=dτ (A^ψ)ψdτ ψ(λψ)=dτ (λψ)ψλdτ ψψ=λdτ ψψλ=λ \int \mathrm{d}\tau\ \psi^* (\hat{A} \psi) = \int \mathrm{d}\tau\ (\hat{A} \psi)^* \psi \ \Downarrow \ \int \mathrm{d}\tau\ \psi^* (\lambda \psi) = \int \mathrm{d}\tau\ (\lambda \psi)^* \psi \ \Downarrow \ \lambda \int \mathrm{d}\tau\ \psi^* \psi = \lambda^* \int \mathrm{d}\tau\ \psi^* \psi \ \Downarrow \ \lambda = \lambda^*

此即表明本征值 λ\lambda 为实数。

厄米算符的平均值

体系的任何状态下,其厄米算符的平均值必为实数;在任何状态下平均值均为实数的算符必为厄米算符。即

厄米算符任何状态下平均值为实数 厄米算符 \Longleftrightarrow 任何状态下平均值为实数

由该定理可知,力学量即可观测量,当然要求在任何状态下平均值都是实数,所以相应的算符一定是厄米算符

证明

\Rightarrow 证明如下:

对于厄米算符 A^\hat{A} 与任意波函数 ψ\psi ,有

Aˉ=(ψ,A^ψ)=(A^ψ,ψ)=(ψ,A^ψ)=Aˉ \bar{A} = (\psi,\hat{A}\psi) = (\hat{A}\psi,\psi) = (\psi,\hat{A}\psi)^* = \bar{A}^*

其中第二个等号为厄米算符的性质,第三个等号为内积的性质。

\Leftarrow 证明如下:

对于任意波函数 ψ\psi ,有 Aˉ=Aˉ\bar{A} = \bar{A}^* ,由 Aˉ=(ψ,A^ψ)\bar{A} = (\psi,\hat{A}\psi)Aˉ=(ψ,A^ψ)=(A^ψ,ψ)\bar{A}^* = (\psi,\hat{A}\psi)^* = (\hat{A}\psi,\psi) 可得

(ψ,A^ψ)=(A^ψ,ψ) (\psi,\hat{A}\psi) = (\hat{A}\psi,\psi)

注意此时并未完成证明,因为厄米算符的定义是对两个独立的波函数 ψ,φ\psi,\varphi 而言的。

取任意独立的波函数 ψ1,ψ2\psi_1,\psi_2 与任意的复数 cc ,令 ψ=ψ1+cψ2\psi = \psi_1 + c\psi_2 ,则

(ψ,A^ψ)=(ψ1+cψ2,A^(ψ1+cψ2))=(ψ1,A^ψ1)+c(ψ1,A^ψ2)+c(ψ2,A^ψ1)+cc(ψ2,A^ψ2)(A^ψ,ψ)=(A^(ψ1+cψ2),ψ1+cψ2)=(A^ψ1,ψ1)+c(A^ψ1,ψ2)+c(A^ψ2,ψ1)+cc(A^ψ2,ψ2) (\psi , \hat{A}\psi) = \left(\psi_1 + c\psi_2 , \hat{A}(\psi_1 + c\psi_2)\right) \ = (\psi_1 , \hat{A}\psi_1) + c(\psi_1 , \hat{A}\psi_2) + c^*(\psi_2 , \hat{A}\psi_1) + c^c(\psi_2 , \hat{A}\psi_2) \ \kern 1em \ (\hat{A}\psi,\psi) = \left(\hat{A}(\psi_1 + c\psi_2) , \psi_1 + c\psi_2\right) \ = (\hat{A}\psi_1 , \psi_1) + c(\hat{A}\psi_1 , \psi_2) + c^(\hat{A}\psi_2 , \psi_1) + c^*c(\hat{A}\psi_2 , \psi_2)

两组式子应该相等,并利用 (ψ1,A^ψ1)=(A^ψ1,ψ1), (ψ2,A^ψ2)=(A^ψ2,ψ2)(\psi_1,\hat{A}\psi_1) = (\hat{A}\psi_1,\psi_1),\ (\psi_2,\hat{A}\psi_2) = (\hat{A}\psi_2,\psi_2) ,可得

c(ψ1,A^ψ2)+c(ψ2,A^ψ1)=c(A^ψ1,ψ2)+c(A^ψ2,ψ1) c(\psi_1 , \hat{A}\psi_2) + c^(\psi_2 , \hat{A}\psi_1) = c(\hat{A}\psi_1 , \psi_2) + c^(\hat{A}\psi_2 , \psi_1)

整理得

c[(ψ1,A^ψ2)(A^ψ1,ψ2)]=c[(A^ψ2,ψ1)(ψ2,A^ψ1)] c\left[(\psi_1 , \hat{A}\psi_2) - (\hat{A}\psi_1 , \psi_2)\right] = c^*\left[(\hat{A}\psi_2 , \psi_1) - (\psi_2 , \hat{A}\psi_1)\right]

分别令 c=1c=1c=ic=i ,可得

(ψ1,A^ψ2)(A^ψ1,ψ2)=+(A^ψ2,ψ1)(ψ2,A^ψ1)(ψ1,A^ψ2)(A^ψ1,ψ2)=(A^ψ2,ψ1)+(ψ2,A^ψ1) (\psi_1 , \hat{A}\psi_2) - (\hat{A}\psi_1 , \psi_2) = + (\hat{A}\psi_2 , \psi_1) - (\psi_2 , \hat{A}\psi_1) \ \kern 1em \ (\psi_1 , \hat{A}\psi_2) - (\hat{A}\psi_1 , \psi_2) = - (\hat{A}\psi_2 , \psi_1) + (\psi_2 , \hat{A}\psi_1)

以上两式分别相加、减,即得

(ψ1,A^ψ2)=(A^ψ1,ψ2),(A^ψ2,ψ1)=(ψ2,A^ψ1) (\psi_1 , \hat{A}\psi_2) = (\hat{A}\psi_1 , \psi_2), \kern 1em (\hat{A}\psi_2 , \psi_1) = (\psi_2 , \hat{A}\psi_1)

此即厄米算符定义的要求。

推论

A^\hat{A} 是厄米算符,则在任意态 ψ\psi 之下,

A2ˉ=(ψ,A^2ψ)=(A^ψ,A^ψ)0 \bar{A^2} = (\psi,\hat{A}^2\psi) = (\hat{A}\psi,\hat{A}\psi) \ge 0

厄米算符本征函数系的正交性

正交性定理

厄米算符的属于不同本征值的本征函数,彼此正交

对于分立谱

(ψm,ψn)=δmn (\psi_m,\psi_n) = \delta_{mn}

对于连续谱

(ψλ,ψλ)=δ(λλ) (\psi_{\lambda},\psi_{\lambda'}) = \delta(\lambda-\lambda')

连续谱本征函数是不能归一化的,可以使用 δ\delta 函数进行规格化(如上),或者使用箱归一化。

证明

以分立谱为例证明,连续谱同理。对于厄米算符 A^\hat{A} 与以归一化的本征函数 ψm,ψn\psi_m,\psi_n ,设

A^ψn=Anψn ,A^ψm=Amψm \hat{A} \psi_n = A_n \psi_n\ , \kern 2em \hat{A} \psi_m = A_m \psi_m

并设 (ψm,ψn)(\psi_m,\psi_n) 存在,对第二个式子取复共轭,得

A^ψm=Amψm \hat{A}^* \psi_m^* = A_m \psi_m^*

ψn\psi_n ,并积分,得

dτA^ψmψn=dτAmψmψn \int \mathrm{d} \tau \hat{A}^* \psi_m^* \psi_n = \int \mathrm{d} \tau A_m \psi_m^* \psi_n

(A^ψm,ψn)=Am(ψm,ψn) (\hat{A}\psi_m,\psi_n) = A_m(\psi_m,\psi_n)

A^\hat{A} 是厄米算符可知 (A^ψm,ψn)=(ψm,A^ψn)=An(ψm,ψn)(\hat{A}\psi_m,\psi_n) = (\psi_m,\hat{A}\psi_n) = A_n(\psi_m,\psi_n) ,故

(AmAn)(ψm,ψn)=0 (A_m - A_n) (\psi_m,\psi_n) = 0

mnm \ne nAmAnA_m \ne A_n 时,则有

(ψm,ψn)=0 (\psi_m,\psi_n) = 0

简并态之间的正交性

在出现简并态时,简并态的选择不是唯一的,而且选择的简并态不一定彼此正交。但可以证明,总可以把它们适当地线性叠加,使之彼此正交

这里给出简并态的选择不是唯一的一个例子,对于本征方程 d2dx2ψ(x)+k2ψ(x)=0\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi(x) + k^2\psi(x) =0 ,既可以选择 e±ikxe^{\pm\mathrm{i}kx} ,也可以选择 sinx\sin xcosx\cos x ,甚至可以选择 sinx+cosx\sin x+\cos xsinxcosx\sin x-\cos x 等。

与线性代数作类比,即选择的子空间的一组基不一定正交,但总可以找到子空间的一组标准正交基(这组正交基也可能不唯一)。

通过正交性定理与该定理,结合后面对完备性的讨论,实际上说明了厄米算符可以选择正交、归一、完备的本征函数系。

证明

设力学量 AA 的本征方程表示为

A^ψnα=Anψnα ,(a=1,2,,fn) \hat{A} \psi_{n\alpha} = A_n \psi_{n\alpha}\ , \kern 2em (a=1,2,\cdots,f_n)

即属于本征值 AnA_n 的本征态有 fnf_n 个(称本征值 AnA_nfnf_n 重简并)。令

ϕnβ=α=1fncβαψnα ,(β=1,2,,fn) \phi_{n\beta} = \sum_{\alpha=1}^{f_n} c_{\beta\alpha} \psi_{n\alpha}\ , \kern 2em (\beta=1,2,\cdots,f_n)

这里得到的 ϕnβ\phi_{n\beta} 仍为 A^\hat{A} 的本征态(不是叠加态),相应的本征值仍为 AnA_n ,因为

A^ϕnβ=A^α=1fncβαψnα=α=1fncβαA^ψnα=α=1fncβαAnψnα=Anα=1fncβαψnα=Anϕnβ \hat{A} \phi_{n\beta} = \hat{A} \sum_{\alpha=1}^{f_n} c_{\beta\alpha} \psi_{n\alpha} = \sum_{\alpha=1}^{f_n} c_{\beta\alpha} \hat{A} \psi_{n\alpha} = \sum_{\alpha=1}^{f_n} c_{\beta\alpha} A_n \psi_{n\alpha} = A_n \sum_{\alpha=1}^{f_n} c_{\beta\alpha} \psi_{n\alpha} = A_n \phi_{n\beta}

可以适当地选择 cβαc_{\beta\alpha} ,使得 ϕnβ\phi_{n\beta} 具有正交性,即

(ϕnβ,ϕnβ)=δββ (\phi_{n\beta} , \phi_{n\beta'}) = \delta_{\beta\beta'}

对于 fn2f_n^2 个系数 cβαc_{\beta\alpha} ,这相当于提出了 12fn(fn1)+fn=12fn(fn+1)\frac12f_n(f_n-1)+f_n = \frac12f_n(f_n+1) 个线性方程(不同的 β,β\beta,\beta'Cfn2=12fn(fn1)C_{f_n}^{2} = \frac12f_n(f_n-1) 个,相同的 β,β\beta,\beta'fnf_n 个),而 fn212fn(fn+1)f_n^2 \ge \frac12f_n(f_n+1) ,则这个线性方程组应该是有解的(并且当 fn2f_n\ge2 时解应该是不唯一的)。具体的系数求解过程可以通过Schmidt正交化实现。

厄米算符本征函数系的完备性

函数系的完备性
完备性的涵义

一个函数系完备,是指任何一个满足适当边界条件和连续性要求的波函数,均可用这个函数系作展开。

从线性代数的角度来看,就是一组向量可以张成一个子空间,这个子空间中的所有向量都可以用这组向量线性表出。

只有那些本征波函数构成完备系的厄米算符所表达的力学量才是可以观测的,才有物理意义。物理上力学量总是可观测的,所以量子力学有理由认为表达力学量的厄米算符的本征函数系是完备的

注:对于常见的势函数体系,其Hamilton量的本征函数系的完备性数学上已经证明,但对任意势函数的情况目前还不能一般地证明。

在完备函数系下的展开

对于分立谱

算符 A^\hat{A} 与本征函数系 {ψn}{\psi_n} 满足

A^ψn=Anψn ,(n=1,2,3,) \hat{A} \psi_n = A_n \psi_n\ , \kern 2em (n=1,2,3,\cdots)

且函数系 {ψn}{\psi_n} 是正交、归一、完备的,(ψm,ψn)=δmn(\psi_m,\psi_n) = \delta_{mn} ,则对任意波函数 ϕ\phi 可以作展开

ϕ(x)=nCnψn(x) \phi(x) = \sum_{n} C_n\psi_n (x)

两边同乘 ψm\psi_m^* ,积分可得

dτψmϕ=dτψmnCnψn=nCndτψmψn=Cm \int \mathrm{d}\tau \psi_m^* \phi = \int \mathrm{d}\tau \psi_m^* \sum_{n} C_n\psi_n = \sum_{n}C_n \int \mathrm{d}\tau \psi_m^* \psi_n = C_m

由此可以求出展开式的系数

Cn=(ψn,ϕ)=dτψnϕ C_n = (\psi_n,\phi) = \int \mathrm{d}\tau \psi_n^* \phi

展开系数 CnC_n 是态矢 ϕ\phi 在本征矢量 ψn\psi_n 上的投影,展开系数的集合

Φ=[C1C2] \Phi = \begin{bmatrix} C_1 \ C_2 \ \vdots \end{bmatrix}

代表态矢 ϕ\phi 在基底 {ψn}{ \psi_n } 上的表示,或称在表象 A^\hat{A} 上的表示。

对于连续谱

算符 A^\hat{A} 与本征函数系 {ψλ}{\psi_{\lambda}} 满足

A^ψλ=λψλ ,(n=1,2,3,) \hat{A} \psi_\lambda = \lambda \psi_\lambda\ , \kern 2em (n=1,2,3,\cdots)

且函数系 {ψλ}{\psi_\lambda} 是正交、归一、完备的,(ψλ,ψλ)=δ(λλ)(\psi_{\lambda},\psi_{\lambda'}) = \delta(\lambda-\lambda') ,则对任意波函数 ϕ\phi 可以作展开

φ(x)=C(λ)ψλ(x)dλ \varphi(x) = \int C(\lambda) \psi_\lambda (x) \mathrm{d}\lambda

两边同乘 ψλ\psi_{\lambda'}^* ,积分可得

dτψλϕ=dτψλC(λ)ψλdλ=dλC(λ)dτψλψλ=dλC(λ)δ(λλ)=C(λ) \int \mathrm{d}\tau \psi_{\lambda'}^* \phi = \int \mathrm{d}\tau \psi_{\lambda'}^* \int C(\lambda) \psi_\lambda \mathrm{d}\lambda = \int \mathrm{d}\lambda C(\lambda) \int \mathrm{d}\tau \psi_{\lambda'}^* \psi_{\lambda} = \int \mathrm{d}\lambda C(\lambda) \delta(\lambda-\lambda') = C(\lambda')

由此可以求出展开式的系数

C(λ)=(ψλ,ϕ)=dτψλϕ C(\lambda) = (\psi_\lambda,\phi) = \int \mathrm{d}\tau \psi_\lambda^* \phi

完备性用封闭关系表示

对于分立谱,本征函数系 {ψn}{ \psi_n } 完备等价于其满足封闭关系

nψn(x) ψn(x)=δ(xx) \sum_{n} \psi_n^*(x')\ \psi_n(x) = \delta(x'-x)

对于连续谱,本征函数系 {ψλ}{ \psi_\lambda } 完备等价于其满足封闭关系

ψλ(x) ψλ(x) dλ=δ(xx) \int \psi_\lambda^*(x')\ \psi_\lambda(x)\ \mathrm{d}\lambda = \delta(x'-x)

证明

对于分立谱

若本征函数系 {ψn}{ \psi_n } 完备,则对任意的波函数 ϕ(x)\phi(x)

ϕ(x)=nCnψn(x)=n[+ψn(x)ϕ(x)dx]ψn(x)=+ϕ(x)[nψn(x)ψn(x)]dx \phi (x) = \sum_n C_n \psi_n(x) = \sum_n \left[ \int_{-\infty}^{+\infty} \psi_n^(x') \phi(x') \mathrm{d}x' \right] \psi_n(x) = \int_{-\infty}^{+\infty} \phi(x') \left[ \sum_n \psi_n^(x') \psi_n(x) \right] \mathrm{d}x'

通过上述积分式可以得到

nψn(x)ψn(x)=δ(xx) \sum_n \psi_n^*(x') \psi_n(x) = \delta(x'-x)

而若本征函数系 {ψn}{ \psi_n } 满足该封闭关系,则

ϕ(x)=+ϕ(x)δ(xx)dx=+ϕ(x)[nψn(x)ψn(x)]dx=n[+ψn(x)ϕ(x)dx]ψn(x) \phi (x) = \int_{-\infty}^{+\infty} \phi(x') \delta(x'-x) \mathrm{d}x' = \int_{-\infty}^{+\infty} \phi(x') \left[ \sum_n \psi_n^(x') \psi_n(x) \right] \mathrm{d}x' = \sum_n \left[ \int_{-\infty}^{+\infty} \psi_n^(x') \phi(x') \mathrm{d}x' \right] \psi_n(x)

这就表现为任意波函数 ϕ(x)\phi(x) 可以对本征函数系 {ψn}{ \psi_n } 作展开,即本征函数系 {ψn}{ \psi_n } 完备。

对于连续谱

若本征函数系 {ψλ}{ \psi_\lambda } 完备,则对任意的波函数 ϕ(x)\phi(x)

ϕ(x)=C(λ)ψλ(x)dλ=[+ψλ(x)ϕ(x)dx]ψλ(x)dλ=+ϕ(x)[ψλ(x)ψλ(x)dλ]dx \phi (x) = \int C(\lambda) \psi_\lambda(x) \mathrm{d}\lambda = \int \left[ \int_{-\infty}^{+\infty} \psi_\lambda^(x') \phi(x') \mathrm{d}x' \right] \psi_\lambda(x) \mathrm{d}\lambda = \int_{-\infty}^{+\infty} \phi(x') \left[ \int \psi_\lambda^(x') \psi_\lambda(x) \mathrm{d}\lambda \right] \mathrm{d}x'

通过上述积分式可以得到

ψλ(x)ψλ(x)dλ=δ(xx) \int \psi_\lambda^*(x') \psi_\lambda(x) \mathrm{d}\lambda = \delta(x'-x)

而若本征函数系 {ψλ}{ \psi_\lambda } 满足该封闭关系,则

ϕ(x)=+ϕ(x)δ(xx)dx=+ϕ(x)[ψλ(x)ψλ(x)dλ]dx=[+ψλ(x)ϕ(x)dx]ψλ(x)dλ \phi (x) = \int_{-\infty}^{+\infty} \phi(x') \delta(x'-x) \mathrm{d}x' = \int_{-\infty}^{+\infty} \phi(x') \left[ \int \psi_\lambda^(x') \psi_\lambda(x) \mathrm{d}\lambda \right] \mathrm{d}x' = \int \left[ \int_{-\infty}^{+\infty} \psi_\lambda^(x') \phi(x') \mathrm{d}x' \right] \psi_\lambda(x) \mathrm{d}\lambda

这就表现为任意波函数 ϕ(x)\phi(x) 可以对本征函数系 {ψλ}{ \psi_\lambda } 作展开,即本征函数系 {ψλ}{ \psi_\lambda } 完备。

力学量用厄米算符表达

这里是课本上关于力学量与算符关系的总结部分,同时补充了课件上关于力学量的平均值部分。

量子力学中力学量用相应的线性的厄米算符来表达,其有以下多个含义:

平均值

在给定状态 ψ\psi 之下,力学量 AA 的平均值 Aˉ\bar{A} 由下式确定:

Aˉ=(ψ,A^ψ)(ψ,ψ) \bar{A} = \frac{(\psi,\hat{A}\psi)}{(\psi,\psi)}

对于叠加态的另一种计算方式

这里直接给出混合谱的计算方式,分立谱与连续谱的计算方式均可从中取出一部分得到。设

ψ(x)=nCnψn(x)+C(λ)ψλdλ \psi(x) = \sum_n C_n\psi_n(x) + \int C(\lambda) \psi_\lambda \mathrm{d}\lambda

其中

Cn=(ψn,ψ) ,C(λ)=(ψλ,ψ) nCn2+C(λ)2dλ=1 C_n = (\psi_n,\psi) \ , \kern 1em C(\lambda) = (\psi_\lambda,\psi) \ \ \ \sum_n \left|C_n\right|^2 + \int \left|C(\lambda)\right|^2 \mathrm{d}\lambda = 1

A^ψn=Anψn ,A^ψλ=λψλ \hat{A} \psi_n = A_n \psi_n \ , \kern 1em \hat{A} \psi_\lambda = \lambda \psi_\lambda

则在状态 ψ\psi 之下,力学量 AA 的平均值

Aˉ=nAnCn2+λC(λ)2dλ \bar{A} = \sum_{n} A_n |C_n|^2 + \int \lambda |C(\lambda)|^2 \mathrm{d}\lambda

本征值

在实验上观测某力学量 AA ,它的可能取值就是算符 A^\hat{A} 的某一个本征值。由于力学量观测值总是实数,所以要求相应的算符为厄米算符。

对易关系

力学量之间的关系也通过相应的算符之间的关系反映出来。例如,两个力学量 AABB ,在一般情况下,可以同时具有确定的观测值的必要条件为 [A^,B^]=0[\hat{A},\hat{B}] = 0 ;反之,若 [A^,B^]0[\hat{A},\hat{B}] \ne 0 ,则一般说来,力学量 AABB 不能同时具有确定的观测值。

3.3 共同本征函数

共同本征态

共同本征函数的定义

若波函数 ψ\psi 同时是至少两个算符 A^,B^,\hat{A},\hat{B},\cdots 的本征函数,即

A^ψ=λψ ,B^ψ=μψ , \hat{A} \psi = \lambda \psi \ , \kern 1em \hat{B}\psi = \mu \psi \ , \kern 1em \cdots

则称 ψ\psi 为算符 A^,B^,\hat{A},\hat{B},\cdots共同本征函数(同时本征函数),也称为共同本征态。

共同本征函数与对易关系

如果算符 A^\hat{A}B^\hat{B} 有一组共同本征函数 {ψn}{ \psi_n } ,而且 {ψn}{ \psi_n } 组成完备系,则算符 A^\hat{A}B^\hat{B} 对易。(连续谱同)

如果算符 A^\hat{A}B^\hat{B} 对易,则 A^\hat{A}B^\hat{B}共同本征函数,即存在 ψ\psi 使得 A^ψ=λψ\hat{A}\psi = \lambda\psiB^ψ=μψ\hat{B}\psi = \mu\psi 同时成立。

要注意第二句话中是存在共同本征函数,也就是不意味着两个算符具有相同的共同本征函数系(二者的简并度甚至都一定不相同),也就是不意味着算符 A^\hat{A} 的本征函数一定是算符 B^\hat{B} 的本征函数;只有对于两者特定的本征函数系,才可从中找到个别相同的本征函数。

第二句话反过来说,如果算符 A^\hat{A}B^\hat{B} 不对易,则一般 A^\hat{A}B^\hat{B} 没有共同本征函数,即不能共同测定,这是不确定度关系的表现。

证明(完备共同本征函数系 \Rightarrow 算符对易)

A^ψn=Anψn ,B^ψn=Bnψn \hat{A} \psi_n = A_n \psi_n \ , \kern 1em \hat{B} \psi_n = B_n \psi_n

(A^B^B^A^)ψn=A^(Bnψn)B^(Anψn)=Bn(A^ψn)An(B^ψn)=BnAnψnAnBnψn=0 (\hat{A}\hat{B}-\hat{B}\hat{A}) \psi_n = \hat{A}(B_n\psi_n) - \hat{B}(A_n\psi_n) = B_n(\hat{A}\psi_n) - A_n(\hat{B}\psi_n) = B_nA_n\psi_n - A_nB_n\psi_n = 0

ψ\psi任意波函数,由于 {ψn}{ \psi_n } 完备,则可以将 ψ\psi{ψn}{ \psi_n } 展开

ψ=nCnψn \psi = \sum_n C_n \psi_n

(A^B^B^A^)ψ=(A^B^B^A^)nCnψn=nCn(A^B^B^A^)ψn=0 (\hat{A}\hat{B}-\hat{B}\hat{A}) \psi = (\hat{A}\hat{B}-\hat{B}\hat{A}) \sum_n C_n\psi_n = \sum_n C_n (\hat{A}\hat{B}-\hat{B}\hat{A}) \psi_n = 0

此即说明了

[A^,B^]=0 [\hat{A},\hat{B}] = 0

证明(算符对易 \Rightarrow 存在共同本征函数)

假设算符 A^\hat{A}B^\hat{B} 中至少有一个有非简并的本征值,不妨设 A^\hat{A} 的本征值 λn\lambda_n 是非简并的,即

A^ψn=λnψn \hat{A} \psi_n = \lambda_n \psi_n

A^\hat{A}B^\hat{B} 对易即 A^\hat{A}B^\hat{B} 可交换可得

A^B^ψn=B^A^ψn=B^(λnψn)=λnB^ψn \hat{A} \hat{B} \psi_n = \hat{B} \hat{A} \psi_n = \hat{B} (\lambda_n \psi_n) = \lambda_n \hat{B} \psi_n

A^(B^ψn)=λn(B^ψn)\hat{A}(\hat{B}\psi_n) = \lambda_n (\hat{B}\psi_n) 可知, B^ψn\hat{B}\psi_nA^\hat{A} 对应于本征值 λn\lambda_n 的本征函数,而A^\hat{A} 的本征值 λn\lambda_n 是非简并的,其本征函数只有 ψn\psi_n (乘某一常数),故

B^ψn=μnψn \hat{B}\psi_n = \mu_n \psi_n

ψn\psi_n 也是 B^\hat{B} 的本征函数,故 ψn\psi_nA^\hat{A}B^\hat{B} 的共同本征函数。

而若算符 A^\hat{A}B^\hat{B} 中所有的本征值都是简并的,不妨设算符 A^\hat{A} 对应于本征值 λn\lambda_n 的简并度为 fnf_n ,本征函数为 ψn1,ψn2,,ψnfn\psi_{n1},\psi_{n2},\cdots,\psi_{nf_n} ,并且根据前面对“简并态之间的正交性”的讨论,可以选取这些本征函数为正交归一完备的。由 A^ψnα=λnψnα(α=1,2,,fn)\hat{A} \psi_{n\alpha} = \lambda_n \psi_{n\alpha} \kern 1em (\alpha=1,2,\cdots,f_n) 以及 A^\hat{A}B^\hat{B} 可交换,得

A^B^ψnα=B^A^ψnα=B^(λnψnα)=λnB^ψnα \hat{A} \hat{B} \psi_{n\alpha} = \hat{B} \hat{A} \psi_{n\alpha} = \hat{B} (\lambda_n\psi_{n\alpha}) = \lambda_n \hat{B} \psi_{n\alpha}

A^(B^ψnα)=λn(B^ψnα)\hat{A}(\hat{B}\psi_{n\alpha}) = \lambda_n (\hat{B}\psi_{n\alpha}) ,并不能得到 ψnα\psi_{n\alpha}B^\hat{B} 的本征函数,而只能得到 B^ψnα\hat{B}\psi_{n\alpha} 可以被 {ψn1,ψn2,,ψnfn}{ \psi_{n1},\psi_{n2},\cdots,\psi_{nf_n} } 线性表出,即

B^ψnα=β=1fnμαβψnβ \hat{B} \psi_{n\alpha} = \sum_{\beta=1}^{f_n} \mu_{\alpha\beta} \psi_{n\beta}

(根据“在完备函数系下的展开”)其中 μαβ=(ψnβ,B^ψnα)\mu_{\alpha\beta} = (\psi_{n\beta},\hat{B}\psi_{n\alpha}) 。为了找到一个函数 ϕ\phi 是算符 A^\hat{A}B^\hat{B} 的共同本征函数,令

ϕ=α=1fnCαψnα \phi = \sum_{\alpha=1}^{f_n} C_\alpha \psi_{n\alpha}

A^ϕ=A^(α=1fnCαψnα)=α=1fnCαA^ψnα=α=1fnCαλnψnα=λn(α=1fnCαψnα)=λnϕ \hat{A} \phi = \hat{A} \left( \sum_{\alpha=1}^{f_n} C_\alpha \psi_{n\alpha} \right) = \sum_{\alpha=1}^{f_n} C_\alpha \hat{A} \psi_{n\alpha} = \sum_{\alpha=1}^{f_n} C_\alpha \lambda_n \psi_{n\alpha} = \lambda_n \left( \sum_{\alpha=1}^{f_n} C_\alpha \psi_{n\alpha} \right) = \lambda_n \phi

ϕ\phiA^\hat{A} 的本征函数,而

B^ϕ=B^(α=1fnCαψnα)=α=1fnCαB^ψnα=α=1fnCα(β=1fnμαβψnβ)=β=1fn(α=1fnCαμαβ)ψnβ \hat{B} \phi = \hat{B} \left( \sum_{\alpha=1}^{f_n} C_\alpha \psi_{n\alpha} \right) = \sum_{\alpha=1}^{f_n} C_\alpha \hat{B} \psi_{n\alpha} = \sum_{\alpha=1}^{f_n} C_\alpha \left( \sum_{\beta=1}^{f_n} \mu_{\alpha\beta} \psi_{n\beta} \right) = \sum_{\beta=1}^{f_n} \left( \sum_{\alpha=1}^{f_n} C_\alpha \mu_{\alpha\beta} \right) \psi_{n\beta}

如果想要 ϕ\phi 也是 B^\hat{B} 的本征函数,则需

B^ϕ=μϕ=μα=1fnCαψnα=β=1fnCβμ ψnβ \hat{B} \phi = \mu \phi = \mu \sum_{\alpha=1}^{f_n} C_\alpha \psi_{n\alpha} = \sum_{\beta=1}^{f_n} C_\beta \mu \ \psi_{n\beta}

将两个式子对比,使得每一个 ψnβ\psi_{n\beta} 前的系数都相同,即

α=1fnCαμαβ=Cβμα=1fn(μαβμδαβ)Cα=0 \sum_{\alpha=1}^{f_n} C_\alpha \mu_{\alpha\beta} = C_\beta \mu \Longrightarrow \sum_{\alpha=1}^{f_n} \left( \mu_{\alpha\beta} - \mu \delta_{\alpha\beta} \right) C_\alpha = 0

对于 fnf_nβ\beta ,就构成了一个关于 Cα (α=1,2,,fn)C_\alpha \ (\alpha=1,2,\cdots,f_n) 的线性方程组,为了使得这个线性方程组由非平凡解(即 CαC_\alpha 不全为零),则需系数矩阵不可逆,即

det([μαβμδαβ]fn×fn)=0 \det\left( \begin{bmatrix} \mu_{\alpha\beta} - \mu \delta_{\alpha\beta} \end{bmatrix} _{f_n \times f_n}\right) = 0

这是一个关于 μ\mufnf_n 次方程,可解得 μ\mufnf_n 重根,任取其中一个根,即可解得一组不全为零的 CαC_\alpha ,也就找到了一个不为零的 ϕ\phi ,其是算符 A^\hat{A}B^\hat{B} 的共同本征函数。

事实上,对于上述行列式,可以化为 det([μαβ]fn×fnμIfn)=0\det\left( \begin{bmatrix} \mu_{\alpha\beta} \end{bmatrix}{f_n \times f_n} - \mu I{f_n} \right) = 0 ,不难发现 μ\mu 就是矩阵 [μαβ]fn×fn\begin{bmatrix} \mu_{\alpha\beta} \end{bmatrix}{f_n \times f_n} 的特征值,则可以解得的线性无关的 CαC\alpha 的组数即为矩阵 [μαβ]fn×fn\begin{bmatrix} \mu_{\alpha\beta} \end{bmatrix}_{f_n \times f_n} 特征子空间的维数。

算符与其函数的共同本征态

通过算符函数的定义(展开式),不难发现算符 A^\hat{A} 与其函数 F(A^)F(\hat{A}) 对易,则二者拥有共同本征态:若算符 A^\hat{A} 的本征值问题的解为 A^ψλ=λψλ\hat{A} \psi_{\lambda} = \lambda \psi_{\lambda} ,则算符 F(A^)F(\hat{A}) 的本征值问题的解为

F(A^)ψλ=F(λ)ψλ F(\hat{A}) \psi_{\lambda} = F(\lambda) \psi_{\lambda}

即本征值为 F(λ)F(\lambda) ,本征函数仍为 ψλ\psi_{\lambda}

对易力学量完全集

对易力学量完全集的定义

设有一组彼此独立而且互相对易的厄米算符 A^1,A^2,\hat{A}_1,\hat{A}2,\cdots ,它们的共同本征态记为 ψα\psi\alpha ,其中 α\alpha 表示一组完备的量子数。设给定一组量子数 α\alpha 之后,就能够确定体系的唯一一个可能状态,则称 A^1,A^2,\hat{A}_1,\hat{A}_2,\cdots 构成体系的一组对易力学量完全集(complete set of commuting observables,简记为CSCO),也称为对易可观测量完全集,或简称为力学量完全集。也可以说:力学量完全集指的是互相对易的能够对一个量子体系全部状态进行彻底(不出现简并)地分类标记的最少数目的力学量算符。

对于此定义,可以这么理解:对于本征方程 A^ψ=λψ\hat{A} \psi = \lambda \psi ,当本征值 λ\lambda 一定时,如果 λ\lambda 对应的本征函数不简并,则由 λ\lambda 就唯一确定了当前的量子态 ψ\psi ,此时 A^\hat{A} 自身就构成了力学量完全集;而如果 λ\lambda 对应的本征函数有简并,则需要更多的本征方程加以限制,找到需要最少数目(这与彼此独立是等价的)的算符的一种限制方式,使得对于这组算符任意的一组本征值,共同本征函数只有一个,从而就可以通过这组本征值(即量子数)来唯一确定当前的量子态 ψ\psi 。譬如原子核的能级,只有 ν,l,j\nu,l,j 三个量子数都确定,才可确定一个能级,这三个量子数对应的算符就构成了力学量完全集.

对易力学量完全集的例子
坐标的力学量完全集

算符为

{x^,y^,z^} { \hat{x},\hat{y},\hat{z} }

共同本征函数系为

ψx0,y0,z0(x,y,z)=δ(xx0)δ(yy0)δ(zz0) \psi_{x_0,y_0,z_0} (x,y,z) = \delta(x-x_0) \delta(y-y_0) \delta(z-z_0)

相应的本征值为

{x0,y0,z0} { x_0,y_0,z_0 }

动量的力学量完全集

算符为

{p^x,p^y,p^z} { \hat{p}_x,\hat{p}_y,\hat{p}_z }

共同本征函数系为

ψpx,py,pz(x,y,z)=1(2π)32exp[i(pxx+pyy+pzz)] \psi_{p_x,p_y,p_z} (x,y,z) = \frac{1}{(2\pi\hbar)^{\frac32}} \exp \left[ \frac{\mathrm{i}}{\hbar} (p_xx+p_yy+p_zz) \right]

相应的本征值为

{px,py,pz} { p_x,p_y,p_z }

转动的力学量完全集

算符为

{L^2,L^z} { \hat{L}^2,\hat{L}_z }

在球坐标系下,算符 L^2\hat{L}^2L^z\hat{L}_z 的正交归一的共同本征函数表示为

Ylm(θ,φ)=(1)m2l+14π(lm)!(l+m)! Plm(cosθ) eimφ \mathrm{Y}_{lm} (\theta,\varphi) = (-1)^m \sqrt{\frac{2l+1}{4\pi} \frac{(l-m)!}{(l+m)!}}\ \mathrm{P}_l^m (\cos\theta)\ \mathrm{e}^{\mathrm{i}m\varphi}

式中 YlmY_{lm} 称为球谐函数PlmP_l^m连带Legendre多项式

算符 L^2\hat{L}^2L^z\hat{L}_z 的本征值都是量子化的, ll 称为轨道角动量量子数,可以取 l=0,1,2,3,l=0,1,2,3,\cdots ,分别为 s,p,d,f,s,p,d,f,\cdots 态, mm 称为磁量子数,可以取 m=l,l1,,l+1,lm=l,l-1,\cdots,-l+1,-lL^z\hat{L}_z 的本征值为 mm\hbarL^2\hat{L}^2 的本征值为 l(l+1)2l(l+1)\hbar^2

只考虑算符 L^2\hat{L}^2 时,对于给定的量子数 ll ,存在 (2l+1)(2l+1) 个简并态,故需要算符 L^z\hat{L}_z 对应的量子数 mm 进一步补充,用来区分这些简并态,从而 L^2\hat{L}^2L^z\hat{L}_z 构成力学量完全集。

角动量的共同本征态的求解

L^z\hat{L}_z 的本征值与本征函数

球坐标下角动量 zz 分量算符 L^z=iφ\hat{L}_z = -\mathrm{i}\hbar\frac{\partial}{\partial\varphi} ,本征方程为

iφψ=Lzψ -\mathrm{i}\hbar\frac{\partial}{\partial\varphi} \psi = L_z \psi

解得

ψ(φ)=CeiLzφ \psi(\varphi) = C \mathrm{e}^{\frac{\mathrm{i}}{\hbar}L_z\varphi}

为保证算符 L^z\hat{L}_z 的厄米性,要求波函数 ψ\psi 满足周期性边界条件

ψ(φ+2π)=ψ(φ) \psi(\varphi + 2\pi) = \psi(\varphi)

由此条件可得

CeiLz(φ+2π)=CeiLz(φ)ei2πLz=12πLz=2mπLz=m C \mathrm{e}^{\frac{\mathrm{i}}{\hbar}L_z(\varphi+2\pi)} = C \mathrm{e}^{\frac{\mathrm{i}}{\hbar}L_z(\varphi)} \Longrightarrow \mathrm{e}^{\frac{\mathrm{i}}{\hbar}2\pi L_z} = 1 \Longrightarrow \frac{2\pi L_z}{\hbar} = 2m\pi \Longrightarrow L_z = m\hbar

式中 mm 为整数,由此可知 L^z\hat{L}_z 的本征值是离散的,为

Lz=m(m=0,±1,±2,) L_z = m\hbar \kern 2em (m=0,\pm1,\pm2,\cdots)

相应的本征函数表示为

ψm(φ)=Ceimφ \psi_m(\varphi) = C \mathrm{e}^{\mathrm{i}m\varphi}

按照归一化条件

(ψm,ψm)=02πψm(φ)2 dφ=2πC2=1 (\psi_m,\psi_m) = \int_{0}^{2\pi} |\psi_m(\varphi)|^2\ \mathrm{d}\varphi = 2\pi|C|^2 = 1

故可取 C=12πC=\frac{1}{\sqrt{2\pi}} ,于是归一化本征函数表示为

ψm(φ)=12πeimφ(m=0,±1,±2,) \psi_m(\varphi) = \frac{1}{\sqrt{2\pi}} \mathrm{e}^{\mathrm{i}m\varphi} \kern 2em (m=0,\pm1,\pm2,\cdots)

容易证明这个本征函数系满足正交归一条件

(ψm,ψn)=δmn (\psi_m,\psi_n) = \delta_{mn}

L^2\hat{L}^2 的本征值与本征函数

球坐标下角动量平方算符

L^2=2[1sinθθ(sinθθ)+1sin2θ2φ2] =21sinθθ(sinθθ)+1sin2θL^z2 \hat{L}^2 = -\hbar^2 \left[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\varphi^2} \right] \ \ \ = -\hbar^2 \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta} \hat{L}_z^2

本征方程为

L^2Y(θ,φ)=λ2Y(θ,φ) \hat{L}^2 Y(\theta,\varphi) = \lambda \hbar^2Y(\theta,\varphi)

这里使用 λ2\lambda\hbar^2 作为本征值是为了使得 λ\lambda 无量纲

考虑到 Y(θ,φ)Y(\theta,\varphi) 应为算符 L^2\hat{L}^2L^z\hat{L}_z共同本征函数,则其与 φ\varphi 有关的部分应该与 ψm(φ)\psi_m(\varphi) 相同,而 ψm(φ)\psi_m(\varphi) 的常数部分由 θ\theta 的函数代替,即设

Y(θ,φ)=Θ(θ)ψm(φ) Y(\theta,\varphi) = \Theta(\theta) \psi_m(\varphi)

可以验证 Y(θ,φ)Y(\theta,\varphi) 仍为 L^z\hat{L}_z 的本征函数

L^zY(θ,φ)=L^z[Θ(θ)ψm(φ)]=Θ(θ)[L^zψm(φ)]=Θ(θ)[Lzψm(φ)]=LzΘ(θ)ψm(φ)=LzY(θ,φ) \hat{L}_z Y(\theta,\varphi) = \hat{L}_z \left[ \Theta(\theta) \psi_m(\varphi) \right] = \Theta(\theta) \left[ \hat{L}_z \psi_m(\varphi) \right] = \Theta(\theta) \left[ L_z \psi_m(\varphi) \right] = L_z \Theta(\theta) \psi_m(\varphi) = L_z Y(\theta,\varphi)

Y(θ,φ)Y(\theta,\varphi) 的分离变量式代入算符 L^2\hat{L}^2 的本征方程,整理得

1sinθddθ(sinθddθΘ)+(λm2sin2θ)Θ=0 \frac{1}{\sin\theta} \frac{\mathrm{d}}{\mathrm{d}\theta} \left( \sin\theta \frac{\mathrm{d}}{\mathrm{d}\theta} \Theta \right) + \left( \lambda - \frac{m^2}{\sin^2\theta} \right) \Theta = 0

其中 0θπ0 \le \theta \le \pi ,设 ξ=cosθ (ξ1)\xi = \cos\theta\ (|\xi|\le1) ,代入上述方程得

ddξ[(1ξ2)ddξΘ]+(λm21ξ2)Θ=0 \frac{\mathrm{d}}{\mathrm{d}\xi} \left[ (1-\xi^2) \frac{\mathrm{d}}{\mathrm{d}\xi} \Theta \right] + \left( \lambda - \frac{m^2}{1-\xi^2} \right) \Theta = 0

整理得

(1ξ2)d2dξ2Θ2ξddξΘ+(λm21ξ2)Θ=0 (1-\xi^2) \frac{\mathrm{d}^2}{\mathrm{d}\xi^2} \Theta - 2\xi \frac{\mathrm{d}}{\mathrm{d}\xi} \Theta + \left( \lambda - \frac{m^2}{1-\xi^2} \right) \Theta = 0

此方程为连带Legendre方程,其求解过程如下:

以下求解过程较为复杂,可以选择性阅读。

先考虑Legendre方程,即 m=0m=0 的情形

(1x2)d2dx2y2xddxy+λy=0 (1-x^2) \frac{\mathrm{d}^2}{\mathrm{d}x^2} y - 2x \frac{\mathrm{d}}{\mathrm{d}x} y + \lambda y = 0

此方程可以通过级数解法求解:在 x=0x=0 附近,用幂级数展开

y(x)=k=0+ckxk y(x) = \sum_{k=0}^{+\infty} c_k x^k

代入Legendre方程,比较同幂项的系数,可得

ck+2=k(k+1)λ(k+2)(k+1)ck(k=0,1,2,) c_{k+2} = \frac{k(k+1)-\lambda}{(k+2)(k+1)} c_k \kern 2em (k=0,1,2,\cdots)

故所有的偶次项系数都可以用 c0c_0 来表示,所有的奇次项系数都可以用 c1c_1 来表示,把 c0c_0c1c_1 作为两个任意常数,就可以得到Legendre方程两个线性无关的解,即级数的偶次项部分与奇次项部分

y1(x)=n=0+c2nx2n=c0+c2x2+c4x4+y2(x)=n=0+c2n+1x2n+1=c1x+c3x3+c5x5+ y_1(x) = \sum_{n=0}^{+\infty} c_{2n} x^{2n} = c_0 + c_2x^2 + c_4x^4 + \cdots \ \kern 1em \ y_2(x) = \sum_{n=0}^{+\infty} c_{2n+1} x^{2n+1} = c_1x + c_3x^3 + c_5x^5 + \cdots

考虑当 x±1x\to\pm1 时的情况,当 k+k\to+\infty 时,

ck+2ck=k(k+1)λ(k+2)(k+1)kk+2=12k+2 \frac{c_{k+2}}{c_k} = \frac{k(k+1)-\lambda}{(k+2)(k+1)} \to \frac{k}{k+2} = 1-\frac{2}{k+2}

对于偶数的情况,即 k=2nk=2n ,有 c2n+2/c2n11/(n+1)c_{2n+2}/c_{2n} \sim 1-1/(n+1) ,这与 ln(1x2)\ln(1-x^2) 的Taylor展开

ln(1x2)=n=0+x2nn \ln(1-x^2) = \sum_{n=0}^{+\infty} \frac{x^{2n}}{n}

相邻两项的系数比相同,因此,

y1(x)ln(1x2) y_1(x) \sim \ln(1-x^2)

同理可得

y2(x)xln(1x2) y_2(x) \sim x\ln(1-x^2)

然而,当 x±1x\to\pm1 时,

y1(x),y2(x) y_1(x) \to \infty , \kern 1em y_2(x) \to \infty

这不是物理上可以接受的解,故 y1y_1y2y_2 两个无穷级数解中,必须至少有一个中断为多项式,也就是要找到合适的 λ\lambda ,使得存在 kNk\in\mathbb{N} 满足 k(k+1)λ(k+2)(k+1)=0\frac{k(k+1)-\lambda}{(k+2)(k+1)} = 0 ,故当

λ=l(l+1)(l=0,1,2,) \lambda = l(l+1) \kern 2em (l=0,1,2,\cdots)

时,级数将中断一个多项式( cl+2=cl+4=cl+6==0c_{l+2} = c_{l+4} = c_{l+6} = \cdots = 0 )。当 ll 为偶时, y1y_1 中断为Legendre多项式 Pl(x)\mathrm{P}_l(x)y2y_2 仍为无穷级数;当 nn 为奇时, y2y_2 中断为Legendre多项式 Pl(x)\mathrm{P}_l(x)y1y_1 仍为无穷级数。其中Legendre多项式表示为

Pl(x)=12ll!dldxl(x21)l =k=0l2(1)k(2l2k)!2lk!(lk)!(l2k)!xl2k \mathrm{P}l(x) = \frac{1}{2^l\cdot l!} \frac{\mathrm{d}^l}{\mathrm{d}x^l} (x^2-1)^l \ \ \ = \sum{k=0}^{\left\lfloor \frac{l}{2} \right\rfloor} \frac{(-1)^k(2l-2k)!}{2^l\cdot k!(l-k)!(l-2k)!} x^{l-2k}

例如

P0(x)=1P1(x)=xP2(x)=12(3x21) \mathrm{P}_0(x) = 1 \ \mathrm{P}_1(x) = x \ \mathrm{P}_2(x) = \frac12 (3x^2 - 1)

Legendre多项式的正交性公式表示为

11Pl(x)Pl(x)dx=22l+1δll \int_{-1}^{1} \mathrm{P}l(x) \mathrm{P}{l'}(x) \mathrm{d}x = \frac{2}{2l+1} \delta_{ll'}

回到连带Legendre方程,先考虑在正则奇点 x=1x=1 邻域的行为,令 z=1xz=1-x ,则连带Legendre方程表示为

z(2z)d2dz2y+2(1z)ddzy+[λm2z(2z)]y=0d2dz2y+2(1z)z(2z)ddzy+[λz(2z)m2z2(2z)2]y=0 z(2-z) \frac{\mathrm{d}^2}{\mathrm{d}z^2} y + 2(1-z) \frac{\mathrm{d}}{\mathrm{d}z} y + \left[ \lambda - \frac{m^2}{z(2-z)} \right] y = 0 \ \Downarrow \ \frac{\mathrm{d}^2}{\mathrm{d}z^2} y + \frac{2(1-z)}{z(2-z)} \frac{\mathrm{d}}{\mathrm{d}z} y + \left[ \frac{\lambda}{z(2-z)} - \frac{m^2}{z^2(2-z)^2} \right] y = 0

z=0(x=1)z=0(x=1) 的邻域, 2(1z)2z1\frac{2(1-z)}{2-z} \sim 1λz(2z)λ2z\frac{\lambda}{z(2-z)} \sim \frac{\lambda}{2z}m2z2(2z)2m24z2\frac{m^2}{z^2(2-z)^2} \sim \frac{m^2}{4z^2} 的无穷大阶数要小,则上述方程可化为

d2dz2y+1zddzym24z2y=0z2d2ydz2+zdydzm24y=0 \frac{\mathrm{d}^2}{\mathrm{d}z^2} y + \frac{1}{z} \frac{\mathrm{d}}{\mathrm{d}z} y - \frac{m^2}{4z^2} y = 0 \ \Downarrow \ z^2\frac{\mathrm{d}^2y}{\mathrm{d}z^2} +z \frac{\mathrm{d}y}{\mathrm{d}z} - \frac{m^2}{4} y = 0

此方程为欧拉方程,解具有 zsz^s 的形式,代入可得

s(s1)+sm24=0 s(s-1) + s - \frac{m^2}{4} = 0

解得 s=±m2s=\pm\frac{|m|}{2} ,但在 z=0z=0 邻域,解 yzm/2y \propto z^{-|m|/2} \to \infty 不满足物理上的要求,因此在 z=0(x=1)z=0(x=1) 邻域,有

yzm2=(1x)m2 y \propto z^{\frac{|m|}{2}} = (1-x)^{\frac{|m|}{2}}

同理,在 x=1x=-1 邻域,有

y(1+x)m2 y \propto (1+x)^{\frac{|m|}{2}}

故可令连带Legendre方程的解为

y(x)=(1x)m2(1+x)m2v(x)=(1x2)m2v(x) y(x) = (1-x)^{\frac{|m|}{2}}(1+x)^{\frac{|m|}{2}}v(x) = (1-x^2)^{\frac{|m|}{2}}v(x)

代入连带Legendre方程

(1x2)d2ydx22xdydx+(λm21x2)y=0 (1-x^2) \frac{\mathrm{d}^2y}{\mathrm{d}x^2} - 2x \frac{\mathrm{d}y}{\mathrm{d}x} + \left( \lambda - \frac{m^2}{1-x^2} \right) y = 0

可得

(1x2)v2(m+1)xv+[λm(m+1)]v=0 (1-x^2)v'' - 2(|m|+1)xv' + [\lambda-|m|(|m|+1)]v = 0

该方程对 xx 求导,可得

(1x2)v2(m+2)xv+[λ(m+1)(m+2)]v=0 (1-x^2)v''' - 2(|m|+2)xv'' + [\lambda-(|m|+1)(|m|+2)]v' = 0

可推得每经过一次求导,发生如下变换: mm+1, vv|m|\to|m|+1,\ v\to v' 。而当 m=0m=0 时,微分方程 (1x2)v2(m+1)xv+[λm(m+1)]v=0(1-x^2)v'' - 2(|m|+1)xv' + [\lambda-|m|(|m|+1)]v = 0 为Legendre方程,而要把 m|m| 变为 00 相当于经过了 m|m| 反求导,则 vv 反求导 m|m| 次可以得到Legendre方程的解 Pl(x)\mathrm{P}_l(x) ,故

v(x)=dmdxmPl(x) v(x) = \frac{\mathrm{d}^{|m|}}{\mathrm{d}x^{|m|}} \mathrm{P}_l(x)

这样就得到了连带Legendre方程物理上允许的解为连带Legendre多项式

Plm(x)=(1x2)m2dmdxmPl(x) \mathrm{P}_{l}^{|m|} (x) = (1-x^2)^{\frac{|m|}{2}} \frac{\mathrm{d}^{|m|}}{\mathrm{d}x^{|m|}} \mathrm{P}_l(x)

对于 m0m\ge0 的情况,

Plm(x)=(1x2)m2dmdxmPl(x) \mathrm{P}_{l}^{m} (x) = (1-x^2)^{\frac{m}{2}} \frac{\mathrm{d}^{m}}{\mathrm{d}x^{m}} \mathrm{P}_l(x)

Pl(x)\mathrm{P}_l(x) 的表达式代入可得

Plm(x)=12ll!(1x2)m2dldxl(x21)l \mathrm{P}_{l}^{m} (x) = \frac{1}{2^l\cdot l!} (1-x^2)^{\frac{m}{2}} \frac{\mathrm{d}^l}{\mathrm{d}x^l} (x^2-1)^l

该式对于 lm<0-l \le m < 0 也有意义,可以证明

Plm(x)=(1)m(lm)!(l+m)!Plm(x) \mathrm{P}{l}^{-m} (x) = (-1)^m \frac{(l-m)!}{(l+m)!} \mathrm{P}{l}^{m} (x)

连带Legendre多项式的正交性公式表示为

11Plm(x)Plm(x)dx=(l+m)!(lm)!22l+1δll \int_{-1}^{1} \mathrm{P}l^m(x) \mathrm{P}^m{l'}(x) \mathrm{d}x = \frac{(l+m)!}{(l-m)!} \frac{2}{2l+1} \delta_{ll'}

这样就基本完成了算符 L^2\hat{L}^2 的本征函数求解,即球谐函数

Ylm(θ,φ)=NlmPlm(cosθ)ψm(φ) \mathrm{Y}{lm} (\theta,\varphi) = N{lm} \mathrm{P}^m_l(\cos\theta) \psi_m(\varphi)

其中归一化系数

Nlm=(1)m2l+14π(lm)!(l+m)! N_{lm} = (-1)^m \sqrt{\frac{2l+1}{4\pi} \frac{(l-m)!}{(l+m)!}}

Ylm(θ,φ)=(1)m2l+14π(lm)!(l+m)! Plm(cosθ) eimφ \mathrm{Y}_{lm} (\theta,\varphi) = (-1)^m \sqrt{\frac{2l+1}{4\pi} \frac{(l-m)!}{(l+m)!}}\ \mathrm{P}_l^m (\cos\theta)\ \mathrm{e}^{\mathrm{i}m\varphi}

其满足

L^2Ylm=l(l+1)2YlmL^zYlm=mYlm \hat{L}^2 \mathrm{Y}{lm} = l(l+1) \hbar^2 \mathrm{Y}{lm} \ \hat{L}z \mathrm{Y}{lm} = m\hbar \mathrm{Y}_{lm}

其中 l=0,1,2,,m=l,l1,,l+1,ll=0,1,2,\cdots, \kern 1em m=l,l-1,\cdots,-l+1,-l ;球谐函数的正交归一化表示为

02πdφ0πsinθdθYlm(θ,φ)Ylm(θ,φ)=δllδmm \int_{0}^{2\pi} \mathrm{d}\varphi \int_{0}^{\pi} \sin\theta \mathrm{d}\theta \mathrm{Y}{lm}^*(\theta,\varphi) Y{l'm'}(\theta,\varphi) = \delta_{ll'} \delta_{mm'}

头两阶球谐函数为

Y0,0(θ,φ)=14π Y1,0(θ,φ)=34πcosθ Y1,±1(θ,φ)=38πsinθe±iφ \mathrm{Y}{0,0}(\theta,\varphi) = \frac{1}{\sqrt{4\pi}} \ \ \ \mathrm{Y}{1,0}(\theta,\varphi) = \sqrt{\frac{3}{4\pi}} \cos\theta \ \ \ \mathrm{Y}_{1,\pm1}(\theta,\varphi) = \mp\sqrt{\frac{3}{8\pi}} \sin\theta \mathrm{e}^{\pm\mathrm{i}\varphi}

3.4 不确定度关系

不确定度关系的表述、含义及导出

不确定度

定义算符 A^\hat{A} 对应的偏差算符

ΔA^=A^Aˉ \Delta\hat{A} = \hat{A} - \bar{A}

容易发现 ΔA^\Delta\hat{A} 的平均值 ΔA^=A^Aˉ=AˉAˉ=0\overline{\Delta\hat{A}} = \overline{\hat{A} - \bar{A}} = \bar{A} - \bar{A} = 0 ,可以知道这个偏差算符的平均值并没有实际意义,因为偏差的正负被抵消为零了,所以应该考虑其平方的平均值。

每次测量的结果围绕平均值有一个涨落,其定义为

(ΔA^)2=(A^Aˉ)2=A^22A^Aˉ+Aˉ2=A^22AˉAˉ+Aˉ2=A^2Aˉ2 \overline{(\Delta\hat{A})^2} = \overline{(\hat{A} - \bar{A})^2} = \overline{\hat{A}^2 - 2\hat{A}\bar{A} + \bar{A}^2} = \overline{\hat{A}^2} - 2\bar{A}\bar{A} + \bar{A}^2 = \overline{\hat{A}^2} - \bar{A}^2

这个量描写了力学量 AA 的测量值的偏差程度。

在态 ψ\psiAA 的取值的不确定度

ΔA=(ΔA^)2=A^2Aˉ2 \Delta A = \sqrt{\overline{(\Delta\hat{A})^2}} = \sqrt{\overline{\hat{A}^2} - \bar{A}^2}

不确定度关系

在任意量子态 ψ\psi 下任意两个力学量 A,BA,B 的不确定度的乘积存在下限,即

ΔAΔB12[A^,B^]=12(ψ,[A^,B^]ψ) \Delta A \cdot \Delta B \ge \frac12 \left| \overline{[\hat{A},\hat{B}]} \right| = \frac12 \left| (\psi, [\hat{A},\hat{B}]\psi) \right|

[A^,B^]0[\hat{A},\hat{B}] \ne 0 时,除了 [A^,B^]=(ψ,[A^,B^]ψ)=0\overline{[\hat{A},\hat{B}]} = (\psi,[\hat{A},\hat{B}]\psi) = 0 的特殊情况外,在任何态下 A,BA,B 不可能同时取确定值;而当 [A^,B^]=0[\hat{A},\hat{B}] = 0 时,A,BA,B 可同时取确定值。

这里的特殊情况,例如氢原子的基态 ψ100(r,θ,φ)\psi_{100}(r,\theta,\varphi) ,其总角动量 L=0L=0zz 方向角动量 Lz=0L_z=0 ,可推得 Lx=Ly=0L_x = L_y = 0 ,此时 Lx,Ly,LzL_x,L_y,L_z 同时取得确定值;然而 Lx,Ly,LzL_x,L_y,L_z 不对易,这里可以同时取得确定值是因为对易关系的平均值为零,如 [L^x,L^y]=iL^z=0\overline{[\hat{L}_x,\hat{L}_y]} = \overline{\mathrm{i}\hbar\hat{L}_z} = 0

说明:

  • 不确定度关系是微观粒子运动的基本规律,在宏观世界不能得到直接的体现;

  • 不确定度关系是微观粒子固有属性决定的,与仪器精度和测量方法的缺陷无关;

  • 不确定度关系不是给物理学带来了不精确性,而正是体现了微观世界的精确性;

  • 不确定度关系给出了微观世界中应用经典粒子的坐标和动量概念时应受到的限制,展示了量子力学与经典力学规律的本质区别。

常见的不确定度关系
坐标与动量的不确定度关系

波动性使微观粒子没有确定的轨道, 即坐标和动量不能同时取确定值, 它们存在一个不确定度关系

ΔxΔpx2,ΔyΔpy2,ΔzΔpz2 \Delta x \Delta p_x \ge \frac{\hbar}{2} ,\kern 1em \Delta y \Delta p_y \ge \frac{\hbar}{2} ,\kern 1em \Delta z \Delta p_z \ge \frac{\hbar}{2}

能量与时间的不确定度关系

ΔEΔt2 \Delta E \Delta t \ge \frac{\hbar}{2}

时间实际上不是力学量,这个不确定度关系不能直接用上述方法得到。

能量与时间的不确定关系式说明了原子的激发态能级都有一定的能级宽度, 实验原子发光都有一定的频率宽度

角动量分量与方位角的不确定度关系

ΔLzΔφ2 \Delta L_z \Delta \varphi \ge \frac{\hbar}{2}

不确定度关系的证明

A^,B^\hat{A},\hat{B} 均为厄米算符、 Aˉ,Bˉ\bar{A},\bar{B} 均为实数,可推得 ΔA^=A^Aˉ, ΔB^=B^Bˉ\Delta\hat{A}=\hat{A}-\bar{A},\ \Delta\hat{B}=\hat{B}-\bar{B} 也均为厄米算符。考虑积分不等式

I(ξ)=(ξΔA^+iB^)ψ2dτ0 I(\xi) = \int \left| (\xi\Delta\hat{A} + \mathrm{i}\hat{B}) \psi \right|^2 \mathrm{d}\tau \ge 0

其中 ξ\xi 为任意实参数,

I(ξ)=(ξΔA^ψ+iB^ψ)(ξΔA^ψ+iB^ψ)dτ =[ξ(ΔA^ψ)i(B^ψ)][ξ(ΔA^ψ)+i(B^ψ)]dτ =ξ2(ΔA^ψ)(ΔA^ψ)dτ+iξ(ΔA^ψ)(ΔB^ψ)dτiξ(ΔB^ψ)(ΔA^ψ)dτ+(ΔB^ψ)(ΔB^ψ)dτ =ξ2(ΔA^ψ,ΔA^ψ)+iξ(ΔA^ψ,ΔB^ψ)iξ(ΔB^ψ,ΔA^ψ)+(ΔB^ψ,ΔB^ψ) =ξ2(ψ,(ΔA^)2ψ)+iξ(ψ,ΔA^ΔB^ψ)iξ(ψ,ΔB^ΔA^ψ)+(ψ,(ΔB^)2ψ) =ξ2(ΔA^)2+iξ(ΔA^ΔB^ΔB^ΔA^)+(ΔB^)2 =ξ2(ΔA^)2+iξ[ΔA^,ΔB^]+(ΔB^)2 I(\xi) = \int \left( \xi\Delta\hat{A}\psi + \mathrm{i}\hat{B}\psi \right)^* \left( \xi\Delta\hat{A}\psi+\mathrm{i}\hat{B}\psi \right) \mathrm{d}\tau \ \ \ = \int \left[ \xi\left(\Delta\hat{A}\psi\right)^* - \mathrm{i}\left(\hat{B}\psi\right)^* \right] \left[ \xi\left(\Delta\hat{A}\psi\right) + \mathrm{i}\left(\hat{B}\psi\right) \right] \mathrm{d}\tau \ \ \ = \xi^2 \int \left(\Delta\hat{A}\psi\right)^* \left(\Delta\hat{A}\psi\right) \mathrm{d}\tau + \mathrm{i}\xi \int \left(\Delta\hat{A}\psi\right)^* \left(\Delta\hat{B}\psi\right) \mathrm{d}\tau - \mathrm{i}\xi \int \left(\Delta\hat{B}\psi\right)^* \left(\Delta\hat{A}\psi\right) \mathrm{d}\tau + \int \left(\Delta\hat{B}\psi\right)^* \left(\Delta\hat{B}\psi\right) \mathrm{d}\tau \ \ \ = \xi^2 \left(\Delta\hat{A}\psi , \Delta\hat{A}\psi\right) + \mathrm{i}\xi \left(\Delta\hat{A}\psi , \Delta\hat{B}\psi\right) - \mathrm{i}\xi \left(\Delta\hat{B}\psi , \Delta\hat{A}\psi\right) + \left(\Delta\hat{B}\psi , \Delta\hat{B}\psi\right) \ \ \ = \xi^2 \left(\psi , (\Delta\hat{A})^2\psi\right) + \mathrm{i}\xi \left(\psi , \Delta\hat{A}\Delta\hat{B}\psi\right) - \mathrm{i}\xi \left(\psi , \Delta\hat{B}\Delta\hat{A}\psi\right) + \left(\psi , (\Delta\hat{B})^2\psi\right) \ \ \ = \xi^2 \overline{(\Delta\hat{A})^2} + \mathrm{i}\xi \left( \overline{\Delta\hat{A}\Delta\hat{B}-\Delta\hat{B}\Delta\hat{A}}\right) + \overline{(\Delta\hat{B})^2} \ \ \ = \xi^2 \overline{(\Delta\hat{A})^2} + \mathrm{i}\xi \overline{[\Delta\hat{A} , \Delta\hat{B}]} + \overline{(\Delta\hat{B})^2}

考虑 [A^,B^][\hat{A},\hat{B}][ΔA^,ΔB^][\Delta\hat{A},\Delta\hat{B}] 之间的关系

[ΔA^,ΔB^]=[A^Aˉ,B^Bˉ]=[A^,B^] [\Delta\hat{A},\Delta\hat{B}] = [\hat{A}-\bar{A},\hat{B}-\bar{B}] = [\hat{A},\hat{B}]

由于上述连等式中所有等号连接的均为实数,引入厄米算符 C^=[A^,B^] / i=C^+\hat{C} = [\hat{A},\hat{B}]\ /\ \mathrm{i} = \hat{C}^+

I(ξ)=ξ2(ΔA^)2+iξ[A^,B^]+(ΔB^)2=(ΔA^)2ξ2C^ξ+(ΔB^)2 =(ΔA^)2(ξC^2(ΔA^)2)2+((ΔB^)2C^24(ΔA^)2)0 I(\xi) = \xi^2 \overline{(\Delta\hat{A})^2} + \mathrm{i}\xi \overline{[\hat{A} , \hat{B}]} + \overline{(\Delta\hat{B})^2} = \overline{(\Delta\hat{A})^2} \xi^2 - \overline{\hat{C}} \xi + \overline{(\Delta\hat{B})^2} \ \ \ = \overline{(\Delta\hat{A})^2} \left( \xi - \frac{\overline{\hat{C}}}{2\overline{(\Delta\hat{A})^2}} \right)^2 + \left( \overline{(\Delta\hat{B})^2} - \frac{\overline{\hat{C}}^2}{4\overline{(\Delta\hat{A})^2}} \right) \ge 0

由于 ξ\xi 为任意实参数,故需要 I(ξ)I(\xi) 的最小值非负,即

(ΔB^)2C^24(ΔA^)20(ΔA^)2(ΔB^)2C^240ΔAΔBC^2=12[A^,B^] \overline{(\Delta\hat{B})^2} -\frac{\overline{\hat{C}}^2}{4\overline{(\Delta\hat{A})^2}} \ge 0 \Longrightarrow \overline{(\Delta\hat{A})^2} \cdot \overline{(\Delta\hat{B})^2} -\frac{\overline{\hat{C}}^2}{4} \ge 0 \Longrightarrow \Delta A \Delta B \ge \frac{\left|\overline{\hat{C}}\right|}{2} = \frac12 \left| \overline{[\hat{A},\hat{B}]} \right|

不确定度关系应用举例

一维谐振子的零点能

一维谐振子的Hamiton算符

H^=12mp^2+12mω2x^2 \hat{H} = \frac{1}{2m} \hat{p}^2 + \frac12 m\omega^2 \hat{x}^2

能量本征函数为

ψn(x)=Aneα2x22Hn(αx) \psi_n(x) = A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x)

其中 α=mω\alpha = \sqrt{\frac{m\omega}{\hbar}} , Hn\mathrm{H}_n 为Hermite多项式, AnA_n 为归一化系数(取正实数)。考虑 xxpp 的平均值

xˉ=(ψ,x^ψ)=+Aneα2x22Hn(αx) x Aneα2x22Hn(αx) dx =An2+eα2x2 Hn2(αx) x dx=0 \bar{x} = (\psi , \hat{x}\psi) = \int_{-\infty}^{+\infty} A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x) \ x \ A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}n(\alpha x) \ \mathrm{d} x \ \ \ = A_n^2 \int{-\infty}^{+\infty} \mathrm{e}^{-\alpha^2x^2} \ \mathrm{H}_n^2(\alpha x) \ x \ \mathrm{d}x = 0

最后一个等号是因为被积函数是奇函数( Hn(x)=(1)nHn(x)\mathrm{H}_n(-x) = (-1)^n \mathrm{H}_n(x) )。

pˉ=(ψ,p^ψ)=+Aneα2x22Hn(αx) iddx[Aneα2x22Hn(αx)]dx =iAn2+eα2x22Hn(αx) ddx[eα2x22Hn(αx)]dx =iAn2eα2x2 Hn2(αx)++ iAn2+eα2x22Hn(αx) ddx[eα2x22Hn(αx)]dx =iAn2+eα2x22Hn(αx) ddx[eα2x22Hn(αx)]dx \bar{p} = (\psi , \hat{p}\psi) = \int_{-\infty}^{+\infty} A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x) \ \frac{\hbar}{\mathrm{i}}\frac{\mathrm{d}}{\mathrm{d}x} \left[ A_n \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}n(\alpha x) \right] \mathrm{d} x \ \ \ = -\mathrm{i}\hbar A_n^2 \int{-\infty}^{+\infty} \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x) \ \frac{\mathrm{d}}{\mathrm{d}x} \left[ \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}n(\alpha x) \right] \mathrm{d} x \ \ \ = \left. -\mathrm{i}\hbar A_n^2 \mathrm{e}^{-\alpha^2x^2} \ \mathrm{H}n^2(\alpha x) \right|{-\infty}^{+\infty} +\ \mathrm{i}\hbar A_n^2 \int{-\infty}^{+\infty} \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x) \ \frac{\mathrm{d}}{\mathrm{d}x} \left[ \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}n(\alpha x) \right] \mathrm{d} x \ \ \ = \mathrm{i}\hbar A_n^2 \int{-\infty}^{+\infty} \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x) \ \frac{\mathrm{d}}{\mathrm{d}x} \left[ \mathrm{e}^{-\frac{\alpha^2x^2}{2}} \mathrm{H}_n(\alpha x) \right] \mathrm{d} x

第三行的通过分部积分法得到

由第二行与第四行相等可得 pˉ=pˉ\bar{p} = -\bar{p} ,即 pˉ=0\bar{p} = 0

ΔA=A^2Aˉ2\Delta A = \sqrt{\overline{\hat{A}^2} - \bar{A}^2} 可知

(Δx)2=x^2,(Δp)2=p^2 (\Delta x)^2 = \overline{\hat{x}^2} ,\kern 1em (\Delta p)^2 = \overline{\hat{p}^2}

Eˉ=12mp^2+12mω2x^2 =12m(Δp)2+12mω2(Δx)2 212m(Δp)212mω2(Δx)2 =ωΔpΔx 12ω \bar{E} = \frac{1}{2m} \overline{\hat{p}^2} + \frac12 m\omega^2 \overline{\hat{x}^2} \ \ \ = \frac{1}{2m} (\Delta p)^2 + \frac12 m\omega^2 (\Delta x)^2 \ \ \ \ge 2 \sqrt{ \frac{1}{2m} (\Delta p)^2 \cdot \frac12 m\omega^2 (\Delta x)^2 } \ \ \ = \omega \Delta p \Delta x \ \ \ \ge \frac12 \hbar\omega

即得一维谐振子的零点能为

Eˉmin=12ω \bar{E}_{\min} = \frac12 \hbar\omega

非零的零点能是不确定度关系的结果。

球谐函数作为本征态时 LxL_xLyL_y 的不确定度关系

求证:在 LzL_z 的本征态 Ylm\mathrm{Y}_{lm} 中, LxL_xLyL_y 的不确定度关系为

(ΔL^x)2(ΔL^y)2m244 \overline{(\Delta\hat{L}_x)^2} \cdot \overline{(\Delta\hat{L}_y)^2} \ge \frac{m^2\hbar^4}{4}

证明

L^z\hat{L}_z 的本征方程为

L^zYlm=mYlm \hat{L}z \mathrm{Y}{lm} = m\hbar \mathrm{Y}_{lm}

根据对易关系 [L^x,L^y]=iL^z[\hat{L}_x,\hat{L}_y] = \mathrm{i}\hbar \hat{L}_z ,可得

ΔLxΔLy12[L^x,L^y]=2L^z=m22 \Delta L_x \Delta L_y \ge \frac12 \left| \overline{[\hat{L}_x,\hat{L}_y]} \right| = \frac{\hbar}{2} \left|\overline{\hat{L}_z}\right| = \frac{m\hbar^2}{2}

(ΔL^x)2(ΔL^y)2m244 \overline{(\Delta\hat{L}_x)^2} \cdot \overline{(\Delta\hat{L}_y)^2} \ge \frac{m^2\hbar^4}{4}